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1 Preface 

We present mathematical details of several cosmological models, whereby the 
topological and the geometrical background will be emphasized. 

This book arose from lectures I read as advanced courses in the following 
universities: 1990 at TU Berlin, 1991 in Minister, 1993 - 1996 and 2002, 2003 
in Potsdam, 1999 and 2000 in Salerno, and 2001 at FU Berlin. 

The reader is assumed to be acquainted with basic knowledge on general 
relativity, e.g. by knowing Stephani's book [250]. As a rule, the deduction is 
made so explicit, that one should be able to follow the details without further 
background knowledge. The intention is to present a piece of mathematical 
physics related to cosmology to fill the gap between standard textbooks like 
MTW 1 , LL 2 , or Weinberg 3 , on the one hand and to the current research 
literature on the other. I hope that this will become useful for both students 
and researchers who are interested in cosmology, mathematical physics or 
differential geometry. 

To the reference list: If more than one source is given in one item, then a dot 
is put between them. If, however, two versions of the same source are given, 
then a semicolon is put between them. The gr-qc number in the reference 
list refers to the number in the preprint archive 

http://xxx.lanl.gov/find/ 
and the numbers in italics given at the end of an item show the pages where 
this item is mentioned. This way, the Bibliography simultaneously serves as 

1 C. Misner, K. Thorne, J. Wheeler Gravitation, Freeman, San Francisco 1973 

2 L. Landau, E. Lifschitz Field Theory, Nauka, Moscow 1973 

3 S. Weinberg Gravitation and Cosmology, J. Wiley, New York 1972 
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author index. 

Parts of this book 4 are composed from revised versions of papers of mine 
published earlier in research journals. They are reprinted with the kind 
permissions of the copyright owners: Akademie-Verlag Berlin for Astron. 
Nachr. and Plenum Press New York for Gen. Rel. Grav. The details about 
the originals can also be seen in 

http: / /www. physik.fu-berlin.de/~hjschmi 

I thank the following colleagues: Luca Amendola, Michael Bachmann, H. 
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v. Borzeszkowski, Salvatore Capozziello, G. Dautcourt, Vladimir Dzhunu- 
shaliev, Mauro Francaviglia, T. Friedrich, V. Frolov, H. Goenner, Stefan 
Gottlober, V. Gurovich, Graham Hall, F. Hehl, Olaf Heinrich, Alan Held, L. 
Herrera, Rainer John, U. Kasper, Jerzy Kijowski, Hagen Kleinert, B. Klotzek, 
Sabine Kluske, O. Kowalski, Andrzej Krasinski, W. Kiihnel, W. Kundt, K. 
Lake, K. Leichtweiss, D.-E. Liebscher, Malcolm MacCallum, W. Mai, R. 
Mansouri, V. Melnikov, P. Michor, Salvatore Mignemi, Volker Miiller, S. 
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Quandt, L. Querella, M. Rainer, J. Reichert, Stefan Reuter, W. Rinow, 
M. Sanchez, R. Schimming, M. Schulz, E. Schmutzer, U. Semmelmann, I. 
Shapiro, A. Starobinsky, P. Teyssandier, H.-J. Treder, H. Tuschik, L. Van- 
hecke and A. Zhuk for valuable comments, and my wife Renate for patience 
and support. 

Hans-Jiirgen Schmidt, Potsdam, July 2004 



The details are mentioned at the corresponding places. 
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2 Cosmological models with changes of the 
Bianchi type 

We investigate such cosmological models which instead of the usual spatial 
homogeneity property only fulfil the condition that in a certain synchronized 
system of reference all spacelike sections t = const, are homogeneous mani- 
folds. 

This allows time-dependent changes of the Bianchi type. Discussing 
differential-geometrical theorems it is shown which of them are permitted. 
Besides the trivial case of changing into type I there exist some possible 
changes between other types. However, physical reasons like energy inequal- 
ities partially exclude them. 

2.1 Introduction to generalized Bianchi models 

Recently, besides the known Bianchi models, there are investigated certain 
classes of inhomogeneous cosmological models. This is done to get a better 
representation of the really existent inhomogeneities, cf. e.g. Bergmann [23], 
Carmeli [51], Collins [55], Spero [240], Szekeres [251] and Wainwright [265]. 

We consider, similar as in Collins [55], such inhomogeneous models V4 
which in a certain synchronized system of reference possess homogeneous 
sections t = const., called V 3 (t). This is analogous to the generalization of 
the concept of spherical symmetry in Krasinski [136], cf. also [137]. 

In this chapter, which is based on [200], we especially investigate which 
time-dependent changes of the Bianchi type are possible. Thereby we im- 
pose, besides the twice continuous differentiability, a physically reasonable 
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2 CHANGES OF THE BIANCHI TYPE 



condition: the energy inequality, 

Too > \T aP \ , (2.1) 

holds in each Lorentz frame. 

We begin with some globally topological properties of spacetime: Under 
the physical condition (2.1) the topology of the sections V 3 (t) is, according 
to Lee [142], independent of t. Hence, the Kantowski-Sachs models, with 
underlying topology S 2 x R or S 2 x S 1 and the models of Bianchi type IX 
with underlying topology S 3 or continuous images of it as SO (3) may not 
change, because all other types are represented by the R 3 -topology, factorized 
with reference to a discrete subgroup of the group of motions. But the 
remaining types can all be represented in R 3 -topology itself; therefore we do 
not get any further global restrictions. All homogeneous models, the above 
mentioned Kantowski-Sachs model being excluded, possess simply transitive 
groups of motion. If the isometry group G has dimension > 4, then this 
statement means: G possesses a 3-dimensional simply-transitive subgroup 
acting transitively on the V^(t). Hence we do not specialize if we deal only 
with locally simply-transitive groups of motions. We consequently do not 
consider here trivial changes of the Bianchi type, e. g. from type III to type 
VIII by means of an intermediate on which a group of motions possessing 
transitive subgroups of both types acts. 

Next, we consider the easily tractable case of a change to type I: For each 
type M there one can find a manifold V4 such that for each t < the section 
V 3 (t) is flat and for each t > it belongs to type M. Indeed, one has simply 
to use for every V 3 (t) such a representative of type M that their curvature 
vanishes as t — > 0+. Choosing exponentially decreasing curvature one can 
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obtain an arbitrarily high differentiable class for the metric, e. g. 



ds 2 = -dt 2 + dr 2 + h 2 ■ {dip 2 + sin 2 ipd(p 2 ) , 



(2.2) 



where 



r 



for 



t < 



h 



exp(t 2 ) • sinh 2 (r • exp(— t 2 )) else . 



(2.3) 



This is a C^-metric whose slices V%(i) belong to type I for t < and to 
type V for t > 0. In t = 0, of course, it cannot be an analytical one. The 
limiting slice belongs necessarily to type I by continuity reasons. Applying 
this fact twice it becomes obvious that by the help of a flat intermediate of 
finite extension or only by a single flat slice, all Bianchi types can be matched 
together. However, if one does not want to use such a flat intermediate the 
transitions of one Bianchi type to another become a non-trivial problem. It 
is shown from the purely differential-geometrical as well as from the physical 
points of view, sections 2.3 and 2.4 respectively, which types can be matched 
together immediately without a flat intermediate. To this end we collect the 
following preliminaries. 

2.2 Spaces possessing homogeneous slices 

As one knows, in cosmology the homogeneity principle is expressed by the 
fact that to a spacetime V4 there exists a group of motions acting transitively 
on the spacelike hypersurfaces V 3 (t) of a slicing of V4. Then the metric is 
given by 



where the g a b(t) are positively definite and u a are the basic 1-forms corre- 
sponding to a certain Bianchi type; for details see e.g. [134]. If the u a are 



ds 2 = -dt 2 + g ab {t)oo a u 



b 



(2.4) 
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2 CHANGES OF THE BIANCHI TYPE 



related to a holonomic basis x\ using type- dependent functions A^(x^), one 
can write: 

u a = A?dx i . (2.5) 

For the spaces considered here we have however: there is a synchronized 
system of reference such that the slices t = const, are homogeneous spaces 
V 3 (t). In this system of reference the metric is given by 

ds 2 = -dt 2 + g ij dx i dx j (2.6) 

where gij(x\t) are twice continuously differentiate and homogeneous for 
constant t. Hence it is a generalization of eqs. (2.4) and (2.5). This is a 
genuine generalization because there is only posed the condition g 0a = —5 0a 
on the composition of the homogeneous slices V 3 (t) to a V4. This means in 
each slice only the first fundamental form g^ is homogeneous; but in the 
contrary to homogeneous models, the second fundamental form Toafe need 
not have this property. Therefore, also the curvature scalar ^R, and with it 
the distribution of matter, need not be constant within a V 3 (t). 

Now let t be fixed. Then one can find coordinates x l t (x j , t) in V 3 (t) such 
that according to eqs. (2.4) and (2.5) the inner metric gets the form 

g ab (t) u? u b t where cu? = A^x\) dx\ . (2.7) 

Transforming this into the original coordinates x % one obtains for the metric 
of the full V 4 



90a = Soa , 9ij = 9ab(t) A%(x$ A\{x\) x\, X l tJ . 



(2.8) 



2.3 Continuous changes of the Bianchi type 
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If the Bianchi type changes with time one has to take such an ansatz 5 like 
eqs. (2.2), (2.3) for each interval of constant type separately; between them 
one has to secure a C 2 -joining. An example of such an inhomogeneous model 
is, cf. Ellis 1967 [70] 

ds 2 = -dt 2 + r 2 / 3 [t + C(x)} 2 dx 2 + t^{dy 2 + dz 2 ) . (2.9) 

The slices t = const, are flat, but only for constant C it belongs to Bianchi 
type I, cf. also [99] for such a model. Later, see eq. (3.19), we will present a 
more general model. 

2.3 Continuous changes of the Bianchi type 

Using the usual homogeneity property the same group of motions acts on 
each slice Vsit). Hence the Bianchi type is independent of time by definition. 
But this fails to be the case for the spaces considered here. However, we 
can deduce the following: completing d/dt to an anholonomic basis which 
is connected with Killing vectors in each V^it) one obtains for the structure 
constants associated to the commutators of the basis: 

Ca/3 = j Cjfe(0 depends continuously on time . (2-10) 

Note: These structure constants are calculated as follows: Without loss of 
generality let x l t (Q : 0, 0,t) = 0. At this point d/dx\ and djdx 1 are taken as 
initial values for Killing vectors within V 3 (t). The structure constants ob- 
tained by these Killing vectors are denoted by C l jk (t) and C l jk (t) respectively, 
where Cj k (t) are the canonical ones. Between them it holds at (0, 0, 0, t): 

x t,jC 3 kl = C l mn x™ k x™ t . (2.11) 

5 The word "ansatz" (setting) stems from the German noun "der Ansatz"; that word 
stems from the German verb "setzen" (to set). 
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2 CHANGES OF THE BIANCHI TYPE 



Using only the usual canonical structure constants then of course no type 
is changeable continuously To answer the question which Bianchi types may 
change continuously we consider all sets of structure constants Cj fc being 
antisymmetric in jk and fulfilling the Jacobi identity. Let Cs be such a 
set belonging to type S. Then we have according to eq. (2.10): Type R is 
changeable continuously into type S, symbolically expressed by the validity 
of R — > S, if and only if to each C s there exists a sequence such that 
one has 



The limit has to be understood componentwise. 

It holds that R — ■> S if and only if there are a Cs and a sequence such 
that eq. (2.12) is fulfilled. The equivalence of both statements is shown by 
means of simultaneous rotations of the basis. In practice one takes as C s the 
canonical structure constants and investigates for which types R there can 
be found corresponding sequences C^: the components of Cs are subjected 
to a perturbation not exceeding e and their Bianchi types are calculated. 
Finally one looks which types appear for all e > 0. Thereby one profits, e.g., 
from the statement that the dimension of the image space of the Lie algebra, 
which equals for type I, . . . , and equals 3 for types VIII and IX, cannot 
increase during such changes; see [189] and [210] for more details about these 
Lie algebras. 

One obtains the following diagram. The validity of R — > S and S — > 
T implies the validity of R — > T, hence the diagram must be continued 
transitively. The statement VIqo — > IV expresses the facts that there exists 
a sequence 



lim Cg> = C s . 



(2.12) 




(2.13) 



2.3 Continuous changes of the Bianchi type 
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and that in each such a sequence the parameter h must necessarily tend to 
infinity. 

Further VI^ — > II holds for every h. Both statements hold analogously 
for type VII^. In MacCallum [146] a similar diagram is shown, but there it 
is only investigated which changes appear if some of the canonical structure 
constants are vanishing. So e.g., the different transitions from type VI^ to 
types II and IV are not contained in it. 




I - V 

Fig. 1, taken from [200] 
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2 CHANGES OF THE BIANCHI TYPE 



To complete the answer to the question posed above it must be added: To 
each transition shown in the diagram, there one can indeed find a spacetime 
V A in which it is realized. Let f(t) be a C°°-function with f(t) — for t < 
and f(t) > else. Then, e.g., 



is a C°°-metric whose slices V 3 (t) belong to type V and IV for t < and 
t > respectively. Presumably it is typical that in a neighbourhood of t — 
the curvature becomes singular at spatial infinity; this is at least the case 
for metric (2.14). And, for the transition R — > S the limiting slice belongs 
necessarily to type S. 

2.4 Physical conditions 

To obtain physically reasonable spacetimes one has at least to secure the 
validity of an energy inequality like (2.1), T 00 > \T a/3 \ in each Lorentz frame. 
Without this requirement the transition II — > I is possible. We prove that 
the requirement Too > alone is sufficient to forbid this transition. To this 
end let V4 be a manifold which in a certain synchronized system of reference 
possesses a flat slice V 3 (0), and for all t > has Bianchi type II-slices V 3 (t). 
Using the notations of eq. (2.11) we have: C), 3 = — C32 = 1 are the only 
non-vanishing canonical structure constants of Bianchi type II. With the 
exception A\ = —x 3 we have Af = 5?. 

By the help of eq. (2.11) one can calculate the structure constants C l jk (t). 
The flat-slice condition is equivalent to 



ds 2 = -dt 2 + dx 2 + e 2x dy 2 + e 2x (dz + xf(t) dyf 



(2.14) 



lim CUt) = 0. 



(2.15) 



2.4 Physical conditions 
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First we consider the special case x\ = a l (t) ■ x % (no sum), i. e. extensions of 
the coordinate axes. It holds x\j = 5* a 1 (no sum). The only non- vanishing 
Cj k (t) are Cl 3 = —C\ 2 = a(t), where a = a 2 ■ a 3 /a 1 . Then eq. (2.15) reads 
lim^o a(t) = 0. The coefficients g a b(t) have to be chosen such that the g^, 
according to eq. (2.8), remain positive definite and twice continuously differ- 
entiable. Let be the scalar curvature within the slices. > appears 
only in Bianchi type IX and in the Kantowski-Sachs models. But changes 
from these types are just the cases already excluded by global considerations. 
Inserting into the Einstein equation by means of the Gauss-Codazzi the- 
orem one obtains 

nT m = R Q0 - ^oo {4) R = \ {3) R +^g~ 1 -H, g = det 9ij , (2.16) 
^R being the scalar curvature within the slices, hence 

{3) R < and H = 0n [022,0033,0 - (023,o) 2 ] 

- 2 • 012 [033,0012,0 - 013,0023,0] + cyclic perm . (2.17) 

In our case if is a quadratic polynomial in x 3 whose quadratic coefficient 
reads 

2 

(2.18) 

Hence, for sufficiently large values x 3 and values t with da/dt ^ we have 
H < 0, and therefore T 00 < 0. 

Secondly we hint at another special case, namely rotations of the coordi- 
nate axes against each other, e.g. 

x\ = x 1 cos uj + x 2 sin oj , x 2 = x 2 cos uj — x 1 sin uj , x 3 = x 3 , uj = u(t) . 



"011 [011033 - (913)' 



da{t) 
dt 
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2 CHANGES OF THE BIANCHI TYPE 



There one obtains in analogy to eq. (2.18) a negative Too- Concerning the 
general case we have: loosely speaking, each diffeomorphism is a composite 
of such extensions and rotations. Hence, for each II — > I-transition one would 
obtain points with negative T o- 

Hence, the energy condition is a genuine restriction to the possible tran- 
sitions of the Bianchi types. Concerning the other transitions we remark: 
equations (2.16), (2.17) keep valid, and the remaining work is to examine the 
signs of the corresponding expressions H. Presumably one always obtains 
points with negative H, i.e., also under this weakened homogeneity presump- 
tions the Bianchi types of the spacelike hypersurfaces at different distances 
from the singularity must coincide. 



19 



3 Inhomogeneous cosmological models with 
flat slices 

A family of cosmological models is considered which in a certain synchronized 
system of reference possess flat slices t = const. They are generated from the 
Einstein-De Sitter universe by a suitable transformation. Under physically 
reasonable presumptions these transformed models fulfil certain energy con- 
ditions. In Wainwright [265] a class of inhomogeneous cosmological models 
is considered which have the following property: there exists a synchronized 
system of reference of such a kind that the slices t = const, are homogeneous 
manifolds. 

Here we consider a special family of such models which possess flat slices. 
Following [201], we use the transformation formalism developed in chapter 
2. In Krasinski [137], e.g. in Fig. 2.1., it is shown how this approach fits 
into other classes of cosmological solutions of the Einstein field equation. 
Additionally, we require that these transformations leave two coordinates 
unchanged; this implies the existence of a 2-dimensional Abelian group of 
motions. A similar requirement is posed in [265], too. 

Starting from a Friedmann universe, we investigate whether the energy 
inequalities are fulfilled in the transformed model, too. In general this fails 
to be the case, but starting from the Einstein-De Sitter universe [68], cf. 
also Tolman [256] , and requiring perfect fluid for the transformed model, the 
energy inequalities in the initial model imply their validity in the transformed 
model. A new review on exact solutions can be found in [26] , and a new proof 
of Birkhoff's theorem is given in [226]. 
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3 INHOMOGENEOUS MODELS WITH FLAT SLICES 



3.1 Models with flat slices 

The transformation formalism of chapter 2 restricted to Bianchi type I reads 
as follows: using the same notations eq. (2.5), we have A? = 5?, u a = dx a 
and 

ds 2 = -dt 2 + g ah (t)dx a dx b (3.1) 

as the initial hypersurface-homogeneous model instead of eq. (2.4). Now let 
us consider the time-dependent transformation x1(x\t), where for each t it 
has to be a diffeomorphism of R 3 . Then one obtains in place of eq. (2.8) 

90a = Soa, gtj = gabtyxliXtj (3.2) 

as the transformed model. It is no restriction to insert g a b{t) = 5 a i> in eq. 
(3.1), i.e. to start from Minkowski spacetime. In the following we consider 
only transformations which leave two coordinates unchanged, i.e., now writ- 
ing t, x, y, z instead of x°, . . .x 3 resp., we restrict to transformations which 
read as follows 

x t (x,t), Dt = V, zt = z. (3.3) 
These transformations we shall call ^-transformations. 

Using ^-transformations, the Killing vectors d/dy and d/dz of the initial 
model remain Killing vectors. They form a 2-dimensional Abelian group of 
motions. All others, including the rotation z(d / dy) — y(d / dz) , may fail to re- 
main Killing vectors. On the contrary to the general case, the ^-transformed 
models depend genuinely on the initial ones. In the following, the spatially 
flat Friedmann universe with power-law behaviour of the cosmic scale factor 
shall be used as initial model, i.e. 

g ab (t)=5 ah K 2 (t) with K(t)=t T . (3.4) 



3.1 Models with Bat slices 
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Together with eq. (3.1) one obtains from the Einstein field equation 

2 

Kfx = kTqo = 3t 2 /t 2 , p = T 2 2 = a fx with a — 1 ; (3.5) 

3r 

where fx is the energy density and p is the pressure. 

Now, inserting eq. (3.4) in eq. (3.1) and transforming to eq. (3.2) with 
restriction eq. (3.3) one obtains for the metric of the x-transformed model 

gu = t 2r ■ (x t ,i) 2 = t 2r ■ h(x, t) , g 22 = 933 = t 2r , = ri^ else . (3.6) 

Here, r\ a p = diag(— 1, 1, 1, 1) is the metric of the flat Minkowski spacetime. 
This metric (3.6) belongs to the so-called Szekeres class, cf. [251]. Defining 
a(x, t) by 

5ii = e Mx ^ , (3.7) 
a coordinate transformation x(x) yields a(x, 1) = 0. If v — dot, then we have 

a(x,t) = J*v(x,i) -t^dt, (3.8) 

v(x,t) being an arbitrary twice continuously differentiable function which 
may be singular at t — and t — oo. The initial model is included by setting 
x t = x, hence a — r In t, v — r. 

For r / different functions v correspond to the same model only if 
they are connected by a translation into x-direction. Inserting eq. (3.6) 
with eq. (3.7) and eq. (3.8) into the Einstein field equation one obtains the 
energy-momentum tensor 

kT 00 = {r 2 + 2vr)/t 2 , 
kT\ = (2r-3r 2 )/t 2 , 
kT 2 = /tT 3 3 = -v fl r l + (r - t 2 - vt - v 2 + v) /t 2 , 
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3 INHOMOGENEOUS MODELS WITH FLAT SLICES 



KT a p = else , and 
kT = kT« = -2u r x - 2(3r 2 - 2r + 2wr + v 2 - v)/t 2 . (3.9) 

The question, in which cases equation (3.6) gives a usual hypersurface- 
homogeneous model, can be answered as follows: metric (3.6) is a Fried- 
mann universe, if and only if h = 0. For this case it is isometric to the 
initial model. Metric (3.6) is a hypersurface-homogeneous model, if and only 
if functions A and B exist for which holds h(x,t) — A(x) • B(t). Because of 
h > this is equivalent to h i01 ■ h = h j0 - h :1 . In this case it is a Bianchi type 
I model. 

3.2 Energy inequalities and perfect fluid models 

In this section it shall be discussed, in which manner the geometrically de- 
fined models described by eqs. (3.6), (3.7), (3.8) and (3.9) fulfil some energy 
conditions. Here we impose the following conditions: each observer measures 
non-negative energy density, time- or lightlike energy flow and spacelike ten- 
sions which are not greater than the energy density. In our coordinate system 
these conditions are expressed by the following inequalities 

Too > l^i, (3.io) 

Too > |T|| , (3.11) 

and T<0. (3.12) 

For the initial model this means r = or r > 1/2, i.e. Minkowski spacetime 
or —1 < a < 1/3. For the limiting case a — — 1 one gets the de Sitter 
spacetime in the original model (3.5) which does not have the assumed power- 
law behaviour for the scale factor. Using the energy-momentum tensor eq. 
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(3.9), eqs. (3.11) and (3.12) read 

v,ot > -v 2 - (3t- 1> + t-2t 2 , (3.13) 

v,o t < -v 2 + (r + l)v + r and (3.14) 
v,o t>-v 2 - (2r - l)u + 2r - 3r 2 . (3.15) 

Now, if r < 0, i.e. a < —1, then eq. (3.10) reads v < r — 1; together with 
eq. (3.14) one obtains i> j0 £ < r — 2. This implies the existence of a t > with 
v(i) > —1 in contradiction to condition (3.10). Therefore, an initial model 
with r < 0, which itself contradicts the energy inequalities, cannot produce 
transformed models which always fulfil them. 

If r = 0, then eqs. (3.13) and (3.14) imply v j0 t — v — v 2 . This equation 
has the solutions v — and v — t(t + C) _1 with arbitrary C(x) > 0. This 
yields a = and a = ln(t + C) — ln(l + C) with eq. (3.7), and gu = 1 and 
gu = (t + C) 2 ■ (1 + Cy 2 resp. with eq. (3.8). Then eq. (3.6) shows that 
this is the Minkowski spacetime itself. Therefore, the Minkowski spacetime 
does not produce any new models. 

Finally, if r > 0, hence a > — 1, eq. (3.10) then reads v > max(r —1,1 — 
2t). A lengthy calculation shows which transformed models fulfil the energy 
inequalities. For each r > models exist which do and models which do not 
fulfil them. 

The situation described above changes if one requires that the trans- 
formed model consists of perfect fluid with an equation of state. The velocity 
vector must be (1, 0, 0, 0) and 

n(Tl — T 2 ) = . \ + * (V) =0 (3.16) 
V 1 2! (l + a)t-h 2Vh \VhJ t0 
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must be fulfilled. Eq. (3.16) expresses the condition that the pressure has to 
be isotropic. T| = T 3 3 is fulfilled anyhow due to our assumption eq. (3.3). 
If / = h fi h~ 1/2 , then eq. (3.16) reads 



(l + a)t 2 

hence f$f~ l does not depend on x, therefore /^i • f — f,o • f,i, and we can 
use the ansatz / = a(t) ■ b(x). Inserting this into eq. (3.16), one obtains 

h = [b(x) ■ t {a ' l){a+1) + c(x)f , where c(x) = yjh(x,0) (3.17) 

with arbitrary non-negative functions b and c fulfilling b(x) + c(x) > for all 
x. For energy density and pressure we then obtain 

4 4(1 -a) -b 



3(1 + a)H 2 3(l + «) 2 t(6t + ct 2 /(i+«)) 



4a 

Kp= 3(l + a)^- (3 - 18) 



An equation of state means that p uniquely depends on /i. This takes 
place if and only if a = or b/c = const. The latter is equivalent to the 
hypersurface-homogeneity of the model and is of lower interest here. For 
a = the initial model is the dust-filled Einstein-De Sitter model [68]. 

With eq. (3.18) we obtain p = and fx > and may formulate: If the 
Einstein-De Sitter universe is ^-transformed into a perfect fluid model, then 
this model also contains dust and fulfils the energy conditions. 

These models have the following form: inserting eq. (3.17) with a = 
into eq. (3.6) we get a dust-filled model 

ds 2 = -dt 2 + t 4 / 3 { [b(x)/t + c(x) f dx 2 + dy 2 + dz 2 } (3.19) 
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with arbitrary 6, c as before. A subcase of eq. (3.19) is contained in Szekeres 
[251] as case (iii), but the parameter e used there may now be non-constant. 
For 6 = 1, eq. (3.19) reduces to eq. (2.9). Of course, locally one can achieve 
a constant 6 by a transformation x = x(x), but, e.g. changes in the sign of 6 
then are not covered. As an illustration we give two examples of this model 
eq. (3.19): 



1. If 6 = 1, c> 0, then 



hence the density contrast at two different values x±, x 2 reads 

fii c 2 1 + Cit 



(3.21) 



fl 2 Ci 1 + c 2 t 

and tends to 1 as t — > oo. This shows that one needs additional presumptions 
if one wants to prove an amplification of initial density fluctuations. 

2. If 6 = 1, c = 0, then eq. (3.19) is the axially symmetric Kasner vacuum 
solution. If now c differs from zero in the neighbourhoods of two values x±, 
x 2 , then the model is built up from two thin dust slices and Kasner-like 
vacuum outside them. The invariant distance of the slices is 

At~ 1/3 + Bt 2/3 (3.22) 

with certain positive constants A and B. In [225] it is argued that the in- 
variant distance along spacelike geodesies need not coincide with what we 
measure as spatial distance. The distinction comes from the fact, that a 
geodesic within a prescribed spacelike hypersurface need not coincide with 
a spacelike geodesic in spacetime. This eq. (3.22) looks like gravitational 
repulsion, because the distance has a minimum at a positive value. But the 
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t 1 / 3 -term is due to the participation of the slices in cosmological expansion 
and the remaining; -term is due to an attractive gravitational force in 
parabolic motion. A more detailed discussion of this behaviour eq. (3.22) is 
given in Griffiths [99]: he argues, "however inhomogeneous the mass distri- 
bution, matter on each plane of symmetry has no net attraction to matter 
on other planes." 

The transformation of a hypersurface-homogeneous cosmological model 
considered here preserves all inner properties, expressed by the first funda- 
mental form, of the slices t = const., and also preserves the property that t 
is a synchronized time, but may change all other outer properties, which are 
expressed by the second fundamental form. The investigations of section 3.2 
show that energy conditions are preserved under very special presumptions 
only. 

One may consider these transformations as a guide in the search for new 
exact solutions of Einstein's field equation. The new models are close to the 
initial hypersurface-homogeneous ones, if the transformation is close to the 
identical one. Thus, one can perturb a Bianchi model with exact solutions 
without use of any approximations. 

3.3 A simple singularity theorem 

What do we exactly know about solutions when no exact solution, in the sense 
of "solution in closed form" , is available? In which sense these solutions have 
a singularity? Here we make some remarks to these questions, see [228]. 

In [221], the following simple type of singularity theorems was discussed: 
The coordinates t, x, y, z shall cover all the reals, and a(t) shall be an 
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arbitrary strictly positive monotonously increasing smooth function defined 
for all real values t, where "smooth" denotes "C°°-differentiable" . Then it 
holds: 

The Riemannian space defined by 

ds 2 = dt 2 + a 2 {t){dx 2 + dy 2 + dz 2 ) (3.23) 

is geodesically complete. This fact is well-known and easy to prove; however, 
on the other hand, for the same class of functions a(t) it holds: 

The spacetime defined by 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ) (3.24) 

is lightlike geodesically complete if and only if 

,o 

/ a(t) dt = oo (3.25) 

J —oo 

is fulfilled. The proof is straightforwardly done by considering lightlike 
geodesies in the x — t— plane. So, one directly concludes that the state- 
ments do not depend on the number n of spatial dimensions, only n > is 
used, and the formulation for n = 3 was chosen as the most interesting case. 
Moreover, this statement remains valid if we replace "lightlike geodesically 
complete" by "lightlike and timelike geodesically complete" . 

Further recent results on singularity theorems are reviewed in Senovilla 
[233]. This review and most of the research concentrated mainly on the 4- 
dimensional Einstein theory. Probably, the majority of arguments can be 
taken over to the higher-dimensional Einstein equation without change, but 
this has not yet worked out up to now. And singularity theorems for F(R)- 
gravity are known up to now only for very special cases only. 
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The most frequently used F(i?)-Lagrangian is 




where / > , 



(3.26) 



here we discuss the non-tachyonic case only. From this Lagrangian one gets 
a fourth-order field equation; only very few closed-form solutions, "exact 
solutions", are known. However, for the class of spatially flat Friedmann 
models, the set of solutions is qualitatively completely described, but not in 
closed form, in [163]. We call them "rigorous solutions". 

One of them, often called "Starobinsky inflation", can be approximated 
with 



This approximation is valid in the region t <C —I. However, this solution 
does not fulfil the above integral condition. Therefore, Starobinsky inflation 
does not represent a lightlike geodesically complete cosmological model as 
has been frequently stated in the literature. 

To prevent a further misinterpretation let me reformulate this result as 
follows: Inspite of the fact that the Starobinsky model is regular, in the 
sense that a(t) > for arbitrary values of synchronized time t, every past- 
directed lightlike geodesic terminates in a curvature singularity, i.e., \R\ — > 
oo, at a finite value of its afline parameter. Therefore, the model is not only 
geodesically incomplete in the coordinates chosen, moreover, it also fails to 
be a subspace of a complete one. 




(3.27) 



In contrast to this one can say: the spatially flat Friedmann model with 
a(t) = exp(Ht), H being a positive constant, is the inflationary de Sitter 
spacetime. According to the above integral condition, it is also incomplete. 
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However, contrary to the Starobinsky model, it is a subspace of a complete 
spacetime, namely of the the de Sitter spacetime represented as a closed 
Friedmann model. 

The application of conformal transformations represents one of the most 
powerful methods to transform different theories into each other. In many 
cases this related theories to each other which have been originally considered 
to be independent ones. The most often discussed question in this context is 
"which of these metrics is the physical one" , but one can, of course, use such 
conformal relations also as a simple mathematical tool to find exact solutions 
without the necessity of answering this question. However, concerning a con- 
formal transformation of singularity theorems one must be cautious, because 
typically, near a singularity in one of the theories, the conformal factor di- 
verges, and then the conformally transformed metric need not be singular 
there. 

Finally, let me again mention the distinct properties of the Euclidean and 
the Lorentzian signature: Smooth connected complete Riemannian spaces are 
geodesically connected. However, smooth connected complete spacetimes, 
e.g. the de Sitter space-time, need not be geodesically connected, see for 
instance Sanchez [196]. The topological origin of this distinction is the same 
as the distinction between the two above statements: for the positive definite 
signature case one uses the compactness of the rotation group, whereas in the 
indefinite case the noncompactness of the Lorentz group has to be observed 
[227]. 
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4 Properties of curvature invariants 

In section 4.1, which is an extended version of [223] we answer the follow- 
ing question: Let A, //, v be arbitrary real numbers. Does there exist a 
3-dimensional homogeneous Riemannian manifold whose eigenvalues of the 
Ricci tensor are just A, fi and ul 

In section 4.2 we discuss, following [224], why all the curvature invariants 
of a gravitational wave vanish. 

In section 4.3 we present the example [222] that it may be possible that 
non-isometric spacetimes with non-vanishing curvature scalar cannot be dis- 
tinguished by curvature invariants. 

4.1 The space of 3— dimensional homogeneous Rieman- 
nian manifolds 

The curvature of a 3-manifold is completely determined by the Ricci tensor, 
and the Ricci tensor of a homogeneous manifold has constant eigenvalues. 
However, it is essential to observe, that, nevertheless, the constancy of all 
three eigenvalues of the Ricci tensor of a 3-manifold does not imply the 
manifold to be locally homogeneous, see [132] for examples. 

Assume A, fi, v to be arbitrary real numbers. First question: Does there 
exist a 3-dimensional homogeneous Riemannian manifold whose eigenvalues 
of the Ricci tensor are just A, \i and v? Second question: If it exists, is 
it uniquely determined, at least locally up to isometries? The answer to 
these questions does not alter if we replace "homogeneous" by "locally ho- 
mogeneous". In principle, they could have been answered by the year 1905 
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already. But in fact, it was not given until 1995: the first of these problems 
have been solved independently by Kowalski and Nikcevic [133] on the one 
hand, and by Rainer and Schmidt [189], on the other hand; the second one 
is solved several times, e.g. in [189]. Curiously enough, the expected answer 
"Yes" turns out to be wrong. 

Of course, both answers [133] and [189] are equivalent. But the formu- 
lations are so different, that it is useful to compare them here. In [133] one 
defines 

pi = max{A,/i, v) , p 3 = min{A, p, u} (4.1) 

and 

p 2 = A + p + v - pi - p 3 (4.2) 
and formulates the conditions as inequalities between the pi. 

In [189] the formulation uses curvature invariants. First step: one ob- 
serves that the answer does not alter, if one multiplies A, p, v with the same 
positive constant. Note: this fails to be the case if we take a negative con- 
stant. Second step: we take curvature invariants which do not alter by this 
multiplication, i.e. we look for suitable homothetic invariants. Third step: 
we formulate the answer by use of the homothetic invariants. 

First step: this can be achieved by a homothetic transformation of the 
metric, i.e., by a conformal transformation with constant conformal factor. 
Then one excludes the trivial case A = p — v — 0. Note: It is due to the 
positive definiteness of the metric that A = p = v = implies the flatness of 
the space. 

Second step: 

R = \ + p + v, N = X 2 + p 2 + u 2 (4.3) 
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are two curvature invariants with N > 0, and 

R = R/y/N (4.4) 

is a homothetic invariant. For denning a second nomothetic invariant we 
introduce the trace-free part of the Ricci tensor 

R 

Sij = Rij - —g-ij (4.5) 

and its invariant 

S = SiS*S* k . (4.6) 
The desired second homothetic invariant reads 

S = S/y/N*. (4.7) 

Of course, one can also directly define R and S as functions of A, p and v, 
then S is the cubic polynomial, 

S = ((2A - p - vf + (2fi-u- A) 3 + {2u - A - uf) /27 , (4.8) 

see [189] for details. Remark: the inequalities S > and p2 < (pi + P3)/2 
are equivalent. 

Third step: The answer is "yes" if and only if one of the following condi- 
tions (4.9) till (4.14) are fulfilled: 

N = 0, flat space (4.9) 

R 2 = 3 and 5 = 0, non-flat spaces of constant curvature (4.10) 

R 2 = 2 and 9^5" = —1 , real line times constant curvature space (4.11) 

R 2 = 1 and 9RS = 2 (4.12) 

18S > R{9 - 5R 2 ) , equivalent to Apz/ > (4.13) 

2 g 5J? 2 

-V3 <R< -V2 and < S < (3 - R 2 ) — . (4.14) 

~ V > 18R 3 
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Remarks: 1. In [133] the last of the four inequalities (4.14) is more elegantly 
written as 

ft >^. (415) 

Pi + P3 

2. There exists a one-parameter set of Bianchi-type VI manifolds, with 
parameter a where < a < 1, having p\ = —a— a 2 , p 2 = —I— a 2 , p 3 = —1—a; 
they represent the case where in (4.14) the fourth inequality is fulfilled with 



Now we turn to the second question: Let the 3 numbers A, p, v be given 
such that a homogeneous 3-space exists having these 3 numbers as its eigen- 
values of the Ricci tensor; is this 3-space locally uniquely determined? The 
first part of the answer is trivial: 

If A = p — v, then "yes" , it is the space of constant curvature. The second 
part of the answer is only slightly more involved: If A, p and v represent 3 
different numbers, then the 3 eigendirections of the Ricci tensor are uniquely 
determined, and then we can also prove: "yes". The third part is really 
surprising: If 

R=l and S = ^ (4.16) 
then the answer is "no" , otherwise "yes" . 

Clearly, eq. (4.16) represents a subcase of (4.12). What is the peculiarity 
with the values (4.16) of the homothetic invariants? They belong to Bianchi 
type IX, have one single and one double eigenvalue, and with these values, 
a one-parameter set of examples with 3-dimensional isometry group, and 
another example with a 4-dimensional isometry group exist. 
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4.2 Why do all the curvature invariants of a gravita- 
tional wave vanish? 

We prove the theorem valid for Pseudo-Riemannian manifolds V n : "Let x G 
V n be a fixed point of a homothetic motion which is not an isometry then 
all curvature invariants vanish at x." and get the Corollary: "All curvature 
invariants of the plane wave metric 

ds 2 = 2dudv + a 2 (u)dw 2 + b 2 (u)dz 2 (4.17) 

identically vanish." 

Analysing the proof we see: The fact that for definite signature flatness 
can be characterized by the vanishing of a curvature invariant, essentially 
rests on the compactness of the rotation group SO(n). For Lorentz signature, 
however, one has the non-compact Lorentz group 5*0(3, 1) instead of it. 

A further and independent proof of the corollary uses the fact, that the 
Geroch limit does not lead to a Hausdorff topology, so a sequence of gravita- 
tional waves can converge to the flat spacetime, even if each element of the 
sequence is the same pp-wave. 

The energy of the gravitational field, especially of gravitational waves, 
within General Relativity was subject of controversies from the very be- 
ginning, see Einstein [66]. Global considerations - e.g. by considering the 
far-field of asymptotically flat spacetimes - soon led to satisfactory answers. 
Local considerations became fruitful if a system of reference is prescribed e.g. 
by choosing a timelike vector field. If, however, no system of reference is pre- 
ferred then it is not a priori clear whether one can constructively distinguish 
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flat spacetime from a gravitational wave. This is connected with the gener- 
ally known fact, that for a pp-wave, see e.g. Stephani [250] especially section 
15.3. and [65] all curvature invariants vanish, cf. Hawking and Ellis [107] 
and Jordan et al. [123], but on the other hand: in the absence of matter or 
reference systems - only curvature invariants are locally constructively mea- 
surable. See also the sentence "R. Penrose has pointed out, in plane- wave 
solutions the scalar polynomials are all zero but the Riemann tensor does not 
vanish." taken from page 260 of [107]. At page 97 of [123] it is mentioned 
that for a pp-wave all curvature invariants constructed from 

Rijki;h...i r (4-18) 
by products and traces do vanish. 

It is the aim of this section to explain the topological origin of this strange 
property. 

4.2.1 Preliminaries 

Let V n be a C°°-Pseudo-Riemannian manifold of arbitrary signature with 
dimension n > 1. The metric and the Riemann tensor have components 
and Rijim resp. The covariant derivative with respect to the coordinate x m 
is denoted by ";m" and is performed with the Christoffel affinity Y\ m . We 
define 

/ is called a generalized curvature invariant of order k if it is a scalar with 
dependence 

I I(.9iji Rijlmi ■ ■ ■ i Rijlm;i\...ik) ■ (4.19) 
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By specialization we get the usual Definition: / is called a curvature invariant 
of order k if it is a generalized curvature invariant of order k which depends 
continuously on all its arguments. The domain of dependence is requested 
to contain the flat space, and 

I( 9ij , 0,...0) = 0. (4.20) 

Examples: Let 

n 

J = sign( Yl \Rijim\), (4.21) 

i,j,l,m=l 

where the sign function is defined by sign(0) = and sign(x) = 1 for x > 0. 
Then I is a generalized curvature invariant of order 0, but it fails to be a 
curvature invariant. It holds: V n is flat if and only if I = 0. Let further 

h = R ij imR ijlm (4-22) 

which is a curvature invariant of order 0. If the metric has definite signature 
or if n = 2 then it holds: V n is flat if and only if I\ = 0. Proof: For definite 
signature Iq = sign^); for n — 2, Ii = implies R = 0, hence flatness, 
q.e.d. 

For all other cases, however, the vanishing of I\ does not imply flatness. 
Moreover, there does not exist another curvature invariant serving for this 
purpose, it holds: 

For dimension n > 3, arbitrary order k and indefinite metric it holds: To 
each curvature invariant / of order k there exists a non-flat V n with 7 = 0. 

Proof: Let n — 3. We use 



ds 2 = Idudv ± a 2 (u)dw 2 



(4.23) 
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with a positive non-linear function a(u). The "±" covers the two possible 
indefinite signatures for n — 3. The Ricci tensor is Rij = R m , ■ and has 
with u = x 1 

R - = -;•*? (424) 

and therefore, eq. (4.23) represents a non-flat metric. Now let n > 3. We use 
the Cartesian product of eq. (4.23) with a flat space of dimension n — 3 and 
arbitrary signature. So we have for each n > 3 and each indefinite signature 
an example of a non-flat V^. It remains to show that for all these examples, 
all curvature invariants of order k vanish. It suffices to prove that at the 
origin of the coordinate system, because at all other points it can be shown 
by translations of all coordinates accompanied by a redefinition of a(u) to 
a(u — u ). Let I be a curvature invariant of order k. Independent of the 
dimension, i.e., how many flat spaces are multiplied to metric eq. (4.23), one 
gets for the case considered here that 

/ = I(a {0 \u),a^(u),...,a {k+2 \u)) (4.25) 

where a^(u) = a(u), a (m+1 \u) = £a^(u), and 

J(a (0) (M),0,..., 0) = 0. (4.26) 

This is valid because each 

Rijlm;ii...i p (4.27) 

continuously depends on a^(u), a^(u), . . . , a^+^-u) and on nothing else; 
and for a = const., metric (4.23) represents a flat space. 



Now we apply a coordinate transformation: Let e > be fixed, we replace 
u by u ■ e and v by v/e. This represents a Lorentz boost in the u — v— plane. 
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Metric eq. (4.23) remains form-invariant by this rotation, only a(u) has to 
be replaced by a(u ■ e). At u = we have 

/ = J(a (0) (0),a (1) (0),..., a (fc+2) (0)) (4.28) 

which must be equal to 

I £ = /( a (°)(0),e.aW(0),..., e k+2 ■ a^ k+2 \0)) (4.29) 

because I is a scalar. By continuity and by the fact that flat space belongs 
to the domain of dependence of /, we have 

lim I e = . (4.30) 

All values I £ with e > coincide, and so I — 0. q.e.d. 

4.2.2 Gravitational waves 

New results on gravitational waves can be found in [72], and [84]. A pp- 
wave is a plane- fronted gravitational wave with parallel rays, see [77]. It is 
a non-flat solution of Einstein's vacuum equation = possessing a non- 
vanishing covariantly constant vector; this vector is then automatically a null 
vector. The simplest type of pp-waves can be represented similar as metric 
eq. (4.23) 

ds 2 = 2dudv + a 2 (u)dw 2 + b 2 (u)dz 2 , (4.31) 

where 

d 2 a d 2 b 

This metric represents flat spacetime if and only if both a and b are linear 
functions. Using the arguments of subsection 4.2.1 one sees that all curvature 
invariants this metric identically vanish. Here we present a second proof of 
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that statement which has the advantage to put the problem into a more 
general framework and to increase the class of spacetimes covered, e.g. to 
the waves 

ds 2 = dx 2 + dy 2 + 2du dv + H(x, y, u)du 2 (4.33) 

Hall [102] considers fixed points of homothetic motions, which cannot exist 
in compact spacetimes, and shows that any plane wave, not only vacuum 
plane waves, admits, for each x, a homothetic vector field which vanishes at 
x. It holds 

Let x G V n be a fixed point of a homothetic motion which is not an 
isometry then all curvature invariants vanish at x. 

Proof: The existence of a homethetic motion which is not an isometry 
means that V n is selfsimilar. Let the underlying differentiable manifold be 
equipped with two metrics gij and 

9a = e 2C g tJ (4.34) 

where C is a no n- vanishing constant. The corresponding Riemannian mani- 
folds are denoted by V n and V n resp. By assumption, there exists an isometry 
from V n to V n leaving x fixed. Let J be a curvature invariant. / can be rep- 
resented as continuous function which vanishes if all the arguments do of 
finitely many of the elementary invariants. The elementary invariants are 
such products of factors g 1 ^ with factors of type 

Rijlm;i 1 ...i p (4.35) 

which lead to a scalar, i.e., all indices are traced out. Let J be such an 
elementary invariant. By construction we have J(x) = e qC J(x) with a non- 
vanishing natural q, which depends on the type of J. Therefore, J(x) = 0. 
q.e.d. 
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From this it follows: All curvature invariants of metric (4.17) identically 
vanish. This statement refers not only to the 14 independent elementary 
invariants of order 0, see Harvey [104] and Lake 1993 [140] for a list of them, 
but for arbitrary order. 

Proof: We have to show that for each point x, there exists a homothetic 
motion with fixed point x which is not an isometry. But this is trivially done 
by suitable linear coordinate transformations of v, w, and z. q.e.d. 

4.2.3 Topological properties 

Sometimes it is discussed that the properties of spacetime which can be lo- 
cally and constructively, i.e., by rods and clocks, measured are not only the 
curvature invariants but primarily the projections of the curvature tensor 
and its covariant derivatives to an orthonormal tetrad, called 4-bein. 6 The 
continuity presumption expresses the fact that a small deformation of space- 
time should also lead to a correspondingly small change of the result of the 
measurement. To prevent a preferred system of reference one can construct 
curvature invariants like 

I 2 = inf ]T \R l3lm \ (4.36) 

i,j,l,m 

where the infimum, here the same as the minimum, is taken over all orthonor- 
mal tetrads. From the first glance one could believe that J 2 = if and only 
if the space is flat. But for indefinite signature this would contradict the 
result of subsection 4.2.1. What is the reason? For definite signature the 
infimum is to be taken about the rotation group 5*0(4), or 0(4) if one allows 

6 The German word "das Bern" denotes both "bone" and "leg", here it is used in the 
second sense. 
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orientation- reversing systems; this group is compact. One knows: A positive 
continuous function over a compactum possesses a positive infimum. So, if 
one of the Riji m differs from zero, then J 2 > at that point. For Lorentz sig- 
nature, however, the infimum is to be taken about the non-compact Lorentz 
group SO(3, 1) and so Riji m ^ does not imply I 2 ^ 0. 

Another topological argument, which underlies our subsection 4.2.1, is 
connected with the Geroch limit of spacetimes [86], we use the version of 
[210]. Theorem 3.1 of reference [210] reads: 

1: For local Riemannian manifolds with definite signature, Geroch's limit 
defines a Hausdorff topology. 

2: For indefinite signature this topology is not even T\. 

Explanation: A topology is Hausdorff if each generalized Moore - Smith 
sequence possesses at most one limit, and it is Ti if each constant sequence 
possesses at most one limit. The main example is a sequence, where each 
element of the sequence is the same pp-wave, and the sequence possesses two 
limits: flat spacetime and that pp-wave. Here the reason is: Only for definite 
signature, geodetic e-balls form a neighbourhood basis for the topology. 

Final remarks: The change from Euclidean to Lorentzian signature of a 
Pseudo- Riemannian space is much more than a purely algebraic duality - an 
impression which is sometimes given by writing an imaginary time coordi- 
nate: One looses all the nice properties which follow from the compactness of 
the rotation group. Lake 1993 [140] pointed out that also for the Robinson- 
Trautman vacuum solutions of Petrov type III all 14 curvature invariants 
vanish. 
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4.3 On spacetimes which cannot be distinguished by 
curvature invariants 

For a positive C°°-function a(u) let 

1 



ds 2 = - 2 



2dudv - a 2 {u)dy 2 - dz 2 . (4.37) 



In the region z > 0, ds 2 represents the anti-de Sitter spacetime if and only if 
a(u) is linear in u. Now, let d 2 a/du 2 < and 



Then 



and 



□0 = <f) yi r = (4.39) 



R 

R H ~ 7f9ij = ^9ij + KTij (4.40) 

with A = — 3andT^ = 4>,i4>,j- So (ds 2 , 0) represents a solution of Einstein's 
equation with negative cosmological constant A and a minimally coupled 
massless scalar field <fi. 

Let / be a curvature invariant of order k, i.e., / is a scalar 

I I (,9iji Rijlmi ■ ■ ■ i ■ R ijlm;ii...ik^) (^-41) 

depending continuously on all its arguments. Then for the metric ds 2 , I does 
not depend on the function a(u). 



This seems to be the first example that non-isometric spacetimes with 
non-vanishing curvature scalar cannot be distinguished by curvature invari- 
ants. The proof essentially uses the non-compactness of the Lorentz group 
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50(3, 1), here of the boosts m->mA, ?;->w/A. One can see this also in the 
representation theory in comparison with the representations of the compact 
rotation group SO (A). 

Let v G R A be a vector and g G 5*0(4) such that g(v) || v, then g(v) = v. 
In Minkowski spacetime M 4 , however, there exist vectors v G M 4 and a /i G 
50(3, 1) with ft,(i>) || f and h(v) ^ v. This is the reason why the null frame 
components of the curvature tensor called Cartan "scalars" are not always 
curvature scalars, see [224]. Further new results on curvature invariants are 
given in [28], [184] and [188]. Limits of spacetimes are considered in [176] 
and [230]. 
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5 Surface layers and relativist ic surface ten- 
sions 

For a thin shell, the intrinsic 3-pressure will be shown to be analogous to 
—A, where A is the classical surface tension: First, interior and exterior 
Schwarzschild solutions will be matched together such that the surface layer 
generated at the common boundary has no gravitational mass; then its intrin- 
sic 3-pressure represents a surface tension fulfilling Kelvin's relation between 
mean curvature and pressure difference in the Newtonian limit. Second, after 
a suitable definition of mean curvature, the general relativistic analogue to 
Kelvin's relation will be proven to be contained in the equation of motion of 
the surface layer. 

5.1 Introduction to surface layers 

In general relativity, an energy-momentum tensor concentrated on a timelike 
hypersurface is called a surface layer. Via Einstein's equations it is related 
to non-spurious jumps of the Christoffel affinities or equivalently to jumps 
of the second fundamental tensor. In [141], [178], [60], [117], [177] and [138] 
there have been given algorithms for their calculating, and in [141], [117], 
and [38], [63], [202], [260], [79] and [140] the spherically symmetric case was 
of a special interest. The discussion of disklike layers as models for accretion 
disks was initiated in [160] and [264], and [143] contains surface layers as 
sources of the Kerr geometry. 

A surface layer has no component in the normal direction, otherwise the 
delta-like character would be destroyed. Hence, an ideal fluid with non- 
vanishing pressure cannot be concentrated on an arbitrarily thin region and 
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therefore in most cases, the layer is considered to be composed of dust. Ad- 
ditionally, an intrinsic 3-pressure is taken into account in [67], where for a 
spherically symmetric configuration each shell r = const, is thought to be 
composed of identical particles moving on circular orbits without a preferred 
direction, and the tangential pressure, which is analogous to the intrinsic 
3-pressure, is due to particle collisions, whereas a radial pressure does not 
appear. See also [202], [140], and [151], where the equation of motion for a 
spherically symmetric layer has been discussed. 

Such an intrinsic 3-pressure as well as surface tensions are both of the 
physical dimension "force per unit length." There the question arises whether 
an intrinsic 3-pressure of a surface layer may be related to a surface tension, 
and this chapter will deal with just this question. Then a general relativistic 
formulation of thermodynamics can be completed by equations for surface 
tensions to answer, e.g., the question how long a drop, say, a liquid comet, 
remains connected while falling towards a compact object. The only paper 
concerned with such questions seems to be [135]. There the influence of 
surface tensions on the propagation of gravitational waves has been calculated 
by perturbation methods, yielding a possibly measurable effect. For the 
background knowledge cf. [118] or [170] for relativistic thermodynamics and 
[174] for non- relativistic surface tensions. 

This chapter, which is based on [203], proceeds as follows: Section 5.2 
contains Kelvin's relation for non-relativistic surface tensions and discusses 
the matching of the interior to the exterior Schwarzschild solution such that 
a surface tension appears, and compares with Kelvin's relation in the New- 
tonian limit, calculations are found in the appendices, sections 5.5, 5.5 and 
5.5. Sections 5.3 and 5.4 are devoted to the non-spherically symmetric case. 
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Section 5.3 contains a suitable definition of mean curvature and section 5.4 
deduces Kelvin's relation from the equation of motion of the surface layer 
without any weak field assumptions. 

5.2 Non-relativistic surface tensions 

Imagine a drop of some liquid moving in vacuo. Its equilibrium configuration 
is a spherical one, and Kelvin's relation [255] between surface tension A, 
pressure difference AP, i.e. outer minus inner pressure, and mean curvature 
H, where H = l/R for a sphere of radius R, reads 



A is a material-dependent constant. Equation (5.1) means, an energy A-AF 
is needed to increase the surface area by AF. This supports our description 
of —A as a kind of intrinsic pressure. But of course, it is a quite different 
physical process: Pressure, say, of an ideal gas, can be explained by collisions 
of freely moving particles, whereas the microphysical explanation of surface 
tensions requires the determination of the intermolecular potential, which 
looks like 



with certain constants fi, N, and p; see Ono [174]. In this context the surface 
has a thickness of about 1CT 7 cm but in most cases this thickness may be 
neglected. 

In addition, for a non-spherically symmetric surface, the mean curvature 
H in equation eq. (5.1) may be obtained from the principal curvature radii 
Ri, R 2 by means of the relation 



AP = -2HA . 





(5.2) 



1 



1 



(5.3) 



H = 



+ 



2P X 



2R 2 
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Now the drop shall be composed of an incompressible liquid. For a general 
relativistic description we have to take the interior Schwarzschild solution 



ds 2 = - 

dr 2 



\{l-rjr,fl 2 -\{l-r g r 2 /rlfl 22 



dt 



2 



H Q . „ + r 2 dtt 2 where dtt 2 = dip 2 + sin 2 tpdcp 2 , (5.4) 

1 — r g r 2 /Tq 

whose energy-momentum tensor represents ideal fluid with energy density 
\i = 3r g /nrQ and pressure 

(l- rgr y r 3)l/2_ (1,^)1/2 

Plrj ^ ' 3(1 - r 9 /r )V2 _ (i _ rgr 2/ r 3)i/2 ■ ^ 
One has p(r ) = 0, and therefore usually < r < r is considered. But now 
we require only a non-vanishing pressure at the inner surface and take eq. 
(5.4) for values r with < r < R and a fixed R < r only. 

The gravitational mass of the inner region equals 

M = n ■ 4?r J R 3 /3 = r g (R/r f/2 (5.6) 

where we assumed unit such that G = c = 1. The outer region shall be 
empty and therefore we have to insert the Schwarzschild solution for r > R. 
Neglecting the gravitational mass of the boundary S C V 4 which is defined 
by r — R = 0, just M of equation (5.6) has to be used as the mass parameter 
of the exterior Schwarzschild metric. Then only delta-like tensions appear 
at S, and AP = —p(R) < 0. The vanishing of the gravitational mass is 
required to single out the properties of an intrinsic 3-pressure. In general, a 
surface layer is composed of both parts, not at least to ensure the validity of 
the energy condition T 00 > jT^ which holds for all known types of matter. 

On S, the energy- momentum tensor is 

T? = t* ■ 5 E (5.7) 
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the non-vanishing components of which are 

r 2 2 = r | = - p (R)R/2(l - 2M/R) 1/2 (5.8) 

cf. section 5.5. 

The 5s distribution is defined such that for all smooth scalar functions / 
the invariant integrals are related by 

/ fd^x = f f.SxdWx. (5.9) 

For a more detailed discussion of distribution- valued tensors in curved space- 
time, cf. [252]. The fact that all components Tq in eq. (5.7) vanish reflects 
the non-existence of a delta-like gravitational mass on E. 

Now consider the Newtonian limit M/R <C 1; then the mean curvature 
becomes again H = 1/R, and together with eq. (5.8) Kelvin's relation eq. 
(5.1) is just equivalent to 

r| = r| = -A[\ + 0(M/R)} . (5.10) 

Therefore: At least for static spherically symmetric configurations and weak 
fields a delta-like negative tangential pressure coincides with the classical 
surface tension. 

In the next two sections we investigate to what extent these presumptions 
are necessary. 

5.3 Mean curvature in curved spacetime 

To obtain a general relativistic analogue to Kelvin's relation eq. (5.1) we have 
to define the mean curvature H of the timelike hypersurface E contained in 
a spacetime V4 such that for weak fields just the usual mean curvature arises. 
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To get the configuration we have in mind, we make the ansatz 

r afS = -A(g af3 + u a up) u a u a = -1 (5.11) 

with A > 0. Thereby again the gravitational mass of £ will be neglected. 
Now the mean curvature shall be defined. But there is a problem: In general, 
there does not exist a surface ScS which can serve as "boundary at a fixed 
moment" for which we are to determine the mean curvature. To circumvent 
this problem we start considering the special case 

U a \\i3 = Up\\ a . (5-12) 

Then there exists a scalar t on £ such that u a = t\\ a , and the surface ScE 
defined by t = may be called "boundary at a fixed moment." Now S has to 
be embedded into a "space at a fixed moment" : We take intervals of geodesies 
starting from points of S in the normal direction rii and the opposite one. 

The union V 3 of these geodetic segments will be called "space at a fixed 
moment," and S C V3 is simply a two-surface in a three-dimensional positive 
definite Riemannian manifold, for which mean curvature has a definite sense: 
Let v a , w a be the principal curvature directions inside S and R±, R2 the 
corresponding principal curvature radii, then equation (5.3) applies to obtain 
H. 

Of course, v a w a = holds, and v a v a = w a w a = 1 shall be attained. 
Then, inserting the second fundamental tensor cf. section 5.5, this becomes 
equivalent to 

H = \{v a v a + w a w a )k aP = \{g^ + u a vP)k aP . (5.13) 

But this latter relation makes sense without any reference to condition eq. 
(5.12). Therefore, we define eq. (5.13) to be the general relativistic analogue 
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to the mean curvature of a surface. For the case H + ^ H we take their 
arithmetic mean 

H= (H + + H-)/2, (5.14) 
and then we obtain from eqs. (5.11) and (5.13) 

-2HA= l -(k + aP + K,)r^. (5.15) 

This choice can be accepted noting that in the Newtonian limit 

\H + -if" | < \H + + H-\ (5.16) 

anyhow. To compare this with Kelvin's relation we have to relate the right- 
hand side of equation (5.15) to the pressure difference AP at S. To this end 
we investigate the equation of motion for the surface layer. 

5.4 Equation of motion for the surface layer 

The equation of motion, T^. k = 0, contains products of 5 distributions and 
#-step functions, where 9(x) — 1 for x > 0, and 6(x) = else, at points where 
F l j k has a jump discontinuity. These products require special care; cf. [53] 
for a discussion of his point. But defining 6 • 5 = |<5 we obtain, cf. section 
5.5 

AP = A ni n k n = ± (k+p + k-p) r a ? and (5.17) 
A?; 1 = An^Tl = -r% . (5.18) 

From equation (5.18) we see the following: The equation = holds only 
under the additional presumption that the regular, i.e., not delta-like, part 
of has no jump on E. 
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This condition is fulfilled e.g., presuming £ to be such a boundary that 
the regular energy flow does not cross it and the four-velocity is parallel to 
E in both V+ and V-. This we will presume in the following. Then AP 
is indeed the difference of the pressures on both sides, and together with 
equations (5.15) and (5.17) we obtain exactly Kelvin's relation eq. (5.1). 
That means, it is the definition of mean curvature used here that enables us 
to generalize Kelvin's formula to general relativity. Furthermore, 0(M/R) 
of equation (5.10) vanishes. 

Finally we want to discuss the equation = 0. Transvection with u a 
and 6" + u a u^ yields 

u\ a = and u a u ll{a + (In A)\\ a (8° + u a u^ = (5.19) 

respectively. But A is a constant here, and therefore u a is an expansion-free 
geodesic vector field in E. But observe that the u a lines are geodesies in V4 
under additional presumptions only. 

Here, we have only considered a phenomenological theory of surface ten- 
sions, and, of course, a more detailed theory has to include intermolecular 
forces. But on that phenomenological level equations (5.11) and (5.19) to- 
gether with Kelvin's relation eq. (5.1), which has been shown to follow from 
the equation of motion, and A = const as a solely temperature-dependent 
equation of state complete the usual general relativistic Cauchy problem for 
a thermodynamical system by including surface tensions. 

5.5 Notation, distributions and mean curvature 

To make the equations better readable, some conventions and formulas shall 
be given. Let £ a , a = 0,2,3, be coordinates in E and x l , i = 0,1,2,3, 
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those for V4. The embedding E C V4 is performed by functions x l (^ a ) whose 
derivatives 

e\ = dx'/d? = x] a (5.20) 

form a triad field in S. E divides, at least locally, V4 into two connected 
components, V + and VI, and the normal rij, defined by 

= 1 , = (5.21) 

is chosen into the V+ direction, which can be thought being the outer region. 
Possibly V + and VI are endowed with different coordinates x % +) x l _ and met- 
rics gik+ and g ik - , respectively. For this case all subsequent formulas had to 
be indexed with +/— , and only the inner metric of E, its first fundamental 
tensor 

9ap = e^ej g ik (5.22) 

has to be the same in both cases. As usual, we require gik to be C 2 - 
differentiable except for jumps of at E. Covariant derivatives within 
V4 and E will be denoted by ; and || respectively. 

The second fundamental tensor on both sides of E is defined by 

k% = (e^n^Y = (e^m)* (5.23) 

and the difference, A/c a/3 = k + af3 — k~p, Ak = g afS Ak al3 , enters the energy- 
momentum tensor via equation (5.7) and the relation 

r ifc = e^r^, where KT aj3 = g a pAk — Ak aj3 (5.24) 

cf. e.g., [117]. From equation (5.21) and the Lanczos equation (5.24) one 
obtains nVjfc = 0, i.e., indeed the absence of a delta-like energy flow in the 
normal direction. 
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Now take a special coordinate system: x a = £ a , and the x 1 lines are 
geodesies starting from £ into rii direction with natural parameter x 1 . Then 
the line element of V4 reads 

ds 2 = -(dx 1 ) 2 + g aP dx a dx 13 (5.25) 

and the only jumps of F l j k are 

F ±1 - -h ±Y - --n ± 

The most natural definition of 9 ■ 5 is |5 being equivalent to the choice 

V )k = \{^ + ^) on E . (5.26) 

But cf. Dautcourt 1963 [60] for another choice of r* fc with the consequence 
that T* k ^ at E. 

Now the 5 part of the equation T* fe = reads 

AT, 1 ee T^ - = i + fc^) . (5.27) 

Analogously one obtains for the other components 

A^ 1 = -rfp • (5.28) 

Reintroducing the original coordinate system, the left-hand sides of equa- 
tions (5.27) and (5.28) have to be replaced by A Pi = An^Tl and AP a = 
Ariie^Tl, respectively; cf. [151]. Thereby, AP is the difference of the energy 
flows on both sides of E. 

To deduce equation eq. (5.8) we take proper time £° and angular coordi- 
nates ip = £ 2 and tp = £ 3 . Then, inside E, 

ds 2 = - (dZ ) 2 + PW . (5.29) 
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Using equation (5.23) and the exterior Schwarzschild solution one obtains 
k£p the non-vanishing components of which are 

k+ = -M/R 2 (l - 2M/R) 1/2 , k^ 2 = R(l - 2M/R) 1/2 (5.30) 

and 

fc 33 = k 22 sin2 ( 5 - 31 ) 

because of spherical symmetry. To avoid long calculations with the metric 
(5.4) one can proceed as follows. By construction, r 00 = 0, and together with 
equation (5.24) and the spherical symmetry k 22 = k 22 , k 33 = k 33 follows. 
From equation (5.27), equation (5.8) follows then immediately without the 
necessity of determining the actual value of fcjjj,. 

And to be independent of the discussions connected with equation (5.26), 
we deduce equation (5.27) another way. First, independent of surface layers, 
for an arbitrary timelike hypersurface and a coordinate system such that eq. 
(5.25) holds, we have 

kT u = X - ( {3) R + k 2 - k aP k af} ) (5.32) 

where ^R is the curvature scalar within that surface. Now turn to a surface 
layer with h*p ^ k~p. Then equation (5.32) splits into a "+" and a "— " 
equation, having in common solely ^R. Inserting all this into equation 
(5.24), one obtains 

\ (*& + Kp) ^ = y k (*& + k *) - MaP ) 

1 \k+) 2 - (AT) 2 - fc+ + k' k-^] = T+ - T n (5.33) 



2k 



i.e., just equation (5.27). 
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6 The massive scalar field in a closed Fried- 
mann universe 

For the minimally coupled scalar field in Einstein's theory of gravitation we 
look for the space of solutions within the class of closed Friedmann universe 
models. We prove D > 1, where D is the dimension of the set of solutions 
which can be integrated up to t — > oo. D > was conjectured by Page in 
1984 [175]. We discuss concepts like "the probability of the appearance of a 
sufficiently long inflationary phase" and argue that it is primarily a proba- 
bility measure \i in the space V of solutions and not in the space of initial 
conditions, which has to be applied. The measure \i is naturally defined for 
Bianchi-type I cosmological models because V is a compact cube. The prob- 
lems with the closed Friedmann model, which led to controversial claims in 
the literature, will be shown to originate from the fact that V has a compli- 
cated non-compact non-Hausdorff Geroch topology: no natural definition of 
fi can be given. We conclude: the present state of our universe can be ex- 
plained by models of the type discussed, but thereby the anthropic principle 
cannot be fully circumvented. We consider a closed Friedmann cosmological 
model, 

ds 2 = g ij dx i dx j = dt 2 - a 2 {t)[dr 2 + sin 2 r(# 2 + sin 2 ^ 2 )] , (6.1) 

with the cosmic scale factor a(t). We apply Einstein's General Relativity 
Theory and take a minimally coupled scalar field <fi without self-interaction 
as source, i.e., the Lagrangian is 

In this chapter we assume units such that H — c — 1; it! is the scalar curvature, 
G Newton's constant, and m the mass of the scalar field. The variational 
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derivative 5 C/ 5(f) = yields the field equation for the scalar field 

(m 2 + D)0 = , □ = g ij ViVj (6.3) 

where □ denotes the covariant D'Alembertian. The variation of the La- 
grangian with respect to the metric 

6£y/-det gkl /6g ij = (6.4) 

yields Einstein's field equation 

Rij - ^9x3 = ^ GT H , T ij = V,0V,0 - l - (v>V fc - m 2 2 ) 9ij , (6.5) 

with Ricci tensor R^. It is the aim of this chapter, which is based on [215], 
to present some rigorous results about the space V of solutions of eqs. (6.3), 
(6.5) with metric (6.1). This is done in sections 6.1, 6.2 and 6.3. We discuss 
them in the context of cosmology in section 6.4: the probability of the ap- 
pearance of a sufficiently long inflationary phase and the anthropic principle. 

Schrodinger [234] already dealt with the massive scalar field in the closed 
Friedmann universe in 1938, but there the back-reaction of the scalar field on 
the evolution of the cosmic scale factor had been neglected. 35 years later the 
model enjoyed a renewed interest, especially as a semi-classical description 
of quantum effects, cf. e.g. Fulling and Parker [82], [83], Starobinsky [243], 
Barrow and Matzner [14], Gottlober [92] and Hawking and Luttrell [108]. 
One intriguing property - the possibility of a bounce, i.e., a positive local 
minimum of the cosmic scale factor - made it interesting in connection with 
a possible avoidance of a big bang singularity a(t) — > 0. The existence of 
periodic solutions (<j>(t) , a(t)) of eqs. (6.1), (6.3) and (6.5) became clear in 
1984 by independent work of Hawking [106] and Gottlober and Schmidt [96]. 
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The existence of an inflationary phase, which is defined by 



\dh/dt\ <C h 2 , where h = d(\na)/dt 



(6.6) 



is the Hubble parameter, in the cosmic evolution is discussed in many papers 
to that model, cf e.g. Belinsky et al. [20], [21], Gottlober and Miiller [93] and 
Page 1987 [175]. Soon it became clear that the satisfactory results obtained 
for the spatially flat model - the existence of a naturally defined measure in 
the space of solutions and with this measure the very large probability to 
have sufficient inflation - cannot be generalized to the closed model easily. 
We shall turn to that point in section 6.4. 

6.1 Some closed-form approximations 

We consider the system (6.3), (6.5) with metric (6.1). Sometimes, one takes 
it as an additional assumption that depends on the coordinate t only, but 
it holds, cf. Turner [258]: the spatial homogeneity of metric (6.1) already 
implies this property. Proof: For i ^ j we have R^ = = and, therefore 
we get Vi4>V j<fi = 0. This means, locally depends on one coordinate x l 
only. Supposed i ^ 0, then g^Tij — 2T 00 = m 2 2 . The l.h.s. depends on 
t only. For m ^ this is already a contradiction. For m = we have 
Rij = Vj0Vj0, which is a contradiction to the spatial isotropy of R^, q.e.d. 

We always require a(t) > 0, for otherwise the metric (6.1) is degenerated, 
leading to a "big bang". Inserting metric (6.1), eqs. (6.3), (6.5) reduce to 



d 2 <p/dt 2 + 3hd(f)/dt + m 2 <p = 



(6.7) 



and 



3 (h 2 + a' 2 ) = 4tvG [m 2 2 + (d<j)/dt) 2 } . 



(6.8) 
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Eq. (6.8) is the 00-component of eq. (6.5), the other components are a conse- 
quence of it. To get the analogous equations for the spatially flat Friedmann 
model, the l.h.s. of eq. (6.8) has to be replaced by 3h 2 . 

For the massless case m = representing stiff matter, i.e., pressure equals 
energy density, eqs. (6.7), (6.8) can be integrated in closed form. First, 
for the spatially flat model, we get solution eq. (3.4) with r = 1/3, i.e. 
a(t) = t 1//3 ; this means a — 1 in eq. (3.5). Second, for the closed model, let 
ij) = d<f>/dt, then eq. (6.7) implies ipa 3 = const. Inserting this into eq. (6.8), 
we get 

da/dt = ±^L 4 /a 4 - 1 (6.9) 

with a constant L > and < a < L. We get a(t) via the inverted function 
up to a t translation from the following equation 

±t = -LaTcsm(a/L) - -yl - a 2 /L 2 . (6.10) 

The minimally coupled massless scalar field in Einstein's theory is confor- 
mally equivalent to the theory of gravity following from the Lagrangian 
L = R 2 , cf. Bicknell [29], or [245] and [16]. Details of this Bicknell the- 
orem are presented in chapter 8, especially section 8.5. 

For a<Lwe have from eq. (6.10) a ~ t 1 ^ 3 and = 0o + 0i Int. The limit 
L — > oo is, only locally of course, the limit from the closed to the spatially 
flat Friedmann model. Thus we recover, as it must be the case, r = 1/3 given 
above. Eq. (6.10) shows that the function a(t) has a maximum at a — L. 
a(t) > is fulfilled for a t- interval of length At = nL/2 only, and there is 
no bounce. In which range can one except the massless case to be a good 
approximation for the massive case? To this end we perform the following 
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substitutions: 

t — t/e, d — a/e, m = me , = 0. (6-11) 

They do not change the differential equations. Therefore, to get a solution 
a(t) with a maximum a max = e < 1 for m fixed, we can transform to a 
solution with d max = 1, rh = em <C m and apply solutions a (tj with rh = 
as a good approximation. Then a(t) > is fulfilled for a t-interval of length 
At ~ 7ra max /2 only. This is in quite good agreement with the estimate in eq. 
(22) of Page 1984 [175]. But in the contrary to the massless case it holds: 
to each e > there exist bouncing solutions which possess a local maximum 

Let henceforth be m > 0. Then we use ^j4nG/3 ■ instead of as scalar 
field and take units such m — 1. With a dot denoting d/dt we finally get 
from eqs. (6.7), (6.8) 

+ 3/10 + = 0, h 2 + a~ 2 = 2 + 2 . (6.12) 

— > —0 is a Z 2 -gauge transformation, where Z 2 denotes the two-point group. 
Derivating eq. (6.12) we can express and as follows 

= ±v / 2 + 2d 2 + aa/v / 3^2, = ± Vl + a 2 - aa /v 7 ^ . (6.13) 

Inserting eq. (6.13) into eq. (6.12) we get 

d 3 a 

a 2 — = 4d(l + d 2 ) - 3add ± 2aV2 + 2d 2 + ad VI + d 2 - ad . (6.14) 

On the r.h.s. we have " + " if 00 > and " — " otherwise. At all points t, 
where one of the roots becomes zero, " + " and " — " have to be interchanged. 
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To get the temporal behaviour for very large values a ^> 1 but small values 
\h\, we make the ansatz 

a(t) = l/e + eA(t), e>0, e^O. (6.15) 

In lowest order of e we get from eq. (6.14) 



d 3 A 
dt 3 



= ±2 



y/2 + A y/l-A. (6.16) 



An additive constant to A(t) can be absorbed by a redefinition of e, eq. (6.15), 
so we require A(0) = 0. After a suitable translation of t, each solution of eq. 
(6.16) can be represented as 

A(t) = at - t 2 /A + I sin(2t) , \a\<n/4. (6.17) 
8 

a = 7r/4 and a = — 7r/4 represent the same solution, so we have a S* 1 space of 
solutions. Here, S n denotes the n-dimensional sphere. a(t) > is fulfilled for 
A(t) > —l/e 2 , i.e. \t\ < 2/e only. In dependence of the value a, A(t) has one 
or two maxima and, accordingly, null or one minimum. The corresponding 
intervals for a meet at two points, a ~ and |a| ~ 7r/4, where one maximum 
and one horizontal turning point 

d 3 a , 

d = d = 0,-^0 (6.18) 

exist. 

6.2 The qualitative behaviour 

Eqs. (6.12) represent a regular system: at t = we prescribe do, sgn do, 0o 5 
0o fulfilling 

aS + (6-19) 
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and then all higher derivatives can be obtained by differentiating: 



a = 



\J a 2 ((j) 2 + 2 ) — 1 , 4> = —3a<j)/a — <p , 



(6.20) 



the r.h.s. being smooth functions. It follows 



d/dt (h 2 + a' 2 ) = -Qh(d(f)/dt) 2 . 



(6.21) 



6.2.1 Existence of a maximum 

The existence of a local maximum for each solution a(t) already follows from 
the "closed universe recollapse conjecture", but we shall prove it for our 
model as follows: if we start integrating with h > 0; otherwise, t — > — t 
serves to reach that; then h 2 + a~ 2 is a monotonously decreasing function as 
long as h > holds, because hcf) = holds at isolated points t only, cf. eq. 
(6.21). We want to show that after a finite time, h changes its sign giving rise 
to a local maximum of a(t). If this is not the case after a short time, then we 
have after a long time ft«l,d»l, and eqs. (6.15), (6.17) become a good 
approximation to the exact solution. The approximate solution (6.15), (6.17) 
has already shown to possess a local quadratic maximum and this property is 
a stable one within C 2 -perturbations. So the exact solution has a maximum, 
too, q.e.d. 

6.2.2 The space of solutions 

We denote the space of solutions for eqs. (6.1), (6.12) by V and endow it 
with Geroch's topology from 1969, see [86] and subsection 6.3.3 for further 
details. Applying the result of subsection 6.2.1, V can be constructed as 
follows: the set of solutions will not be diminished if we start integrating 
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with do = 0. We prescribe 

/ = Og 1 and g = (1 - a a ) ■ sgn(0 o 0o) , (6.22) 

then all other values are fixed. / and g are restricted by / > 0; \g\ < 3, 
where g — 3 and g = — 3 describe the same solution. Therefore 

V = (K x S X )IQ, (6.23) 

where 7£ denotes the real line, considered as topological space, and Q is 
an equivalence relation defined as follows: Some solutions a(t) have more 
than one, but at most countably many, extremal points, but each extremum 
defines one point in 1Z x S 1 ; these are just the points being Q-equivalent. 

Considering eqs. (6.12), (6.22) in more details, we get the following: for 
\g\ < 1 we have a minimum and for \g\ > 1 a maximum for a(t). 

% > ^ 

holds for < g < 3. a(t) is an even function for g — and g — ±3 only. For 
g — it is a symmetric minimum of a(t) and <p is even, too. For g — ±3 it is 
a symmetric maximum of a(t) and is odd. For \g\ = 1 one has a horizontal 
turning point for a(t). Time reversal t — > —t leads to g — > — p. 

6.2.3 From one extremum to the next 

To prove our result, D > 1, where Z) is the dimension of the set of solutions 
which can be integrated up to t — > oo, we need a better knowledge of the 
equivalence relation Q. To this end we define a map 

p-.TZxS 1 ^TZxS 1 (6.25) 
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as follows: For x = (f, g) we start integrating at t — 0. Let 

ti = min{t|t > 0, d(t) = 0} . (6.26) 

For ti < oo we define 

p(x)=x=(f,g) = (f(t 1 ),g(t 1 )) (6.27) 

and otherwise p is not defined. In words: p maps the initial conditions from 
one extremal point of a(t) to the next one. It holds: p is injective and xQy if 
and only if there exists an integer m such that p m (x) = y. Let V m C TZ x S 1 
be that subspace for which p m is defined. Then the inclusion 

{(/,0)|-1<0<1}cVi (6.28) 

means: a minimum is always followed by a maximum; \g\ < 1 implies |g| > 1. 
V 2 7^ follows from the end of section 6.1. For x = (/, g) we define — x = 
(/, — With this notation it holds, see the end of subsection 6.2.2, 

p m (V m )=V- m = -V m , (6.29) 

and for x G V we have 

p m (-p m (a;)) = -:r and p~ m (V m fl Vl m ) = Vi m . (6.30) 

A horizontal turning point can be continuously deformed to a pair of extrema; 
such points give rise to discontinuities of the function p, but for a suitably 
defined non-constant integer power m the function p m is a continuous one. If 
x G Vi\int(Vi) where int denotes the topological interior, then a(t) possesses 
a horizontal turning point. To elucidate the contents of these sentences we 
give an example: 
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For very small values / and / we may apply eqs. (6.15), (6.17) to calculate 
the function p m . In the approximation used, no more than 3 extremal points 
appear, so we have to consider 

p m , me {-2, -1, 0, 1, 2} (6.31) 

only. / and / approximately coincide, so we concentrate on the function 
g{g). The necessary power m is sketched at the curve. To come from m to 
—m, the curve has to he reflected at the line g = g. Time reversal can be 
achieved, if it is reflected at g = —g. V\ is the interval — 1 < g < 2.8853. The 
jump discontinuity of p at g — 1 and the boundary value g = 2.8853 . . . are 
both connected with horizontal turning points g, g — ±1. The shape of the 
function g(g) for m — 1 can be obtained by calculating the extrema of A(t) 
eq. (6.17) in dependence of a and then applying eq. (6.22). For m^l one 
applies eq. (6.30). 

6.2.4 The periodic solutions 

The periodic solutions a(t) are characterized by the fixed points of some p m , 
m > 2, whereas p itself has no fixed points. The notation p m means, that 
the map p, eqs. (6.25), (6.27), will be applied m times. The existence of 
the fixed points can be proved as follows: we start integrating at t — with 
a symmetric minimum, g — 0, of a(t) and count the number m of zeros of 
4>(t) in the time interval < t < t 2 , where a (£2) is the first local maximum 
of a(t). The number m depends on / = 1/ao and has jump discontinuities 
only at points where ^(£2) = 0. For initial values / = /&, where the number 
m jumps from k to k + 1, the function a(t) is symmetric about t = t 2 . But 
a function, which is symmetric about two different points, is a periodic one. 
We call periodic solutions obtained by this procedure, periodic solutions of 
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the first type, they represent fixed points of the map p 2 . Page 1984 [175] 
gives the numerically obtained result: 

/i = 1/ao with a = 0.76207. . . (6.32) 

That fk exists also for very large values k can be seen as follows: take solution 
eqs. (6.15), (6.17) with a very small value e and insert it into eq. (6.13); we 
roughly get m = 1/eir. 

Fixed points of p 4 , which are not fixed points of p 2 , are called periodic 
solutions of the second type. They can be constructed as follows: we start 
with g = but f ^ fk and count the number n of zeros of <fi(t) in the interval 
< t < £3, where a(ts) > is next local minimum of a(t). For / fk, £3 
is defined by continuity reasons. For f — f l , the number n jumps from I to 
/ + 1, and we have a periodic solution. Numerical evaluations yield 

f 1 = 1/ao , a = 0.74720 . . . (6.33) 

The existence of /' for very large values / is again ensured by solution (6.15), 
(6.17). 

6.2.5 The aperiodic perpetually oscillating solutions 

Now we look at the solutions in the neighbourhood of the periodic ones. Let 
x o — {fk, 0) G V2 be one of the fixed points of p 2 representing a periodic 
solution of the first type; ao(t), the corresponding solution has no horizontal 
turning points. Therefore, the function p is smooth at xq. Let us denote the 
circle with boundary of radius e around x by K(e). By continuity reasons 
there exists an e > such that K{e) C int(V r 2 nl / _ 2 ) and the first two extrema 
of the functions a(t) corresponding to points of K(e) are either maxima or 
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minima. Let Ro — K{ £ ) and for n > 1 

R n := Ro n p 2 (i? n _!) ; i?_ n := R fl p~ 2 (i?!_ n ) . (6.34) 

By assumption, p ±2 is defined in R , so i? m is a well-defined compact set 
with x G int i? m for each integer m. It holds 

R n +i c R n , -R- n = —R n (6.35) 

and 

R n = {x\p 2m (x) G R for m = 0, . . . n} . (6.36) 

Then the set 

oo 

Roo ■= fl R n (6-37) 

?1=0 

is a non-empty compact set with x G -Rqo. In words: R n represents the set 
of all these solutions which have n different maxima. R^ then corresponds 
to the set of solutions possessing infinitely many maxima. For each x G -Roo , 
the corresponding solution a(t) can be integrated up to t — > oo. The last 
statement follows from the fact that the time from one extremum to the next 
is bounded from below by a positive number within the compact set R , i.e., 
an infinite number of extrema can be covered only by an infinite amount of 
time. Analogous statements hold for R-oo = —Roo and t — > — oo. 

Let us fix an integer m > 1. We start integrating at x$ := (/&, 5) G Ro, 
0<5<e. 

S(m) := mzx{5\5 < e, p 2k (x s ) G R for k = 0, . . . m} (6.38) 

exists because of compactness, i.e., xs( m ) £ Rm and there is an integer k{m) < 
m such that 

y{m) := p 2k ^ (x 5{m) ) e5R , (6.39) 
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5Rq being the boundary of Rq. 



j(m) : 



1 , if k{m) < m/2 



(6.40) 







otherwise . 



The sequence 



(y(m), j(m)) C S 1 x Z 2 



(6.41) 



possesses a converging subsequence with 



y(oo) := lim y{m^) , j(oo) := lim j(™i) • 



(6.42) 



Now let us start integrating from y(oo) forward in time for j(oo) = 1 and 
backwards in time for j(oo) = 0. We consider only the case j(oo) = 1, the 
other case will be solved by t — > —t. For each / with mj > 2m we have 
yirrii) e i? m , therefore, y(oo) G i? m for all m, i.e., y(oo) G Roo. y(oo) G 5i? 
and for all n > 1, p 2n (y(oo)) G i?o- 

By continuously diminishing e we get a one-parameter set of solutions 
a £ (t), which can be integrated up to t — > oo. To this end remember that 
one solution a(t) is represented by at most countably many points of Rq. 
Supposed a £ (t) is a periodic function. By construction this solution has a 
symmetric minimum and there exist only countably many such solutions. Let 
the set M of solutions be defined as follows: M = {a(t)\a(0) is a minimum 
parametrized by a point of Rq, a(ti) is the next minimum, and < t < t{\\ 
then we get M = [a m i n , a max ] with a min > 0, a max < oo, and for each e 
and each t > it holds a min < a e (t) < a mSuX . So we have proven: in each 
neighbourhood of the periodic solutions of the first type there exists a set 
of Hausdorff dimension D > 1 of uniformly bounded aperiodic perpetually 
oscillating solutions which can be integrated up to t — > oo. 
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6.3 Problems with die probability measure 

To give concepts like "the probability p of the appearance of a sufficiently 
long inflationary phase" a concrete meaning, we have to define a probability 
measure p in the space V of solutions. Let us suppose we can find a hy- 
persurface H in the space G of initial conditions such that each solution is 
characterized by exactly one point of H . Then V and H are homeomorphic 
and we need not to make a distinction between them. Let us further suppose 
that H C G is defined by a suitably chosen physical quantity ip to take the 
Planckian value. Then we are justified to call H the quantum boundary. By 
construction, H divides G into two connected components; ip < ?/>pi defines 
the classical region. All classical trajectories start their evolution at H and 
remain in the classical region forever. 

Let us remember the situation for the spatially flat Friedmann model 
Belinsky et al. [20] subsection 6.3.1, and for the Bianchi type I model, Lukash 
and Schmidt [145] subsection 6.3.2, before we discuss the closed Friedmann 
model in subsection 6.3.3. 

6.3.1 The spatially flat Friedmann model 

For this case, eq. (6.21) reduces to hh = —3h(j) 2 , i.e., each solution crosses 
the surface h = hp\ exactly once, the only exception is the flat Minkowski 
spacetime h = 0. Reason: For <fi ^ we get h < 0. The corresponding 
physical quantity 

$ = h 2 = 2 + 2 (6.43) 

is the energy density. The space of non-flat spatially flat Friedmann models 
is topologically S 1 , and equipartition of initial conditions gives a natural 
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probability measure there. With this definition it turned out that for m <C 
topi it holds 

p ?s 1 - 8m/m P i , (6.44) 

cf. e.g. Miiller and Schmidt [163], i.e., inflation becomes quite probable. 
If, on the other hand, equipartition is taken at some ho <C h?\m/m-p\ then 
inflation is quite improbable. The total space V of solutions has Geroch 
topology aS 1 , i.e., V = S 1 U {a}, and the space itself is the only neighbour- 
hood around the added point a which corresponds to h = 0, because each 
solution is asymptotically flat for t — > oo. 

6.3.2 The Bianchi-type I model 

With the metric 

the analogue to eq. (6.43) is 

if, = h 2 = <j) 2 + 2 + r 2 + s 2 , (6.46) 

and ft, = ftpi defines a sphere S* 3 in eq. (6.46). Here all solutions cross this 
sphere exactly once, even the flat Minkowski spacetime: it is represented as 
(0 1)-Kasner solution a, so the space of solutions is V = S 3 /Q, where Q 
is a 12-fold cover of S 3 composed of the Z 2 -gauge transformation <fi — > — <fi 
and of the six permutations of the three spatial axes. Eq. (6.46) induces a 
natural probability measure on the space ifj = ip-p\ by equipartition, and the 
equivalence relation Q does not essentially influence this. As in subsection 
6.3.1, a E V has only one neighbourhood: V itself. Up to this exception, 
V is topologically a 3-dimensional cube with boundary. One diagonal line 
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through it represents the Kasner solution = 0. The boundary 5V of V has 
topology S 2 and represents the axially symmetric solutions and one great 
circle of it the isotropic ones. As usual, the solutions with higher symmetry 
form the boundary of the space of solutions, p turns out to be the same as in 
subsection 6.3.1 and how ip eq. (6.46) can be invariantly defined, is discussed 
in [145]. 

6.3.3 The closed Friedmann model 

Now we come to the analogous questions concerning the closed Friedmann 
model. Before defining a measure, one should have a topology in a set. I 
feel it should be a variant of Geroch [86]. The Geroch topology, cf. [210], 
is defined as follows: let Xi = (aj(i), <f>i(t)) be a sequence of solutions and 
x = (a(t), 4>(t)) a further solution. Then Xi — > x in Geroch's topology, if there 
exist suitable gauge and coordinate transformations after which a,i(t) — > a(t) 
and <f>i(t) — > (f>(t) converge uniformly together with all their derivatives in 
the interval t G [— e, e] for some e > 0. Because of the validity of the field 
equations, "with all their derivatives" may be substituted by "with their first 
derivatives" . 

With this definition one gets just the same space V as in subsection 6.2.2, 
eq. (6.23). The existence of aperiodic perpetually oscillating solutions, which 
go right across the region of astrophysical interest subsection 6.2.5 shows 
that for a subset of dimension D > 1 of R x S 1 , Q identifies countably many 
points. All these points lie in a compact neighbourhood of the corresponding 
periodic point (/*,, 0) and possess therefore at least one accumulation point 
z. At these points z, V has a highly non-Euclidean topology. Further, V 
has a non-compact non-Hausdorff topology. So there is no chance to define a 
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probability measure in a natural way and no possibility to define a continuous 
hypersurface in the space of initial conditions which each trajectory crosses 
exactly once. 

6.4 Discussion of the inflationary phase 

Supposed, we had obtained a result of the type: "Each solution a(t) has at 
least one but at most seven local maxima." Then one could define — up to 
a factor 7 = 0(1) — a probability measure. So it is just the existence of 
the perpetually bouncing aperiodic solutions which gives the problems. We 
conclude: it is not a lack of mathematical knowledge but an inherent property 
of the closed Friedmann model which hinders to generalize the convincing 
results obtained for the spatially flat model. So it is no wonder that different 
trials led to controversial results, cf. Belinsky et al. [20], [21] and Page 1987 
[175]. One of these results reads "inflationary and non- inflationary solutions 
have both infinite measure", hence, nothing is clear. 

Let us now discuss the results of Starobinsky [243] and Barrow and 
Matzner [14] concerning the probability of a bounce. They have obtained 
a very low probability to get a bounce, but they used equipartition at some 
ho <C hp\m/mp\. As seen in subsection 6.3.1 for the spatially flat Fried- 
mann model concerning inflation, this low probability does not hinder to get 
a considerable large probability if equipartition is applied at h — hp\. We 
conclude, the probability of bouncing solutions is not a well-defined concept 
up to now. Well-defined is, on the other hand, some type of conditional 
probability. If we suppose that some fixed value a, say 10 28 cm or so, and 
there some fixed value h, say 50 km/sec ■ Mpc or so, appear within the 
cosmic evolution, then the remaining degree of freedom is just the phase of 
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the scalar field, which is the compact set S* 1 as configuration space. But this 
solves not all problems, because the perpetually bouncing solutions discussed 
in subsection 6.2.5 cross this range of astrophysical interest infinitely often. 
Calculating conditional probabilities instead of absolute probabilities, and, if 
this condition is related to our own human existence, then we have already 
applied the anthropic principle. Cf. similar opinions in Singh and Padman- 
abhan [238] concerning the so far proposed explanations of the smallness of 
the cosmological constant. 

The massive scalar field in a closed Friedmann model with Einstein's the- 
ory of gravity cannot explain the long inflationary stage of cosmic evolution 
as an absolutely probable event and so some type of an anthropic principle 
has to be applied; see the well-balanced monograph by Barrow and Tipler 
[17] to this theme. We have discussed the solutions for the minimally coupled 
scalar field, but many results for the conformally coupled one are similar, see 
e.g. Turner and Widrow [259]; this fact in turn can be explained by the ex- 
istence of a conformal transformation relating between them, see [213]. The 
problems in the case of defining a probability measure in the set of not neces- 
sarily spatially flat Friedmann models are also discussed in Madsen and Ellis 
[149]. They conclude that inflation need not to solve the flatness problem. 
The Gibbons-Hawking-Stewart approach [89] gives approximately the same 
probability measure as the equipartition of initial conditions used here. The 
Wheeler-de Witt equation for the massive scalar field in a closed Friedmann 
universe model is also discussed by Calzetta [46]. Possibly, this approach is 
the route out of the problems mentioned here. 



75 



7 The metric in the superspace of Rieman- 
nian metrics 

The space of all Riemannian metrics is infinite-dimensional. Nevertheless a 
great deal of usual Riemannian geometry can be carried over. The superspace 
of all Riemannian metrics shall be endowed with a class of Riemannian met- 
rics; their curvature and invariance properties are discussed. Just one of this 
class has the property to bring the Lagrangian of General Relativity into the 
form of a classical particle's motion. The signature of the superspace metric 
depends on the signature of the original metric in a non-trivial manner, we 
derive the corresponding formula. Our approach, which is based on [217], is 
a local one: the essence is a metric in the space of all symmetric rank-two 
tensors, and then the space becomes a warped product of the real line with 
an Einstein space. 

7.1 The superspace 

Let n > 2, n be the dimension of the basic Riemannian spaces. Let M be 
an n-dimensional differentiable manifold with an atlas x of coordinates x l , 
i = 1, ... The signature s shall be fixed; s is the number of negative 
eigenvalues of the metric. Let V be the space of all Riemannian metrics 
gij(x) in M with signature s, related to the coordinates x l . 

This implies that isometric metrics in M are different points in V in 
general. The V is called superspace, its points are the Riemannian metrics. 
The tangent space in V is the vector space 



T = {hij(x)\x E M, hij = hji} , 



(7.1) 
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the space of all symmetric tensor fields of rank 2. All considerations are local 
ones, so we may have in mind one single fixed coordinate system in M. 

7.2 Coordinates in superspace 

Coordinates should possess one contravariant index, so we need a transfor- 
mation of the type 

y A = ^ gi] {x) (7.2) 

such that the y A are the coordinates for V. To have a defined one-to-one 
correspondence between the index pairs (i, j) and the index A we require 

A=1,...N = n(n + l)/2, (7.3) 

and A = 1, . . . N corresponds to the pairs 

(1,1), (2,2),...(n,n), (1,2), (2, 3), . . . (n - 1, n), (1,3), 

...(n-2,n),...(l,n) (7.4) 

consecutively, (i, j) and (j, i) correspond to the same A. We make the 
ansatz 

b for i 7^ j 

c for i — j (7.5) 
if does not correspond to A 
with certain real numbers b and c to be fixed later and require the usual 
inversion relations 

li Aii tiBn = &i and ^^ Akl = 5 { :6p. (7.6) 

Bracketed indices are to be symmetrized, which is necessary because of sym- 
metry of the metric g^. Inserting ansatz (7.5) into eq. (7.6) gives c 2 = 1, 
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b 2 = 1/2. Changing the sign of b or c only changes the sign of the coordinates, 
so we may put 

c = l, b = l/V2. (7.7) 

The object ji Al ^ is analogous to the Pauli spin matrices relating two spinorial 
indices to one vector index. 

7.3 Metric in superspace 

The metric in the superspace shall be denoted by H AB , it holds 

Hab = H BA (7.8) 

and the transformed metric is 

G ijkl = H AB pt Ai ^ Bkl , H AB = pL Aij piBkiG ijkl . (7.9) 
From eqs. (7.5) and (7.8) it follows that 

Qijkl _ Qjikl _ Qklij j^Q-j 

The inverse to H A b is H AB , and we define 

Gijki = H AB HMjUBki (7-H) 

which has the same symmetries as eq. (7.10). We require G ljkl to be a tensor 
and use only the metric gij(x) to define it. 



Then the ansatz 



Qijkl = Z g i(k g l)j + a ^ l2 j 

Gijki = vg i{k gi)j + pgijg k i (7.13) 
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where v, z, a and (3 are constants, is the most general one to fulfil the 
symmetries eq. (7.10). One should mention that also curvature-dependent 
constants could have heen introduced. 

The requirement that Hab is the inverse to H AB leads via eqs. (7.9) and 
(7.11) to 

Gijki G klmp = <5- m 8f . (7.14) 

The requirement that G ljkl is a tensor can be justified as follows: Let a curve 
y A (t), < t < 1 in V be given, then its length is 



i.e., with eqs. (7.2) and (7.9) 



-r( G,ja f d -wT iu <7Ae) 

A coordinate transformation in M : x l — > ex 1 changes ^ — > s~ 2 gij. 

We now require that a shall not be changed by such a transformation. 
Then a and z are constant real numbers. Inserting eqs. (7.12) and (7.13) 
into eq. (7.14) gives vz — 1, hence z ^ 0. By a constant rescaling we get 

v = z=l (7.17) 

and then eqs. (7.12), (7.13) and (7.14) yield the conditions 

a^— , P = T^~- (7-18) 
n 1 + an 

So we have got a one-parameter set of metrics in V. 
Eq. (7.18) fulfils the following duality relation: with 

f(a) = -«/(! + cm), (7.19) 
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the equation /(/(a)) = a holds for all a ^ —1/n. 

Inserting eq. (7.17) into eqs. (7.12) and (7.13), we finally get for the 
superspace metric 

G m = gi (k g i)j + a g ij g ki and 

Gijki = 9i(k 9i)j + P gij 9ki ■ (7.20) 
It holds: For a = —1/n, the metric Hab is not invertible. 

Indirect proof: G^ kl depends continuously on a, so it must be the case 
with the inverse. But 

lim 

o^ — 1/n 

applied to G^i gives no finite result. Contradiction. 

7.4 Signature of the superspace metric 

Let S be the superspace signature, i.e., the number of negative eigenvalues 
of the superspace metric Hab- S depends on a and s, but a = —1/n is 
excluded. For convenience we define 







0, a > —1/n 

(7.21) 

1, a < —1/n 



From continuity reasons it follows that S is a function of 6 and s: S — 
S(Q, s). If we transform — > —g^ i.e., s — > n — s, then Hab is not changed, 
i.e., 

S{Q,s) = S{Q,n-s). (7.22) 
We transform g^ to diagonal form as follows 

9u = #22 = • • • = g ss = -1, &j = ^ otherwise . (7.23) 
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Let us now give the signature for O = 0. To calculate S(0, s) we may put 
a = and get with eqs. (7.9), (7.12) and (7.17) 

H AB = pLAijpLBkig ik g jl (7.24) 

which is a diagonal matrix. It holds Hu = ... = H nn = 1 and the other 
diagonal components are ±1. A full estimate gives in agreement with eq. 
(7.22) 

S(0, s) = s(n - s) . (7.25) 



Now, we look for the signature for = 1. To calculate S(l,s) we may- 
put a = — 1 and get 

H ab = hmj HBM (g lk g jl - g 13 g kl ) ■ (7.26) 

For A < n < B, Hab = 0, i.e., the matrix Hab is composed of two blocks. 
For A, B < n we get 

f for A = B 
Hab = { (7.27) 
[ 1 for A^B 

a matrix which has the (n — l)-fold eigenvalue 1 and the single eigenvalue 
1 — n. For A, B > n we have the same result as for the case a — 0, i.e., we 
get 5(1, s) = l + s(n- s). 

Result: The signature of the superspace metric is 



S=e + s(n-s). 



(7.28) 



7.5 Supercurvature and superdeterminant 
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7.5 Supercurvature and superdeterminant 

We use exactly the same formulae as for finite-dimensional Riemannian ge- 
ometry to define Christoffel affinities Tg C and Riemann tensor R^ CD in su- 
perspace. We even omit the prefix "super" in the following. Using eq. (7.5) 
we write all equations with indices i,j — l,...n. 

Then each pair of covariant indices i,j corresponds to one contravariant 
index A. The following formulae appear: 



dg 



■'j 



dg 



km 



= -g i{k g m)j (7.29) 



i i 

^ijklmp gi(k gl)(m g P )j _ (xgiJgHmgP)* gi( k g l )( m gP) i (7.30) 

and, surprisingly independent of a we get 

I Mmp = _ 5 (k g l)(m 5 P) _ ( 7 31 ) 

Consequently, also Riemann- and Ricci tensor do not depend on a: 

R klmpij = 1 g l)(m g p)(i 6 J) _ 5 {k gWgjXm ^ (7.33) 

Summing over r = m and s = p we get 

R kiij = 1 ^gijgki _ ngWgj)^ _ ( 7>33 ) 

The Ricci tensor has one eigenvalue 0. Proof: It is not invertible because 

it is proportional to the metric for the degenerated case a = —1/n, cf. eq. 
(7.20) in section 7.3. 

The co-contravariant Ricci tensor reads 



R\\ = G klm pR m ™ = \ (g ij 9ki - nS^Sp) , (7.34) 



82 



7 THE SUPERSPACE OF RIEMANNIAN METRICS 



and the curvature scalar is 

R= -\n(n- l)(n + 2). (7.35) 
8 

The eigenvector to the eigenvalue of the Ricci tensor is gij. All other eigen- 
values equal — n/A, and the corresponding eigenvectors can be parametrized 
by the symmetric traceless matrices, i.e. the multiplicity of the eigenvalue 
— n/A is (n — \){n + 2)/2. This is another view to the well-known split of 
metric perturbations into conformal transformations, i.e., one bulk degree of 
freedom on the one hand, and the remaining volume-preserving degrees of 
freedom on the other. 

We define the superdeterminant 

H = detH AB . (7.36) 

H is a function of g, a and n which becomes zero for a = —1/n, cf. section 
7.3. We use eqs. (7.9) and (7.20) to look in more details for the explicit value 
of H . Though we assume n > 2, here the formal calculation for n — 1 makes 
sense; it leads to 

H = H 11 = G 1111 = g ll g 11 + ag ll g 11 = (1 + a)g- 2 . (7.37) 

Let us return to the general case n > 2. Multiplication of with e gives 
g -> e n g, H AB -> £- 2 H AB and H -> £-< n+1 ^H. 

So we get in an intermediate step 

H = H 1 g- n - 1 (7.38) 

where ifi is the value of H for c/ = 1. Hi depends on a and n only. To 
calculate i?i we put ^ = 5^ and get via H^ = 5y + a, H Ai = for A > n, 
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and Hab = &ab for A, B > n finally 

H x = \ + an. (7.39) 

This result is in agreement with the n = 1-calculation eq. (7.37) and also 
with the fact that a = —1/n gives H = 0. 

7.6 Gravity and quantum cosmology 

Now, we come to the main application: The action for gravity shall be ex- 
pressed by the metric of superspace. The purpose is to explain the math- 
ematical background of quantum cosmology, see [9], [59], [73], [101], [126], 
[148], [158], [179], [182] and [198] for new papers on that topic. 

We start from the metric 

ds 2 = dt 2 - gij dx { dx j (7.40) 

i,j = 1, . . .n with positive definite and x° = t. We define the second 
fundamental form by 

Kii = \9iifi. (7.41) 
The Einstein action for metric (7.40) is 

/ = -/ \ *R Vgd n+1 x (7.42) 

where g = detgij and *R is the (n + l)-dimensional curvature scalar for eq. 
(7.40). Indices at K tj will be shifted with g tj , and K = K\. With eq. (7.41) 
we get 

(K^) fi =(K fi + K 2 )^. (7.43) 
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This divergence can be added to the integrand of eq. (7.42) without changing 
the field equations. It serves to cancel the term K of /. So we get 

1 = I \ ( RijKij - R2 + R ) v^ n+1 * (7-44) 
where R is the n-dimensional curvature scalar for gij. 

Applying eq. (7.20), we now make the ansatz for the kinetic 

W= l - G^ mp K l3 K mp = l - (K ij Kij + aK 2 ) . 

Comparing eq. (7.45) with eq. (7.44) we see that for a = — 1 

holds. 

Surprisingly, this value for a does not depend on n. Because of n > 2 
this value a gives a regular superspace metric. For n — 1, eq. (7.42) is a 
divergence, and a = — 1 does not give an invertible superspace metric. This 
is another form of the result that Einstein gravity does not lead to a local 
field equation in 1+1-dimensional spacetime. 

Using the fx Alj and the notations z A = fi Alj / 2 and v A = dz A /dt we 
get from eqs. (7.45) and (7.46) 

I = J ^(H AB v A v B + R(z A ))^gd n+1 x (7.47) 

i.e., the action has the classical form of kinetic plus potential energy. The 
signature of the metric H AB is S — 1. This can be seen from eqs. (7.21) and 
(7.28). 



energy 

(7.45) 



(7.46) 



7.7 The Wheeler-DeWitt equation 
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7.7 The Wheeler-DeWitt equation 

In eq. (7.47), Einstein gravity is given in a form to allow canonical quantiza- 
tion: The momentum v A is replaced by —id/dz A in units where h — 1, and 
then the Wheeler-DeWitt equation for the world function ip(z A ) appears as 
Hamiltonian constraint in form of a wave equation: 



After early attempts in [7], the Wheeler - DeWitt equation has often been 
discussed, especially for cosmology, see e.g. [30], [89], [100] and [103]. Besides 
curvature, matter fields can be inserted as potential, too. It is remarkable 
that exactly for Lorentz and for Euclidean signatures in eq. (7.40), i.e., 
positive and negative definite respectively, the usual D'Alembert operator 
with S — 1 in eq. (7.48) appears. For other signatures in eq. (7.40), eq. 
(7.48) has at least two timelike axes. 




(7.48) 
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8 Comparison of scalar fields and f(R) for 
cosmology 

Following [209], we generalize the well-known analogies between m 2 2 and 
R + R 2 theories to include the self-interaction A(/> 4 -term for the scalar field. It 
turns out to be the R + R 3 Lagrangian which gives an appropriate model for 
it. Considering a spatially flat Friedmann cosmological model, common and 
different properties of these models are discussed, e.g., by linearizing around 
a ground state the masses of the corresponding spin 0-parts coincide. Then 
we prove a general conformal equivalence theorem between a Lagrangian 
L = L(R), L'L" 7^ 0, and a minimally coupled scalar field in a general 
potential in section 8.5. This theorem was independently deduced by several 
persons, and it is now known as Bicknell theorem [29]. In the final section 8.6, 
which is based on [212] we discuss Ellis' programme, on which length scale 
the Einstein field equation is valid, on microscipic or on cosmic distances? 

8.1 Introduction to scalar fields 

For the gravitational Lagrangian 

L = (R/2 + (3R 2 )/8irG , (8.1) 

where (3 is some free but constant parameter, the value 

R = R crit = -1/4(3 (8.2) 

is the critical value of the curvature scalar, cf. Nariai [169] and Schmidt 
[206], [207]. It is defined by 

OL/dR = . (8.3) 
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In regions where 



R/ -Rcrit < 1 



(8.4) 



holds, we can define 



i> = ln(l - Rj -R cr it) 



(8.5) 



and 



9ij — (1 — R/ -Rcrit) 9ij ■ 



(8.6) 



In units where 8nG = 1 we now obtain from the Lagrangian eq. (8.1) via 
the conformal transformation eq. (8.6) the transformed Lagrangian 



being equivalent to L, cf. Whitt [270]; see [207] for the version of this equiv- 
alence used here. 

For p < 0, i.e., the absence of tachyons in L eq. (8.1), we have massive 
gravitons of mass 



in L, cf. Stelle [248]. For the weak field limit, the potential in eq. (8.7) 
can be simplified to be ?/> 2 /(16 • f3), i.e., we have got a minimally coupled 
scalar field whose mass is also m . The superfluous factor 3/2 in eq. (8.7) 
can be absorbed by a redefinition of ip. Therefore, it is not astonishing, that 
all results concerning the weak field limit for both R + i? 2 -gravity without 
tachyons and Einstein gravity with a minimally coupled massive scalar field 
exactly coincide. Of course, one cannot expect this coincidence to hold for 
the non-linear region, too, but it is interesting to observe which properties 
hold there also. 




(8.7) 




(8.8) 



8.2 The Higgs Geld 
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We give only one example here: we consider a cosmological model of the 
spatially flat Friedmann type, start integrating at the quantum boundary, 
which is obtained by 

RijkiR tjkl 

on the one hand, and T 0Q on the other hand, to have Planckian values, with 
uniformly distributed initial conditions and look whether or not an inflation- 
ary phase of the expansion appears. In both cases we get the following result: 
The probability p to have sufficient inflation is about p — 1 — VXmo/mpi, 
i.e., p = 99.992% if we take m = Kr 5 m P1 from GUT and A = 64, where e A 
is the linear multiplication factor of inflation. Cf. Belinsky et al. [19] for the 
scalar field and Schmidt [207] for R + R 2 , respectively. 

From Quantum field theory, however, instead of the massive scalar field, 
a Higgs field with self-interaction turns out to he a better candidate for 
describing effects of the early universe. One of the advances of the latter 
is its possibility to describe a spontaneous breakdown of symmetry. In the 
following, we try to look for a purely geometric model for this Higgs field 
which is analogous to the above mentioned type where L = R + R 2 modelled 
a massive scalar field. 

8.2 The Higgs field 

For the massive scalar Field we have the mater Lagrangian 

£ m = -(0 ;i ;i -m 2 2 )/2, (8.9) 
and for the Higgs field to be discussed now, 

L\ — + /i 2 2 - A0 4 /12) /2 . (8.10) 
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The ground states are defined by = const, and dL/d(f) = 0. This means 

= for the scalar field, and the three ground states = O — 0, and 

= <j>± = ±VV 2 /A (8.11) 

for the Higgs field. 
The expression 

(9 2 L/90 2 ) 1/2 (8.12) 

represents the effective mass at these points. This gives the value m for the 
scalar field eq. (8.9), so justifying the notation. Further, eq. (8.12) give mass 

1 ii at = and ^/2fi at = <f>± for the Higgs field eq. (8.10). The imaginary 
value of the mass at the ground state = shows the instability met there, 
and in the particle picture, this gives rise to a tachyon. 

To have a vanishing Lagrangian at the ground state (f>± eq. (8.11) we add 
a constant 

A = -3/i 4 /2A (8.13) 
to the Lagrangian eq. (8.10). The final Lagrangian reads 

L = R/2 + L x + A (8.14) 
with L\ eq. (8.10) and A eq. (8.13). 

8.3 The non-linear gravitational Lagrangian 

Preliminarily we direct the attention to the following fact: on the one hand, 
for Lagrangians (8.9), (8.10) and (8.14) the transformation — > —0 is a pure 
gauge transformation, it does not change any invariant or geometric objects. 
On the other hand, 

Rijki — > —Rijki (8.15) 



8.3 The non-linear gravitational Lagrangian 
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or simpler 

R^-R (8.16) 
is a gauge transformation at the linearized level only: taking 

9ik = Vik + eh ik , (8.17) 

where 

77 ifc = diag(l, -1,-1,-1), (8.18) 

then e — > — e implies curvature inversion eq. (8.15) To be strict at the 
linearized level in e. On the other hand, the curvature inversion eq. (8.15), 
and even its simpler version eq. (8.16), fails to hold quadratic in e. This 
corresponds to fact that the e 2 -term in eq. (8.7), which corresponds to the 
ip 3 -term in the development of L in powers of ip, is the first one to break the 
ip — > — ip symmetry in eq. (8.7). 

In fact, the potential is essentially 

(l-e ( - x) ) 2 . (8.19) 

Calculating this to Order = 14, the mapleresult reads 

x 2 - x 3 + — x 4 - - x 5 + — x 6 - — x 7 + 
12 4 360 40 

127 8 17 9 73 10 31 n 
■ x 8 - — — x 9 + — — x 10 - — — — x 11 + 



20160 12096 259200 604800 
2047 1 

x 12 — x 13 + 0(x 14 ) (8.20) 

239500800 760320 K ' y 1 

Now, let us introduce the general non-linear Lagrangian L = L(R) which 
we at the moment only assume to be an analytical function of R. The ground 
states are defined by R = const., i.e., 

L'Rik - g ik L/2 = . (8.21) 

Here, V = dL/dR. 
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8.3.1 Calculation of the ground states 

From eq. (8.21) one immediately sees that dL/dR = defines critical values 
of the curvature scalar. For these values R = R ciit it holds: For L(R crit ) ^ 
no such ground state exists, and for L(R CT i t ) = 0, we have only one equation 
R = R cr it to be solved with 10 arbitrary functions g^. We call these ground 
states degenerated ones. For L = R 2 , R cr[t = 0, this has been discussed by 
Buchdahl [43]. Now, let us concentrate on the case dL/dR ^ 0. Then Ry is 
proportional to with a constant proportionality factor, i.e., each ground 
state is an Einstein space 

R ij = Rg ij /4, (8.22) 

with a prescribed constant value R. Inserting eq. (8.22) into eq. (8.21) we 
get as condition for ground states 

RL' = 2L. (8.23) 

As an example, let L be a third order polynomial 

L = A + R/2 + (3R 2 + XR 3 /12 . (8.24) 

We consider only Lagrangians with a positive linear term as we wish to 
reestablish Einstein gravity in the A — > weak field limit, and (3 < to 
exclude tachyons there. 

We now solve eq. (8.23) for the Lagrangian eq. (8.24). For A = we have, 
independently of /3, the only ground state R = — 4A. It is a degenerated one 
if and only if (3 A = 1/16. That implies that for usual R + R 2 gravity eq. 
(8.1), i.e. A = A = 0, we get R — as the only ground state; it is a 
non-degenerated one. 
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Now, let A 7^ and A = 0. To get non-trivial ground states we need the 
additional assumption A > 0. Then, besides R = 0, the ground states are 

R = R ± = ±^6/\ (8.25) 

being quite analogous to the ground states eq. (8.11) of the Higgs field 
eq. (8.10). The ground state R = is not degenerated. Of course, this 
statement is independent of A and holds true, as one knows, for A = 0. To 
exclude tachyons, we require (5 < 0, then R_ is not degenerated and R + 
is degenerated if and only if (5 — — \JQ/\- The case AA ^ will not be 
considered here. 

8.3.2 Definition of the masses 

For the usual R + R 2 theory eq. (8.1), the mass is calculated at the level 
of the linearized theory. Then the equivalence to the Einstein field equation 
with a scalar field applies, and we use eqs. (8.12) and (8.2) to calculate 

m = ( J Rcrit/3) 1/2 = (-12/3)- 1 / 2 . (8.26) 

Thus, we recover the value eq. (8.8). But how to define the graviton's masses 
for the Lagrangian eq. (8.24)? To give such a definition a profound meaning 
one should do the following: linearize the full vacuum field equation around 
the ground state, preferably de Sitter- or anti-de Sitter spacetime, respec- 
tively, decompose its solutions with respect to a suitably chosen orthonormal 
system, which is a kind of higher spherical harmonics, and look for the prop- 
erties of its single modes. For L eq. (8.1) this procedure just gave m . 

A little less complicated way to look at this mass problem is to consider 
a spatially flat Friedmann cosmological model and to calculate the frequency 
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with which the scale factor oscillates around the ground state, from which 
the mass m turned out to be the graviton's mass for L eq. (8.1), too. 

Keep in mind, 1. that all things concerning a linearization around flat 
vacuous spacetime do not depend on the parameter A neither for the Higgs 
field nor for the L(R) model, and 2. that a field redefinition R — > R* + R± is 
not possible like — > 0* + <j>± because curvature remains absolutely present. 

8.4 The cosmological model 

Now we take as Lagrangian eq. (8.24) and as line element 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ) . (8.27) 
The dot denotes d/dt and h — a/a. We have 

R = -6h-12h 2 , (8.28) 
and the field equation will be obtained as follows. 

8.4.1 The field equation 

For L = L(R) the variation 

5(Lv^)/^ = (8.29) 

gives with L' = dL/dR the following fourth-order gravitational field equation 

L'Rij - gijL/2 + 9ij UL' - L' ;ij = , (8.30) 

cf. e.g., Novotny (1985) [172]; see also the ideas presented by Kerner [125]. 
It holds 

L '-,ij = L"R-ij + U" R-iRj ■ (8.31) 
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With eq. (8.31), the trace of eq. (8.30) reads 

L'R -2L + 3L"DR + 3L"' R. k R k = , (8.32) 
i.e., with L eq. (8.24) 

3A 

-2A - R/2 + \R 3 /12 + 6/3DR + —(RVR + R. k R} k ) = . (8.33) 

Inserting eqs. (8.27), (8.28) and (8.31) into the 00-component of eq. (8.30) 
we get the equation 

= h 2 /2 - A/6 - Q(3(2hh - h 2 + Qh 2 h) 
+3X(h + 2h 2 ) (Qhh + 19h 2 h - 2h 2 - 2h A ) . (8.34) 

The remaining components are a consequence of this one. 
8.4.2 The masses 

Linearizing the trace equation (8.32) around the flat spacetime, hence A = 0, 
gives independently of A of course, R = 12/3Di?, and the oscillations around 
the flat spacetime indeed correspond to a mass m = (— 12/?) -1 / 2 . This once 
again confirms the evaluation eq. (8.8). 

Now, let us linearize around the ground states eq. (8.25) by inserting 
A = and R = ±^/o/X + Z into eq. (8.32). It gives 

Z = ^—6/3 =F y / 27A/2^j UZ , (8.35) 

and, correspondingly, comparing eq. (8.35) with the equation (□ + m±)Z = 
0, we get 

m± = (6/3 ± V / 27A/J) . (8.36) 
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For < -V% m± is imaginary, and its absolute value differs by a factor 
\/2 from m . This is quite analogous to the A0 4 -theory, cf. section 8.2. 
Therefore, we concentrate on discussing this range of parameters. 

For the ground state for A ^ 0, A = we get with R = — 4A + Z just 
Z = 120OZ, i.e., mass m eq. (8.8) just as in the case A = A = 0. 

Let us generalize this estimate to L — L(R); according to eq. (8.23), 
R = Rq = const, is a ground state if 

L'(R )R = 2L(R ) (8.37) 

holds. It is degenerated if L'(R ) = 0. Now, we linearize around R = Rq: 
R = Ro + Z . For L"(Rq) = 0, only Z = solves the linearized equation, and 
R = R is a singular solution. For L"(R ) ^0 we get the mass 

m = (i? /3 - L'(R )/3L"(R )) 1/2 (8.38) 

meaning the absence of tachyons for real values m. Eq. (8.38) is the analogue 
to eq. (8.12) for the general Lagrangian L(R). 

8.4.3 The Friedmann model 

Here we only consider the spatially flat Friedmann model eq. (8.27). There- 
fore, we can discuss only de Sitter stages with R < 0, especially the ground 
state R + eq. (8.25) representing an anti-de Sitter spacetime does not enter 
our discussion, but i?_ does. 

Now, let A = 0. Solutions of eq. (8.34) with constant values h are h = 
representing flat spacetime and in the case that A > also 
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representing the de Sitter spacetime. These are the non-degenerated ground 
states R = and R = itL = — y/E/X, respectively. Eq. (8.34) can be written 

as 

= h 2 (l - 24A/i 4 )/2 + hh (l/m 2 + 18X(h + 2h 2 )) 
-QXh 3 + h 2 (A5Xh 2 - l/2m 2 ) + 3h 2 h(l/m 2 + 3QXh 2 ) . (8.40) 

First, let us consider the singular curve defined by the vanishing of the coef- 
ficient of h in eq. (8.40) in the h — /i-phase plane. It is, besides h = 0, the 
curve 

h = -2h 2 - l/18Amg (8.41) 

i.e., just the curve 

R = l/3Xm 2 = -A(3/X (8.42) 

which is defined by L" = 0, cf. eq. (8.32). This value coincides with R + if 
(3 = - fix/8, this value we need not discuss here. Points of the curve eq. 
(8.41) fulfil eq. (8.40) for 



h = ±l/ 18Xm 3 V3 yjl - l/18Amg (8.43) 

only, which is not real because of A TOq. Therefore, the space of solutions 
is composed of at least two connected components. 

Second, for h — we have h = or 

h=-l/12Xm 2 . (8.44) 

From the field equation we get: h — h — implies hh > 0, i.e. h does not 
change its sign. We know such a behaviour already from the calculations in 
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[163], where the same model with A = was discussed. In a neighbourhood 
of eq. (8.44) we can make the ansatz 

oo 

h = -t/\2\ml + £ a n t n (8-45) 

n=2 

which has solutions with arbitrary values 02- This means: one can change 
from expansion to subsequent recontraction, but only through the "eye of 
a needle" eq. (8.44). On the other hand, a local minimum of the scale 
factor never appears. Further, eq. (8.44) does not belong to the connected 
component of flat spacetime. 

But we are especially interested in the latter one, and therefore, we restrict 
to the subset h > h( eq. (8.41)) and need only to discuss expanding solutions 
h > 0. Inserting h — 0, 

h = h{24Xh 4 - l)/{2/ml + 72\h 2 ) (8.46) 

turns out, i.e., h > for h > 1/ \ / 2A\ only. All other points in the h — h phase 
plane are regular ones, and one can write dh/dh = h/h — F(h, h) which can 
be calculated by eq. (8.40). 

For a concrete discussion let A lO 2 /^ and m = 10 _5 m P1 . Then both 
conditions (3 — V\ and \R_\ < /p 2 are fulfilled. Now the qualitative be- 
haviour of the solutions can be summarized: There exist two special solutions 
which approximate the ground state R_ for t — > — 00. All other solutions 
have a past singularity h — > 00. Two other special solutions approximate the 
ground state R_ for t — > +00. Further solutions have a future singularity 
h — > 00, and all other solutions have a power-like behaviour for t — > 00, a(t) 
oscillates around the classical dust model a(t) ~ t 2 ^ 3 . But if we restrict the 
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initial conditions to lie in a small neighbourhood of the unstable ground state 
only one of the following three cases appears: 

1. Immediately one goes with increasing values h to a singularity. 

2. As a special case: one goes back to the de Sitter stage 

3. The only interesting one: One starts with a finite /pi-valued inflationary 
era, goes over to a GUT-valued second inflation and ends with a power-like 
Friedmann behaviour. 

In the last case to be considered here, let A = 0, A > and (3 < 0. The 
analogue to eq. (8.40) then reads 

= h 2 /2 - A/6 + (2hh -h 2 + 6h 2 h)/2m 2 . (8.47) 

Here, always h ^ holds, we consider only expanding solutions h > 0. For 
h = we have 

h = {Am 2 /3 - m 2 h 2 )/2h . (8.48) 
For h = we have h > ml/6 and 

h = (A/3 + h 2 /m 2 y/ 2 (l + Qh/mlY 1 ' 2 . (8.49) 

Using the methods of [163], where the case A = has been discussed, we 
obtain the following result: All solutions approach the de Sitter phase h 2 = 
A/3 as t — > oo. There exists one special solution approaching h = — ml/ 6 for 
h — > oo, and all solutions have a past singularity h — > oo. For a sufficiently 
small value A we have again two different inflationary eras in most of all 
models. 
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8.5 The generalized Bicknell theorem 

In this section we derive a general equivalence theorem between a non-linear 
Lagrangian L(R) and a minimally coupled scalar field with a general po- 
tential with Einstein's theory. Instead of <j) we take 

^ = ^2/3 0. (8.50) 

This is done to avoid square roots in the exponents. Then the Lagrangian 
for the scalar field reads 

L — R/2 — 3S%^/4 + Vty) . (8.51) 

At ground states ip = ipo, defined by dV/dt/j = the effective mass is 

m = ^2/3^d 2 V/dip 2 , (8.52) 

cf. eqs. (8.12) and (8.50). The variation = 5L/5tp gives 

= dV/dip + 3 g ij ViVj i)/2 (8.53) 

and Einstein's equation is 

En = KTij (8.54) 

with 

KTij = 3^/2 + gij (y(if>) - l~g ab tatb) • (8.55) 

Now, let 

~9i3 = ^9i3- (8-56) 

The conformal transformation eq. (8.56) shall be inserted into eqs. (8.53), 
(8.54) and eq. (8.55). One obtains from eq. (8.53) with 

t k : 

Urj} + = -2(e^dV/dif;)/3 (8.57) 
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and from eqs. (8.54), (8.55) 

Eij = ^ + W-j + 9ij (e*Vty) -Uil>- . (8.58) 
Its trace reads 

-R = Ae^Vty) - 3D-0 - 3^ ;a . (8.59) 
Comparing with eq. (8.57) one obtains 

R = R{$) = -2e-^d (e 2 ^V(^)) /dip . (8.60) 

Now, let us presume OR/ dip ^ 0, then eq. (8.60) can be inverted as 

i> = F(R) . (8.61) 

In the last step, eq. (8.61) shall be inserted into eqs. (8.57), (8.58), (8.59). 
Because of 

F(R). ij = OF /OR ■ R. yij + d 2 F/dR 2 ■ R ;i R. yj (8.62) 

and dF/dR ^ 0, eq. (8.58) is a fourth-order equation for the metric g^. We 
try to find a Lagrangian L = L(R) such that the equation 5Ly^g/5g lj = 
becomes just eq. (8.58). For V = dL/dR ^ 0, eq. (8.30) can be solved to 
be 

Eij = -QijR/2 + gij L/2L' - gijDL'/L' - L'.^/L' . (8.63) 

We compare the coefficients of the R ;i j terms in eqs. (8.58) and (8.63), this 
gives 

dF/dR = L"/U , hence 

L(R) =fx [ R e F{x) dx + A (8.64) 
JRo 

with suitable constants A , /i, and R , \x ^ 0. We fix them as follows: We 
are interested in a neighbourhood of R — R and require L'(R ) = 1/2. 
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Otherwise L should be multiplied by a constant factor. Further, a constant 
translation of tp can be used to obtain F(R ) = 0, hence fj, — 1/2, L(R ) = 
A , and 

L'(Rq) = 0F/dR(R )/2 ^ . (8.65) 

With eq. (8.64) being fulfilled, the traceless parts of eqs. (8.58) and (8.64) 
identically coincide. Furthermore, we have 

UL'/L' = DF + F ;i F ;i (8.66) 

and it suffices to test the validity of the relation 

e F V(F(R)) = -R/2 + L/2L' . (8.67) 

It holds 

2L' = e F , i.e., 
e 2F V(F(R)) = L — Re F /2 . (8.68) 

At R = R , this relation reads V(0) = A - R /2. Applying d/dR to eq. 
(8.68) gives just eq. (8.57), and, by the way, V'(0) = Ro/2 — 2A . In sum, 

L(R) = V(0) + R /2 + [ R e F{x) dx/2 , (8.69) 

JRa 

where F(x) is defined via F(R ) = 0, 

■if) = F (-2e-^9(e 2 ^T/(V'))/9V') . (8.70) 

Now, let us go the other direction: Let L = L(R) be given such that at 
R = R , L'L" ^0. By a constant change of L let L'(R ) = 1/2. Define 
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A = L(Rq), ip = F(R) = ln(2L'(i?)) and consider the inverted function 
R = F _1 (^). Then 

V(4>) = (A - Ro/2)e~ 2lp - e" 2 ^ f e x F~ l (x)dx/2 (8.71) 

Jo 

is the potential ensuring the above mentioned conformal equivalence. This 
procedure is possible at all i?-intervals where V L" ^ holds. For analytical 
functions L(R), this inequality can be violated for discrete values R only or 
one has simply a linear function L(R) being Einstein gravity with A-term. 

Eq. (8.71) is given here in the form published in [209]. Later it turned 
out that this integral can be evaluated in closed form as follows: 

Vt$) = L{F- l m ■ - \ F-\i,) ■ e-+ . (8.72) 

Examples: 1. Let L = A + R 2 , R = 1/4, then AR = and 

V(ip) = Ae~ 2,/, -l/16. (8.73) 

For A = 0, this is proven in Bicknell [29] and Starobinsky and Schmidt [245]. 

2. Let L = A + R/2 + f3R 2 + XR 3 /12, Rq = 0, hence f3 ^ is necessary. 
We get 

- 1 = 4(3 R + XR 2 /2 and 
V(4>) = Ae~ 2 ^ + 

2(3\- l e~ 2lp (f-l- 16/? 2 (3A)- 1 ((1 + A(e^ - /)/8/3 2 ) 3/2 - 1)) . (8.74) 
The limit A — > in eq. (8.74) is possible and leads to 

V{$) = Ae- 2 ^ - (1 - e^) 2 /16p , (8.75) 
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so we get for A = again the potential from eq. (8.7). 

Now, let R be a non-degenerated ground state, hence 

L(R) =A + (R- Rq)/2 + L"(Ro)(R - R ) 2 /2 + ... (8.76) 

with £"(#0) + and A = Rq/A, cf. subsection 8.3.1. Using eq. (8.71) we 
get V'(0) = and 

V"(0)=R /2-l/(AL"(R o )). (8.77) 

Inserting this into eq. (8.52) we exactly reproduce eq. (8.38). This fact once 
again confirms the estimate eq. (8.38) and, moreover, shows it to be a true 
analogue to eq. (8.12). To understand this coincidence one should note that 
at ground states, the conformal factor becomes a constant = 1. 

8.6 On Ellis' programme within homogeneous world 
models 

For the non-tachyonic curvature squared action we show that the expanding 
Bianchi-type I models tend to the dust-filled Einstein-de Sitter model for 
t tending to infinity if the metric is averaged over the typical oscillation 
period. Applying a conformal equivalence between curvature squared action 
and a minimally coupled scalar field, which holds for all dimensions > 2, 
the problem is solved by discussing a massive scalar field in an anisotropic 
cosmological model. 

8.6.1 Ellis' programme 

Ellis [70] has asked in 1984 on which length scale the Einstein field equation 
is valid. On laboratory or on cosmic distances? Here we extend this question 
to cover microscopic scales, too. This is a real problem as one knows: 
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1. An averaging procedure does not commute with a non-linear differential 
operator as the Einstein tensor is. 

2. Einstein's theory is well tested at large distances 3> 1 cm. 

3. The ultraviolet divergencies of Einstein gravity can be removed by adding 
curvature squared terms to the action, see Weinberg [267]; this is a micro- 
scopic phenomenon and forces to prefer a curvature squared action already 
on the level of classical field theory, as we are concerned here. 

Now, we propose a synthesis of 1., 2., and 3. as follows: microscopically, 
we take 



the minus sign before the i? 2 -term shows that we consider the non-tachyonic 
case only. By an averaging procedure we get Einstein gravity on large scales 
^> I PS 10~ 28 cm. As we are dealing with homogeneous cosmological mod- 
els, the average is taken over the time coordinate only: We will show: the 
curvature squared contribution represents effectively dust in the asymptotic 
region t — > oo. For spatially flat Friedmann models this was already proven in 
Miiller and Schmidt [163], here we generalize to models with less symmetry. 
In the present paper we consider only the vacuum case 



and therefore, we interpret the effectively obtained dust as invisible gravitat- 
ing matter necessary to get a spatially flat universe. 

We conjecture that additional matter contributions, that means usual 
dust plus radiation, do not alter the result qualitatively. The general ex- 
pectation for Lagrangian eq. (8.78) is the following: starting at die Planck 



L g = (R/2-l 2 R 2 )/8nG 





(8.79) 
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era 



RijkiR ijM ~ 10 131 cm" 4 , 



(8.80) 



the i? 2 -term is dominant, the inflationary de Sitter phase is an attractor; its 
appearance becomes very probable, see [207] and [245]; the i?-term in eq. 
(8.78) yields a parametric decay of the value of R and one turns over to the 
region t — > oo where R 2 gives just dust in the mean. 

In Carfora and Marzuori [50] another approach to Ellis' programme was 
initiated: a smoothing out of die spatially closed 3-geometry in the direction 
of die spatial Ricci-tensor leads to different values of mean mass density 
before and after the smoothing out procedure. 

8.6.2 The massive scalar field in a Bianchi-type I model 

Consider a minimally coupled scalar field in a potential V((f>), where 87rG = 
1 is assumed, 



We suppose V(<t>) to be a C 3 -function and V(0) = to be a local quadratic 
minimum of the potential V, i.e., V(0) = 0, V"(0) =m 2 ,m> 0. To describe 
the asymptotic behaviour — > as t — > oo it suffices to use V{4>) = m 2 2 /2 
because the higher order terms do not affect it. This statement can be proved 
as follows: For each e > 0, e < m, there exists a O > such that for all 
with |0| < 0o it holds 



Then all further development is enclosed by inequalities with m ± e, e — > 0. 



L — R/2 — l -4>^ + V{4>) . 



(8.81) 



(m - e) 2 2 < 21/(0) < (m + e) 2 <p 2 . 



(8.82) 
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Here, we concentrate on a cosmological model of Bianchi-type I. In Misner 
parametrization it can be written as 

The functions 0, a, r, s depend on t only, a = da/dt = h is the Hubble 
parameter and n = u/h, u = (r 2 + s 2 ) 1 / 2 the anisotropy parameters. It holds 
< i] < 1, i] = represents the isotropic model, and rj = 1 gives = 0, 
therefore, we consider only the case < rj < 1 in the following, h = 
is possible for Minkowski spacetime only, and the time arrow is defined by 
h > 0, i.e., we restrict to expanding solutions. 

In Lukash and Schmidt [145] it was shown that all such solutions tend 
to Minkowski spacetime for t — > oo. More detailed: At t — 0, let a — r — 
s = by a coordinate transformation. Then for each prescribed quadruple 
(r, s, 0, 0) at t — the integration of the relevant system up to t — > oo is 
possible with r, s, 0, 0, /i, h tending to zero in that limit. The relevant 
equations are 

(0 2 + 2 + m 2 ) 1 / 2 = h (8.84) 
+ 3/10 + = (8.85) 
f = C r e~ 3a , s = C s e~ 3a . (8.86) 

We put the mass m — 1; C r , C s are constants. Now, we consider this 
limit in more details. Up to now the following is known: For the isotropic 
models u = 0, e a = a one knows that the asymptotic behaviour is given 
by oscillations around a ~ t 2 / 3 i.e., we get the Einstein-de Sitter model 
in the mean, and the effective equation of state is that of dust, cf. e.g. 
Starobinsky [243]. In Gottlober 1977 [92] this is generalized to a special class 
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of inhomogeneous, nearly isotropic models. After averaging over space and 
oscillation period one gets as effective equation of state 



i.e., also dust in the limit t — > oo. 

The question for the Bianchi-type I model is now: does one get the energy- 
momentum tensor of an ideal fluid in the mean, or will there appear a strongly 
anisotropic pressure? The first question to be answered is about the averaging 
procedure: From eq. (8.85) and the fact that h — > as t — > oo one gets a 
fixed oscillation period: let t — t n be the n-th local extremum of (pit), then 



Therefore, we average the metric about this period of time. Let us denote 
r(t n ) by r n , ... For a fixed finite time it is ambiguous how to perform this 
average, but in the limit t — > oo one can describe the system by adiabatic 
or parametric deviations from pure ~ cost-oscillations, and the averaging 
procedure consists of constructing monotonous smooth curves a(t) fulfilling 
a{t n ) = a n , . . . , and the effective energy-momentum tensor is obtained from 
eq. (8.78) with a — a(t), ... using the Einstein equation. 

Derivating eq. (8.84) one gets with eq. (8.85) 




(8.87) 



lim t 



■n+l 



(8.88) 



> h = 30 2 - 3h 2 > -3h 2 



(8.89) 



and therefore, 



K > K+i > K - 3nhl 



(8.90) 



Let 



h > 



1 



(8.91) 



3n(n + no) 



8.6 On Ellis' programme within homogeneous world models 



109 



be valid for one value n — ri\ then, by induction, this holds true also for all 
larger values n, i.e., 

ft w £ mhj (8 - 92) 

for some to, let to = subsequently. Derivating eq. (8.86) we get with eq. 
(8.84) 

rj = _ 3?7 (i _ r?) h . . (8.93) 

The last factor can be substituted by its mean value, which equals 1/2. From 
eq. (8.92) one can see that for initial conditions < rj < 1 as we met here, eq. 
(8.93) leads to r] — > as t — > oo. Knowing this, we can perform a stronger 
estimate for h, because in the mean, h = —3h 2 /2, i.e., in the leading order 
we have h(t) = 2/3t and, using eq. (8.93), we get 

fi = -ri/t, V = Vo/t. (8.94) 

Remember: the inflationary era diminishes n exponentially with increasing 
time, but here we consider only the asymptotic region. In sum we get for 
metric (8.83) 

e a = t 2/3 , r = -C r /t , s = -C 8 /t . (8.95) 

Inserting eq. (8.95) into eq. (8.83) we get, via the Einstein equation, a 
diagonal energy- momentum tensor with p = Too — ^/3t 2 and pressure of the 
order (C 2 + C 2 )/t 4 . The effective equation of state is p ~ p/t 2 ~ p/a 3 , i.e., 
also dust in the limit t — > oo as eq. (8.87). From the details of the paper 
[145] one can see that for almost all models the end of the inflationary stage 
is just the beginning of the Einstein-de Sitter stage. A similar result occurs 
if not the metric, but p and p will be averaged, cf. Gottlober 



Result: For the minimally coupled scalar field in a potential V which 
has a single quadratic minimum at V = all expanding Bianchi-type I 
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cosmological solutions tend to the Einstein-de Sitter model for t — > oo if the 
metric is averaged over the oscillation period. 

8.6.3 The generalized equivalence 

Some types of a conformal equivalence theorem between fourth-order gravity 
and minimally coupled scalar fields are obtained in [209], [211]. This theo- 
rem was also independently obtained by Ferraris [76], see also Jakubiec and 
Kijowski [120], Goenner [90]. All of them are restricted to 4-dimensional 
spacetimes. On the other hand, both i? 2 -terms and scalar fields have been 
discussed for higher-dimensional spacetimes, cf. e.g. Ishihara [116]. 

Therefore, it is worth mentioning that this conformal equivalence theorem 
can be formulated for arbitrary dimensions n > 2: Let 

C = R/2- + V(<f>) (8.96) 

and 

5. . = e^ga , A = - £ (8.97) 

yj yj ' (n-l)(n-2) V ; 

be the conformally transformed metric. Then the solutions of the variation 

of eq. (8.96) are transformed by eq. (8.97) to the solutions of the variation 

of L = L(R), where 



is supposed to be locally, i.e. near R = Rq, invertible as 

<f> = F(R), F(R )=0, F'(R )^0 

C(R) = Ir + V(0) + I l R e F ^^dx . (8.99) 
2 2 Jr 
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The inverse direction is possible provided L'(R) ■ L"(R) ^ 0, cf. [209], [211] 
for details with n — 4. Other cosmological models with higher dimensions 
are in [32], [119], and [190]. Two-dimensional models are discussed in [154] 
and [229]. 

8.6.4 The fourth-order gravity model 

Now we come to the question posed in the introduction: Let L(R) be a 
C 3 -function fulfilling L(0) = 0, L'(0) • L"(0) < 0. Then we can write 

C(R) = - + [3R 2 + 0(R 3 ), (3<0. (8.100) 

We consider the Bianchi-type I vacuum solutions which start in a neigh- 
bourhood of the Minkowski spacetime and ask for the behaviour as t — > oo. 
Applying the equivalence theorem cited in subsection 8.6.3 we arrive at the 
models discussed in subsection 8.6.2, and this is applicable for \R \ being 
small enough. The conformal factor depends on t only, and therefore, the 
space of Bianchi-type I models will not be leaved, and we can formulate the 
following: In a neighbourhood of Minkowski spacetime, all Bianchi-type I 
models which represent a stationary point of the action eq. (8.100), can be 
integrated up to t — > oo or — oo, let it be +oo. One singular solution is the 
Kasner solution and all other solutions undergo isotropization and have an 
averaged equation of state p = for t — > oo. 

Let us conclude: Anderson [6] discussed the possibility that curvature 
squared terms give an effective contribution to the energy density of a Fried- 
mann model. This could explain the discrepancy between the observed mean 
mass density of about 1/10 the critical one and the predicted one, from infla- 
tionary cosmology, nearly the critical one. Here, we have shown that also for 
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a large initial anisotropy the oscillating curvature squared contributions give 
just dust in the mean and not an equally large anisotropic pressure as one 
could have expected. The next step would be to look for a generalization of 
this fact to inhomogeneous cosmological models. Further fourth-order grav- 
ity models can be seen in [11], [48], [49], [61], [74], [110], [111], [115], [124], 
[128], [139], [193], [214] and [246]. Another approach to averaging procedures, 
i.e., to Ellis' programme, is due to Zalaletdinov [275]. 
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9 Friedmann models with de Sitter 
and power-law inflation 

Following [216], we consider the spatially flat Friedmann model 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ) . (9.1) 

For a ~ t p , especially, if p > 1, this is called power-law inflation. For the 
Lagrangian 

L = R m with p=-(m-l)(2m-l)/(m-2) (9.2) 

power-law inflation is an exact solution, as it is for Einstein gravity with a 
minimally coupled scalar field $ in an exponential potential 

V($) = exp(ji$) (9.3) 

and also for the higher-dimensional Einstein equation with a special Kaluza- 
Klein ansatz. The synchronized coordinates are not adapted to allow a 
closed-form solution, so we write 

ds 2 = a 2 (Q 2 (a) da 2 - dx 2 - dy 2 - dz 2 ) . (9.4) 

The general solutions reads 

Q(a) = (a b + CY' b (9.5) 

with free integration constant C, where C = gives exact power-law inflation, 
and m-dependent values b and /: 

/ = -2 + l/p, 6=(4m-5)/(m-l). (9.6) 

Further, special solutions for the closed and open Friedmann model are found. 
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The de Sitter spacetime 



ds 2 = dt 2 - e 2Ht {dx 2 + dy 2 + dz 2 ), H^O 



(9.7) 



is the spacetime being mainly discussed to represent the inflationary phase 
of cosmic evolution. However, a spacetime defined by 



enjoys increasing interest for these discussions, too. Especially, eq. (9.8) with 
p > 1, t > is called power-law inflation; and with p < 0, t < it is called 
polar inflation. 

We summarize some differential-geometrical properties of both de Sitter 
and power-law/polar inflation in section 9.1 and show, from which kind of 
scale-invariant field equations they arise in section 9.2; we give the complete 
set of solutions for the spatially flat and special solutions for the closed and 
open Friedmann models in closed form for field equations following from the 
Lagrangian L = R m in section 9.3, and discuss the results in section 9.4 
under the point of view that power-law inflation is an attractor solution of 
the corresponding field equations. 

9.1 Differential— geometrical properties 

Eq. (9.8) defines a self-similar spacetime: if we multiply the metric ds 2 by 
an arbitrary positive constant a 2 , then the resulting ds 2 = a 2 ds 2 is isometric 
to ds 2 . Proof: We perform a coordinate transformation t — at, x — b(a,p)x 
. . . On the other hand, the de Sitter spacetime eq. (9.7) is not self-similar, 
because it has a constant non-vanishing curvature scalar. 



ds 2 = dt 2 - \t\ 2p (dx 2 + dy 2 + dz 2 ) , p + 



(9.8) 



9.2 Scale-invariant Geld equations 
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Power-law inflation is intrinsically time-oriented. Proof: The gradient of 
the curvature scalar defines a temporal orientation. On the other hand, the 
expanding (H > 0) and the contracting (if < 0) de Sitter spacetime can be 
transformed into each other by a coordinate transformation, because both of 
them can be transformed to the closed Friedmann universe with scale factor 
cosh(iJt), which is an even function of t. This property is connected with 
the fact that eq. (9.8) gives a global description, whereas eq. (9.7) gives only 
a proper subset of the full de Sitter spacetime. 

Forp^oo, eq. (9.8) tends to eq. (9.7). Such a statement has to be taken 
with care, even for the case with real functions. Even more carefully one has 
to deal with spacetimes. The most often used limit — the Geroch-limit [86] 
of spacetimes — has the property that a symmetry, here we take it to be 
self-similarity, of all the elements of the sequence must also be a symmetry 
of the limit. 

From this it follows that the Geroch limit of spacetimes eq. (9.8) with 
p — > oo cannot be unique, moreover, it is just the one-parameter set eq. (9.7) 
parametrized by arbitrary values H > 0. 

9.2 Scale— invariant field equations 

A gravitational field equation is called scale-invariant, if to each solution 
ds 2 and to each positive constant c 2 the resulting homothetically equivalent 
metric ds 2 = c 2 ds 2 is also a solution. One example of such field equations is 
that one following from eq. (9.2). Moreover, no Lagrangian L = L(R) gives 
rise to a scale-invariant field equation which is not yet covered by eq. (9.2) 
already. 
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Secondly, for 



L = R/16nG - -g^Q&j + ^exp(^) , 



(9.9) 



we assume 8tvG = 1 henceforth, the homothetic transformation has to be 
accompanied by a suitable translation of $ to ensure scale-invariance. For 
/i^O, the value of V can be normalized to 1, or —1. 

A third example is the following: for the Kaluza-Klein ansatz eq. (9.10) 
the iV-dimensional Einstein equation Rab = is scale-invariant. We take 
the ansatz as 



where i,j = 0, ... 3; a, (3 = 4, . . . N — 1; A, B = ... iV — 1, and restrict 
to the warped product of the (N — 4)-dimensional internal space dr 2 with 
4-dimensional spacetime ds 2 . 

Let us consider the limits m — > oo and m — > of eq. (9.2). One gets 
L = exp(_R/i? ) and L = \n(R/R ), respectively. Both of them give rise to 
a field equation which is not scale- invariant: a homothetic transformation 
changes also the reference value R . For the second case this means a change 
of the cosmological constant. So we have a similar result as before: the limits 
exist, but they are not unique. 

9.3 Cosmological Friedmann models 

We consider the closed and open model in subsection 9.3.1. and the spatially 
flat model in subsection 9.3.2. 



dS 2 = ds 2 (x i ) + W{x l )dT 2 (x a ) , 



(9.10) 
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9.3.1 The closed and open models 

Let us start with the closed model. We restrict ourselves to a region where one 
has expansion, so we may use the cosmic scale factor as timelike coordinate: 

ds 2 = a 2 (<Q 2 (a) do 2 -da 2 ) , Q(a) > , (9.11) 

where da 2 , defined by 

da 2 = dr 2 + sin 2 rdVl 2 , dtf = d9 2 + sin 2 6dip 2 (9.12) 

is the positively curved 3-space of constant curvature. It holds: if the 00- 
component of the field equation, here: eq. (9.15) below, is fulfilled, then 
all other components are fulfilled, too. Such a statement holds true for all 
Friedmann models and "almost all sensible field theories" . With ansatz eqs. 
(9.11), (9.12) we get via 

B° Q = 3(a~ 4 Q- 2 + a- 3 Q- 3 dQ/da) (9.13) 

and 

R = Q( a - 3 Q- 3 dQ/da - a' 2 ) (9.14) 

the result: the field equation following from the Lagrangian (9.2) is fulfilled 
for metric (9.11), (9.12), if and only if 

mRR° - R 2 /2 + 3m(m - l)Q- 2 a - 3 dR/da = (9.15) 

holds. For m — 1, which represents Einstein's theory, no solution exists. 
For all other values m, eqs. (9.13), (9.14) and (9.15) lead to a second order 
equation for Q(a). 

We look now for solutions with vanishing Rq , i.e., with eq. (9.13) we 

get 

Q = C/a, C = const. >0. (9.16) 
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By the way, metric (9.11), (9.12) is self-similar if and only if eq. (9.16) holds. 
From eqs. (9.14), (9.16) we get 

R = D/a 2 , D = —6(1 + 1/C 2 ) . (9.17) 

We insert eqs. (9.16), (9.17) into eq. (9.15) and get 

C = (2m 2 - 2m - l) 1 ' 2 , (9.18) 

which fulfils eq. (9.16) if 

m> (l + V$)/2 or m<(l-Vs)/2 (9.19) 

holds. Inserting eq. (9.16) and eq. (9.18) into eq. (9.11) and introducing 
synchronized coordinates we get as a result: if eq. (9.19) holds, then 

ds 2 = dt 2 - t 2 da 2 /(2m 2 - 2m - 1) (9.20) 

is a solution of the fourth-order field equation following from Lagrangian 
(9.2). It is a self-similar solution, and no other self-similar solution describing 
a closed Friedmann model exists. 

For the open Friedmann model all things are analogous, one gets for 

(l-v / 3)/2<m< (l + VS)/2 (9.21) 

and with sinh r instead of sin r in eq. (9.12) the only self-similar open 
solution, which is flat for m G {0, 1}, 



ds 2 = dt 2 - t 2 da 2 /(-2m 2 + 2m + 1) . 



(9.22) 
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9.3.2 The spatially flat model 

The field equation for the spatially flat model can be deduced from that 
one of a closed model by a limiting procedure as follows: we insert the 
transformation r — > er, a — > a/e, Q — > Q^ 2 and apply the limit £ — > 
afterwards. One gets via 

lim sin(er) = r (9.23) 

the metric 

ds 2 = a 2 (Q 2 (a)da 2 - dx 2 - dy 2 - dz 2 ) (9.24) 
with unchanged eqs. (9.13) and (9.15), whereas eq. (9.14) yields 

R = 6a- 3 Q- 3 dQ/da. (9.25) 

The trivial solutions are the flat Minkowski spacetime and the model with 
constant value of Q, i.e., R — 0, a valid solution for m > 1 only, which has 
the geometry of Friedmann's radiation model. 

Now, we consider only regions with non-vanishing curvature scalar. For 
the next step we apply the fact that the spatially flat model has one symme- 
try more than the closed one: the spatial part of the metric is self-similar. 
In the coordinates eq. (9.24) this means that each solution Q(a) may be 
multiplied by an arbitrary constant. To cancel this arbitrariness, we define 
a new function 

P(a) = d(lnQ)/da. (9.26) 

We insert eqs. (9.13), (9.25) into eq. (9.15), and then eq. (9.26) into the 
resulting second order equation for Q. We get the first order equation for P 

= m (m - l)dP/da + (m - 1)(1 - 2m)P 2 + m(4 - 3m)P/a . (9.27) 
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As it must be the case, for m — 0, 1, only P = is a solution. For m = 1/2, 

P~a 5 and Q = exp(ca 6 ) , (9.28) 

where c denotes an integration constant. For the other values m we define 
z = aP as new dependent and t = In a as new independent variable. Eq. 
(9.27) then becomes 

= dz/dt + gz 2 -bz, g = l/m-2, b = (4m - 5)/(m - 1) . (9.29) 
For m = 5/4 we get 

2 = -5/(6t-c), i.e., P = -5/(6aln(a/c)), Q = (ln(a/c))^ 5/6 . (9.30) 

For the other values we get 

z = f/(eM~bt + c) - 1) , / = , P = //(eV- 6 - a) 

Q = {±a b + c) 1/g . (9.31) 

For m — > 1/2 we get 6 — > 6 and 1/p — > oo; for m — > 5/4 we get 6^0 and 
l/(? — > —5/6, so the two special cases could also have been obtained by a 
limiting procedure from eq. (9.31). 

Metric (9.24) with eq. (9.31) can be explicitly written in synchronized 
coordinates for special examples only, see e.g. Burd and Barrow [45]. 

9.4 Discussion of power-law inflation 

We have considered scale-invariant field equations. The three examples men- 
tioned in section 9.2 can be transformed into each other by a conformal trans- 
formation of the four-dimensional spacetime metric. The parameters of eqs. 
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(9.2) and (9.9) are related by y/3 p = y/2(2 - m)/(m - 1), cf. [214], a similar 
relation exists to the internal dimension in eq. (9.10), one has 

m = l + l/{l + V3[l + 2/(N - 4)] ±1/2 } , (9.32) 

for details see [31]. The necessary conformal factor is a suitable power of the 
curvature scalar eq. (9.2). A further conformal transformation in addition 
with the field redefinition 9 = tanh $ leads to the conformally coupled scalar 
field 9 in the potential 

(1 + 9) 2+ » (1 - 9) 2 -» , (9.33) 

the conformal factor being cosh 4 $, cf. [213]. So, equations stemming from 
quite different physical foundations are seen to be equivalent. We have looked 
at them from the point of view of self-similar solutions and of limiting pro- 
cesses changing the type of symmetry. 

The general solution for the spatially flat Friedmann models in fourth- 
order gravity eq. (9.2), eqs. (9.24), (9.31) can be written for small c in 
synchronized coordinates as follows 

ds 2 = dt 2 -a 2 (t)(dx 2 + dy 2 + dz 2 ), a(t)=t p (l+er bp + 0(e 2 )). (9.34) 

e = gives exact power-law inflation with p = (m — l)(2m — l)/(2 — m) 
and b = (4m — 5)/(m — 1). We see: in the range 1 < m < 2, power-law 
inflation is an attractor solution within the set of spatially flat Friedmann 
models for m > 5/4 only, because otherwise b would become negative. The 
generalization to polar inflation with m > 2 and p < is similar. 

9.5 The cosmic no hair theorem 

Following [129], we now discuss the cosmic no hair theorem, which tells under 
which circumstances the de Sitter spacetime represents an attractor solution 



122 9 MODELS WITH DE SITTER AND POWER-LAW INFLATION 



within the set of other nearby solutions. This property ensures the inflation- 
ary model to be a typical solution. The notion "cosmic no hair theorem" is 
chosen because of its analogous properties to the "no hair theorem" for black 
holes, for which one can cf. [194]. 

After a general introduction we restrict our consideration to spatially flat 
Friedmann models. In this section, we choose gravitational Lagrangians 



and linear combinations of them. They are motivated from trials how to 
overcome the non-renormalizability [248] of Einstein's theory of gravity. Re- 
sults are: For arbitrary k, i.e., for arbitrarily large order 2 k + 4 of the field 
equation, one can always find examples where the attractor property takes 
place. Such examples necessarily need a non-vanishing i? 2 -term. The main 
formulas do not depend on the dimension, so one gets similar results also for 
1+1-dimensional gravity and for Kaluza-Klein cosmology. 

9.5.1 Introduction to no hair theorems 

Over the years, the notion "no hair conjecture" drifted to "no hair theo- 
rem" without possessing a generally accepted formulation or even a complete 
proof. Several trials have been made to formulate and prove it at least for 
certain special cases. They all have the overall structure: "For a geometri- 
cally defined class of spacetimes and physically motivated properties of the 
energy-momentum tensor, all the solutions of the gravitational field equation 
asymptotically converge to a space of constant curvature." 

The paper Weyl [269] is cited in [13] with the phrase: "The behaviour 
of every world satisfying certain natural homogeneity conditions in the large 




(9.35) 
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follows the de Sitter solution asymptotically." to be the first published ver- 
sion of the no hair conjecture. Barrow and Gotz [13] apply the formulation: 
"All ever-expanding universes with A > approach the de Sitter spacetime 
locally." The first proof of the stability of the de Sitter solution, in this case 
within the steady-state theory, is due to Hoyle and Narlikar 1963 [112]. 

The probability of inflation is large if the no hair theorem is valid, cf. Alt- 
shuler [1]. Peter, Polarski, Starobinsky [181] and Kofman, Linde, Starobinsky 
[130] compared the double-inflationary models with cosmological observa- 
tions. Brauer, Rendall, Reula [40] discussed the no hair conjecture within 
Newtonian cosmological models. Hiibner and Ehlers [113] and Burd [44] 
considered inflation in an open Friedmann universe and have noted that in- 
flationary models need not to be spatially flat. 

Gibbons and Hawking 1977 [88] have found two of the earliest strict 
results on the no hair conjecture for Einstein's theory, cf. also Hawking, 
Moss [109] and Demianski [62]. Barrow 1986 [12] gave examples that the no 
hair conjecture fails if the energy condition is relaxed and points out, that 
this is necessary to solve the graceful exit problem. He uses the formulation of 
the no hair conjecture "in the presence of an effective cosmological constant, 
stemming e.g. from viscosity, the de Sitter spacetime is a stable asymptotic 
solution" . 

Usually, energy inequalities are presumed for formulating the no hair 
conjecture. Nakao et al. [168] found some cases where it remains valid also 
for negative Abbott-Deser mass. The latter goes over to the well-known 
ADM-mass [7] for A — > 0. In Murphy [165], viscosity terms as source are 
considered to get a singularity-free cosmological model. In the eighties, these 
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non-singular models with viscosity where reinterpreted as inflationary ones, 
cf. Oleak [173]. 

In the three papers [186], Prigogine et al. developed a phenomenological 
model of particle and entropy creation. It allows particle creation from space- 
time curvature, but the inverse procedure, i.e. particle decay into spacetime 
curvature is forbidden. This breaks the t — > — t-invariance of the model. 
Within that model, the expanding de Sitter spacetime is an attractor solu- 
tion independently of the initial fluctuations; this means, only the expanding 
de Sitter solution is thermodynamically possible. 

Vilenkin 1992 [262] and Borde [33], [34] discussed future-eternal inflating 
universe models; they must have a singularity if the condition D: "There is 
at least one point p such that for some point q to the future of p the volume 
of the difference of the pasts of p and q is finite" is fulfilled. Mondaini and 
Vilar [157] have considered recollapse and the no hair conjecture in closed 
higher-dimensional Friedmann models. 

In Shiromizu et al. [237], Shibata et al. [235] and Chiba, Maeda [52] 
the following argument is given: If the matter distribution is too clumpy, 
then a large number of small black holes appears. Then one should look for 
an inflationary scenario where these black holes are harmless. They cannot 
clump together to one giant black hole because of the exponential expansion 
of the universe; this explains the existing upper bound of the mass of black 
holes in the quasi-de Sitter model: above 

M crit = ^= (9.36) 

there do not exist horizons; this restriction is called cosmic hoop conjecture. 
The paper [52] from 1994 is entitled "Cosmic hoop conjecture?", and that 
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conjecture is formulated as: when an apparent horizon forms in a gravita- 
tional collapse, the matter must be sufficiently compactified such that the 
circumference C satisfies the condition 

C < AnM < 47rM crit , (9.37) 

where M is the Abbott-Deser mass of the collapsed body. 

Xu, Li and Liu [272] proved the instability of the anti-de Sitter spacetime 
against gravitational waves and dust matter perturbations in 1994. Coley and 
Tavakol [54] discussed the robustness of the cosmic no hair conjecture under 
using the concept of the structural stability. Sirousse-Zia [239] considered 
the Bianchi type IX model in Einstein's theory with a positive A-term and 
got an asymptotic isotropization of the mixmaster model. She uses methods 
of Belinsky et al. [22]. Miiller [162] used L = R 2 and discussed the power- 
asymptotes of Bianchi models. Barrow and Sirousse-Zia [16] discussed the 
mixmaster i? 2 -model and the question, under which conditions the Bianchi 
type IX model becomes asymptotic de Sitter. In Yokoyama, Maeda [273] 
and Cotsakis et al. [57] the no hair conjecture for Bianchi type IX models 
and Einstein's theory with a cosmological term are discussed; they get R 2 — 
inflation in anisotropic universe models and typically, an initial anisotropy 
helps to enhance inflation. They got some recollapsing solutions besides those 
converging to the de Sitter solution. In Spindel [241] also general Bianchi 
type I models in general dimensions are considered. 

Breizman et al. 1970 [41] considered 

L = R + ai? 4/3 (9.38) 

to get a singularity-free model, the solutions are of a quasi de Sitter type. 
The papers Barrow [12] from 1987, Pollock [185], Miiller et al. [164] and 
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Mijic et al. [155] consider the no hair conjecture for R 2 models, they use 
the formulation "asymptotical de Sitter, at least on patch". The restriction 
"on patch" is not strictly defined but refers to a kind of local validity of the 
statement, e.g., in a region being covered by one single synchronized system 
of reference in which the spatial curvature is non-positive and the energy 
conditions are fulfilled. The Starobinsky model is one of those which does 
not need an additional inflaton field to get the desired quasi de Sitter stage. 

One should observe a notational change: Originally, 

L = R + aR 2 In R (9.39) 

was called Starobinsky model, whereas 

L = R + aR 2 (9.40) 

got the name "improved Starobinsky model" - but now the latter carries sim- 
ply the name "Starobinsky model", see e.g. Starobinsky 1983 [244], Hwang 
[114], Gottlober et al. [95] and Amendola et al. [3] and [4]. For the infla- 
tionary phase, both versions are quite similar. A further result of the papers 
[12] and [185] is that by the addition of a cosmological term, the Starobinsky 
model leads naturally to double inflation. Let us comment this result: It is 
correct, but one should add that this is got at the price of getting a "graceful 
exit problem" , by this phrase there is meant the problem of how to finish the 
inflationary phase dynamically - in the Starobinsky model this problem is 
automatically solved by the fact that the quasi de Sitter phase is a transient 
attractor only. 

The papers Maeda [150], Maeda et al. [152] and Barrow, Saich [15] discuss 
the no hair conjecture within i? 2 -models and found inflation as a transient 
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attractor in fourth-order gravity. The papers Feldman [75], Rogers, Isaacson 
[192] and [220] discuss the stability of inflation in i? 2 -gravity. The papers 
Bicak, Podolsky [27], Cotsakis, Flessas [58], Borzeszkowski, Treder [36] and 
Contreras et al. [56] discuss generalized cosmic no hair theorems for quasi 
exponential expansion. 

In Starobinsky 1983 [244] the no hair theorem for Einstein's theory with 
a positive A-term is tackled by using a power series expansion as ansatz to 
describe a general spacetime. This includes a definition of an asymptotic de 
Sitter spacetime even for inhomogeneous models. However, the convergence 
of the sequence is not rigorously proven. In [245], the ansatz [244] was gen- 
eralized to consider also the cosmic no hair theorem for L = R 2 . Shiromizu 
et al. [237] discussed an inflationary inhomogeneous scenario and mentioned 
the open problem how to define asymptotic de Sitter spacetimes. 

Jensen, Stein-Schabes [122] consider the no hair theorem for a special 
class of inhomogeneous models and give partial proofs. Morris [161] considers 
inhomogeneous models for R + i? 2 -cosmology. In Calzetta, Sakellariadou [47] 
inflation in inhomogeneous but spherically symmetric cosmological models 
is obtained only if the Cauchy data are homogeneous over several horizon 
lengths. The analogous problem is considered in Nakao et al. [166], [167]; 
in Shinkai, Maeda [236] also the inclusion of colliding plane gravitational 
waves is considered, they give a numerical support of the no hair conjecture 
by concentrating on the dynamics of gravitational waves. Berkin [24] gets 
as further result, that for L = f(R), a diagonal Bianchi metric is always 
possible. Barrow and Sirousse-Zia [16] and Spindel [241] also worked on 
the diagonalization problem. They apply the diagonalization condition of 
MacCallum et al. [146], [147]. 
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Kottler 1918 [131] found the spherically symmetric vacuum solution for 
Einstein's theory with A-term 

dr 2 

ds 2 = A{r)dt 2 - — - r 2 {d9 2 + sin 2 6 d<f) 2 ) (9.41) 

with 

A(r) = l- — -K 2 . (9.42) 
r 3 

At the horizon A(r) = the Killing vector changes its sign and one gets 
by interchanging the coordinates t and r the corresponding Kantowski-Sachs 
model. The cosmic no hair conjecture within Kantowski-Sachs models and 
A > is discussed by Moniz [159]. To get a theory renormalizable to one- 
loop order one needs at least curvature squared terms in the Lagrangian, 
and they lead to a fourth-order theory, Stelle [248]. But there exist several 
further motivations to consider fourth-order theories; to get an overview for 
early papers on fourth-order gravity see e.g. [199]. 

Let us now discuss sixth-order equations. Stelle [248] 1977 considers 
mainly fourth-order i? 2 -models; in the introduction he mentioned that in the 
next order, terms like 

R 3 + /,',,:/, /,"' : ' (9.43) 

become admissible, but the pure i? 3 -term is not admissible. The paper Buch- 
dahl 1951 [42] deals with Lagrangians of arbitrarily high order. Its results 
are applied in [218], Quandt et al. [187], and Battaglia Mayer et al. [18] to 
general Lagrangians F(R, □). In the paper Vilkovisky [263], the Sacharov- 
approach was generalized. The original idea of A. Sakharov [195] from 1967, 
see the commented reprint from 2000, was to define higher order curvature 
corrections to the Einstein action to get a kind of elasticity of the vacuum. 
Then the usual breakdown of measurements at the Planck length, such a 
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short de Broglie wave length corresponds to such a large mass which makes 
the measuring apparatus to a black hole, is softened. 

Vilkovisky [263] discusses the effective gravitational action in the form 
Rf(n)R, where 



Martin and Mazzitelli [153] discuss the non-local Lagrangian R^R as con- 
formal anomaly in two dimensions. Treder 1991 [257] used higher-order La- 
grangians, especially i? 2 -terms, and he mentioned that for 

R + R,R, k g lk (9.45) 

a sixth-order field equation appears. Remark: This Lagrangian leads to the 
same field equation as R—RnR. Lu and Wise [144] consider the gravitational 
Lagrangian as a sequence S — S + Si + S 2 + ■ ■ . ordered with respect to 
physical dimension. So, S = R and Si sums up the _R 2 -terms. Kirsten et 
al. [127] consider the effective Lagrangian for self-interacting scalar fields; in 
the renormalized action, the term 

Un (9.46) 



c + R 



appears. Wands [266] classifies Lagrangians of the type F(R,<&)nR and 
mentions that not all of them can be conformally transformed to Einstein's 
theory. Amendola [2] considers the Lagrangian $ 2 Di?, Gottlober et al. [97] 
the Lagrangian i?Di?, Gottlober et al. [94] and Schmidt [219] discuss double 
inflation from $ and i? 2 -terms. Besides Berkin [24] and Berkin, Maeda 

[25] consider the de Sitter spacetime as attractor solution for field equations 
where the variational derivative of the term CijkiC tjH is included. 
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This section is organized as follows: subsection 9.5.2 compares several 
possible definitions of an asymptotic de Sitter spacetime. Subsection 9.5.3 
deals with the Lagrangian and corresponding field equations for higher-order 
gravity, and 9.5.4 gives the no hair theorem for it. In subsection 9.5.5, we 
determine under which circumstances the Bianchi models in higher-order 
gravity can be written in diagonal form without loss of generality, and in 
subsection 9.5.6 we discuss and summarize the results. 

9.5.2 Definition of an asymptotic de Sitter spacetime 

In this subsection we want to compare some possible definitions of an asymp- 
totic de Sitter spacetime. To this end let us consider the metric 

n 

ds 2 = dt 2 -e 2a ^J2 d ( xt ) 2 ( 9 - 47 ) 
i=i 

which is the spatially flat Friedmann model in n spatial dimensions. We 
consider all values n > 1, but then concentrate on the usual case n = 3. The 
Hubble parameter is 



H = a = ^. (9.48) 
at 



We get 



Roo = -n (^j- + H 2 ^j , R = -2nH - n{n + l)H 2 . (9.49) 

Then it holds: The following 4 conditions for metric (9.47) are equivalent. 
1: The spacetime is flat. 2: It holds R = R 00 = 0. 3: The curvature 
invariant RijR 1 ^ vanishes. 4: Either a = const, or 

n = l and a = In \t — t \ + const. (9.50) 

For the last case with n — 1 one has to observe that 

ds 2 = dt 2 - (t - t ) 2 dx 2 (9.51) 



9.5 The cosmic no hair theorem 



131 



represents flat spacetime in polar coordinates. For the proof we use the 
identity 

RijW = (R 00 ) 2 + \{R~ Roo) 2 ■ (9.52) 

An analogous characterization is valid for the de Sitter spacetime. The 
following 4 conditions for metric (9.47) are equivalent. 1: It is a non-flat 
spacetime of constant curvature. 2: R 00 = R/(n+ 1) = const. 7^ 0. 3: 
(n + l)RijR ij = R 2 = const. ^ 0. 4: Either H = const. ^ or 

n = 1 and (ds 2 = dt 2 - sm 2 (\t)dx 2 or ds 2 = dt 2 - sinh 2 (At)cfe 2 ) . (9.53) 

For n = 1, the de Sitter spacetime and anti-de Sitter spacetime differ by the 
factor —1 in front of the metric only. For n > 1, only the de Sitter spacetime, 
having R < 0, is covered, because the anti-de Sitter spacetime cannot be 
represented as spatially flat Friedmann model. Our result shows that within 
the class of spatially flat Friedmann models, a characterization of the de 
Sitter spacetime using polynomial curvature invariants only, is possible. 

Next, let us look for isometries leaving the form of the metric (9.47) 
invariant. Besides spatial isometries, the map a — > a defined by 

a(t) = c + a(±t + 1 ), c, t = const. (9.54) 

leads to an isometric spacetime. The simplest expressions being invariant by 
such a transformation are H 2 and H. We take a as dimensionless, then H is 
an inverse time and H an inverse time squared. Let H 7^ in the following. 
The expression 

e := EE- 2 (9.55) 

is the simplest dimensionless quantity defined for the spatially flat Friedmann 
models and being invariant with respect to this transformation. Let n > 1 in 
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the following, then it holds: Two metrics of type (9.47) are isometric if and 
only if the corresponding functions a and a are related by this transforma- 
tion. It follows: Metric (9.47) with H ^ represents the de Sitter spacetime 
if and only if e = 0. 

All dimensionless invariants containing at most second order derivatives 
of the metric can be expressed as f(e), where / is any given function. But if 
one has no restriction to the order, one gets a sequence of further invariants 

e 2 = HH~ 3 , . . . , e p = -^H'^ 1 . (9.56) 

Let H > in metric (9.47) with n > 1. We call it an asymptotic de Sitter 
spacetime if 

lim = const. ^ (9.57) 

or 

lim R 2 = const. > and lim (n + l)RijRij - R 2 = (9.58) 

t^oo t^oo 

or for some natural number p one has 

hm Ej = for l<j<p. (9.59) 

In general, all these definitions are different. Using the identity 

RijW =n 2 {H + H 2 ) 2 + n{H + nH 2 ) 2 (9.60) 

we will see that all these definitions lead to the same result if we restrict our- 
selves to the set of solutions of the higher-order field equations. The problem 
is that none of the above definitions can be generalized to inhomogeneous 
models. One should find a polynomial curvature invariant which equals a 
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positive constant if and only if the spacetime is locally the de Sitter space- 
time. To our knowledge, such an invariant cannot be found in the literature, 
but also the non-existence of such an invariant has not been proven up to 
now. For the anti-de Sitter spacetime, a partial result is given in section 4.3. 

This situation is quite different for the positive definite case: For signature 
(+ + ++) and Sij = Rij - it holds: 

I=(R- R ) 2 + C ljkl C^ kl + = (9.61) 

if and only if the is a space of constant curvature R . So I — > is a 
suitable definition of an asymptotic space of constant curvature. One possi- 
bility exists, however, for the Lorentz signature case, if one allows additional 
structure as follows: An ideal fluid has an energy-momentum tensor 

T ij = (p + p)uiUj - P9ij (9.62) 

where Ui is a continuous vector field with UiU % = 1. For matter with equation 
of state p = —p, the equation T lJ .j = implies p = const., and so every 
solution of Einstein's theory with such matter is isometric to a vacuum so- 
lution of Einstein's theory with a cosmo logical term. The inverse statement, 
however, is valid only locally: 

Given a vacuum solution of Einstein's theory with a A-term, one has 
to find continuous timelike unit vector fields which need not to exist from 
topological reasons. And if they exist, they are not at all unique. So, it 
becomes possible to define an invariant J which vanishes if and only if the 
spacetime is de Sitter by transvecting the curvature tensor with u l v? and/or 
g % i and suitable linear and quadratic combinations of such terms. Then time 
t becomes defined by the streamlines of the vector u l . If one defines the 
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asymptotic de Sitter spacetime by J — > as t — > oo, then it turns out, 
that this definition is not independent of the vector field u\ 

9.5.3 Lagrangian F(R, OR, B 2 R, . . . , D k R) 

Let us consider the Lagrangian density L given by 

L = F(R, OR, a 2 R,...,n k R)^=g~ (9.63) 

where R is the curvature scalar, □ the D'Alembertian and gij the metric 
of a Pseudo-Riemannian Vd of dimension D > 2 and arbitrary signature; 
g — — | det \ . The main application will be D = 4 and metric signature 
(H ). F is supposed to be a sufficiently smooth function of its argu- 
ments, preferably a polynomial. Buchdahl [42] already dealt with such kind 
of Lagrangians in 1951, but then it became quiet of them for decades. 

The variational derivative of L with respect to the metric yields the tensor 

P i] = ^= (9-64) 

The components of this tensor read 

Pij = GRij - - G, tJ + 9ij UG + X tJ (9.65) 

where the semicolon denotes the covariant derivative, Rij the Ricci tensor, 
and ^ 

*a = E IgiAM^RrU - F^p^R].^ (9.66) 

A=l Z 

having the round symmetrization brackets in its last term. For k — 0, i.e. 
F = F(R), a case considered in subsection 9.5.4, the tensor identically 
vanishes. It remains to define the expressions Fa, A — 0, ... ,k . 
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and for A — k — 1 , . . . , 

dF 

Fa = UF A+1 + (9.68) 

and finally G = F . The brackets are essential, for any scalar $ it holds 

□ (<&,«) - (□$), = R* $ ; . (9.69) 

Inserting $ = O m R into this equation, one gets identities to be applied in 
the sequel without further notice. The covariant form of energy- momentum 
conservation reads 

P) ;i = (9.70) 

and Pij identically vanishes if and only if F is a divergence, i.e., locally there 
can be found a vector v l such that F = v 1 ., holds. Remark: Even for com- 
pact manifolds without boundary the restriction "locally" is unavoidable, for 
example, let D = 2 and V 2 be the Riemannian two-sphere S 2 with arbitrary 
positive definite metric. R is a divergence, but there do not exist continuous 
vector fields v l fulfilling R = v 1 ^ on the whole S 2 . Example: for m, n > 
it holds 

U m R O n R - R U m+n R = divergence . (9.71) 

So, the terms O m R U n R with naturals m and n can be restricted to the case 
m = without loss of generality. 

9.5.4 No hair theorems for higher-order gravity 

For n > 1, the n + 1-dimensional de Sitter spacetime is an attractor solution 
for the field equation derived from the Lagrangian 



R (n+l)/2 _ (£L72) 
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It is not an attractor solution for the Lagrangian RD k R and k > 0. There 
exist combinations of coefficients q, such that the de Sitter spacetime is 
an attractor solution for the field equation derived from the more general 
Lagrangian 

m 

L = c i? (n+1)/2 + c k RO k R. (9.73) 

fc=0 

Idea of Proof: Concerning fourth-order gravity this method was previously 
applied e.g. by Barrow 1986 [12]. The de Sitter spacetime is an exact solution 
for the field equation, if and only if 2RG — (n + l)F. If we make the ansatz 

a(t) = l + P(t) (9.74) 

n+1 

we get the linearized field equation = (3 + n/3 for the Lagrangian R~ . 
For the Lagrangian L = RO k R we get the linearized field equation O h R = 
(□ fe i?) . For the characteristic polynomial we get a recursive formula such 
that the next order is received from the previous one by multiplying with 
• x • (x + n). We get the roots x\ = —n — 1, x 2 = —n (fc-fold), x 3 = 
(fc-fold) and x 4 = +1. Because of the last root the de Sitter spacetime is not 
an attractor solution. For the Lagrangian (9.73) we get the characteristic 
polynomial 



P(x) = x(x + n) 



(9.75) 



c + CkX^ix + n) k ~ l (x - l)(x + n + 1) 

k=l 

for the linearized field equation. The transformation z = x 2 + nx + ^ 
gives a polynomial Q(z) which can be solved explicitly. So one can find 
those combinations of the coefficients such that the de Sitter spacetime is an 
attractor solution. 



It turned out that all the variants of the definition of an asymptotic de 
Sitter solution given in subsection 9.5.2 lead to the same class of solutions, 
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i.e., the validity of the theorems written below does not depend on which 
of the variants of definition of an asymptotic de Sitter spacetime listed is 
applied. For the 6th-order case we can summarize as follows: Let 

L = R 2 + ci R aR (9.76) 

and 

L E = R - l -L (9.77) 
with length I > 0. Then the following statements are equivalent. 

1. The Newtonian limit of L E is well-behaved, and the potential consists 
of terms e~ ar /r with a > only. 

2. The de Sitter spacetime with H = l/l is an attractor solution for L 
in the set of spatially flat Friedmann models, and this can already be seen 
from the linearized field equation. 

3. The coefficient c\ > 0, and the graceful exit problem is solved for the 
quasi de Sitter phase H < l/l of L E . 

4. I 2 = /q + 1\ such that I 2 c\ = /q l\ has a solution with < l < h- 

5. < ci < / 2 /4. 

A formulation which includes also the marginally well-behaved cases reads 
as follows: Let L and L E as in the previous result, then the following state- 
ments are equivalent: 

I. The Newtonian limit of L E is well-behaved, for the potential we 
allow 1/r and terms like 

^l e ~ ar with a>0 (9.78) 
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and a polynomial P. 

2. The de Sitter spacetime with H = l/l cannot be ruled out to be an 
attractor solution for L in the set of spatially flat Friedmann models if one 
considers the linearized field equation only. 

3. Le is tachyonic-free. 

4. I 2 = /q + 1\ with I 2 d = /g l\ has a solution with < l < h- 

5. For the coefficients we have < c\ < j. 

Of course, it would be interesting what happens in the region where the 
linearized equation does not suffice to decide. 

9.5.5 Diagonalizability of Bianchi models 

A Bianchi model can always be written as 

ds 2 = dt 2 - g a/3 (t)a a a^ (9.79) 

where g a/ 3 is positive definite and cr a are the characterizing one-forms. It 
holds 

da^ = - l -Cl^N^ (9.80) 

with structure constants CT a p of the corresponding Bianchi type. It belongs 
to class A if = 0. The Abelian group, Bianchi type I, and the rotation 
group, Bianchi type IX, both belong to class A. In most cases, the g a p are 
written in diagonal form; it is a non-trivial problem to decide under which 
circumstances this can be done without loss of generality. For Einstein's 
theory, this problem is solved. One of its results reads: If a Bianchi model 
of class A, except Bianchi types I and II, has a diagonal energy-momentum 
tensor, then the metric g a /3(t) can be chosen in diagonal form. Here, the 
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energy-momentum tensor is called diagonal, if it is diagonal in the basis 
(dt, a\ a 2 , a 3 ). 

This result rests of course on Einstein's theory and cannot be directly 
applied to higher-order gravity For fourth-order gravity following from a 
Lagrangian L = f(R) considered in an interval of R- values where 

m ■ m * o 

one can do the following: The application of the conformal equivalence the- 
orem, cf. e.g. [29] or [218], is possible, the conformal factor depends on t 
only, so the diagonal form of metric does not change. The conformal pic- 
ture gives Einstein's theory with a minimally coupled scalar field as source; 
the energy-momentum tensor is automatically diagonal. So, in this class of 
fourth-order theories of gravity, we can apply the above cited theorem of 
MacCallum et al. [146], [147], where the initial- value formulation of General 
Relativity is applied. As example we formulate: All solutions of Bianchi 
type IX of fourth-order gravity following from L = R 2 considered in a region 
where R ^ can be written in diagonal form. Consequently, the ansatz used 
by Barrow and Sirousse-Zia [16] for this problem is already the most general 
one, cf. Spindel [241]. 

For fourth-order gravity of a more complicated structure, however, things 
are more involved; example: Let 

L = R + aR 2 + bC ijkl C ijkl (9.82) 

with ah 7^ 0. Then there exist Bianchi type IX models which cannot be 
written in diagonal form. This is a non-trivial statement because its proof 
needs a careful analysis of the correspondingly allowed initial values in the 
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Cauchy problem. To understand the difference between the cases 5 = and 
b 7^ it proves useful to perform the analysis independently of the above cited 
papers [46]. For simplicity, we restrict to Bianchi type I. Then the internal 
metric of the hypersurface [t = 0] is flat and we can choose as initial value 
daffi®) — $a/3- Spatial rotations do not change this equation, and we can take 
advantage of them to diagonalize the second fundamental form -|<7 ai g(0). 

First case: 6 = 0. As additional initial conditions one has only R(0) and 
■^R(0). The field equation ensures g a p{t) to remain diagonal for all times. 

Second case: 6^0. Then one has further initial data 

^fctf(O) • ( 9 - 83 ) 

In the generic case, they cannot be brought to diagonal form simultaneously 
with 4^ aj g(0). This excludes a diagonal form of the whole solution. To 
complete the proof, one has of course to check that these initial data are 
not in contradiction to the constraint equations. This case has the following 
relation to the Bicknell theorem: Just for this case & ^ 0, the conformal 
relation to Einstein's theory breaks down, and if one tries to re-interpret the 
variational derivative of C^mC 1 ^ 1 as energy-momentum tensor then it turns 
out to be non-diagonal generically, and the theorem cannot be applied. 

For higher-order gravity, the situation becomes even more involved. For 
a special class of theories, however, the diagonalizability condition is exactly 
the same as in Einstein's theory: If 

m 

L = R + Y, a k Rn k R, a m ^0, (9.84) 

then in a region where 2L ^ R the Cauchy data are the data of General 
Relativity, R(0), and the first 2m + 1 temporal derivatives of R at t — 0. All 
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terms with the higher derivatives behave as an energy-momentum tensor in 
diagonal form, and so the classical theorem applies. Let us comment on the 
restriction 2L ^ R supposed above: F = G = represents a singular point 
of the differential equation; and for the Lagrangian given here G = 2L/R — 1. 
For fourth-order gravity defined by a non-linear Lagrangian L(R) one has 

dL 

and G = defines the critical value of the curvature scalar. 



9.5.6 Discussion of no hair theorems 

In subsection 9.5.4 we have shown: The results of the Starobinsky model 
are structurally stable with respect to the addition of a sixth-order term 
~ RUR ) if the coefficients fufil certain inequalities. For the eighth-order case 
we got: For 

L = R 2 + Cl RDR + c 2 RanR, c 2 + (9.85) 

and the case n > 1 the de Sitter spacetime with H = 1 is an attractor solution 
in the set of spatially flat Friedmann models if and only if the following 
inequalities are fulfilled: 

< ci < — < c 2 < * 2 (9.86) 
n + 1 (n + iy 

and 



ci > -{n 2 + n + l)c 2 + ^(n 4 + 4n 3 + 4n 2 )c| + 4c 2 (9.87) 

These inequalities define an open region in the C\ — c 2 — plane whose boundary 
contains the origin; and for the other boundary points the linearized equation 
does not suffice to decide the attractor property. 
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This situation shall be called "semi-attractor" for simplicity. In a gen- 
eral context this notion is used to describe a situation where all Lyapunov 
coefficients have non-positive real parts and at least one of them is purely 
imaginary. 

In contrary to the 6th-order case, here we do not have a one-to-one 
correspondence, but a non-void open intersection with that parameter set 
having the Newtonian limit for L E well-behaved. 

To find out, whether another de Sitter spacetime with an arbitrary Hubble 
parameter if > is an attractor solution for the eighth-order field equation 
following from the above Lagrangian, one should remember that H has the 
physical dimension of an inverted time, c\ is a time squared, c 2 is a time to 
power 4. So, we have to replace C\ by ciH 2 and c 2 by c 2 H 4 in the above 
dimensionless inequalities to get the correct conditions. Example: 

< Cl H 2 < \ . 

4 

Let us summarize: Here for a theory of gravity of order higher than four 
the Newtonian limit and the attractor property of the de Sitter spacetime are 
systematically compared. It should be noted that the details of the theory 
sensibly depend on the numerical values of the corresponding coefficients. So, 
no general overall result about this class of theories is ever to be expected. 

The facts, contradicting each other from the first glance, field equations of 
higher than second order follow from quantum gravity considerations Mirz- 
abekian et al. [156] on the one hand but such equations are known to be 
unstable in general and are therefore unphysical on the other hand, have now 



9.5 The cosmic no hair theorem 



143 



considered in more details: It depends on the special circumstances whether 
such theories are unstable. 

We have found out that for the class of theories considered here, one of 
the typical indicators of instability - cosmological runaway-solutions - need 
not to exist, even for an arbitrarily high order of the field equation. It is 
an additional satisfactory result that this takes place in the same range of 
parameters where the Newtonian limit is well behaved. 
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10 Solutions of the Bach-Einstein equation 

For field equations of 4th order, following from a Lagrangian "Ricci scalar 
plus Weyl scalar", it is shown using methods of non-standard analysis that 
in a neighbourhood of Minkowski spacetime there do not exist regular static 
spherically symmetric solutions. With that, besides the known local expan- 
sions about r = and r = oo respectively, a global statement on the existence 
of such solutions is given, see [205]. Then, this result will be discussed in 
connection with Einstein's particle programme. Finally, in section 10.6 it will 
be shown, that cosmological Bianchi type I solutions of the Bach equation 
exist which fail to be conformally related to an Einstein space. 

10.1 Introduction to the Bach equation 

General Relativity Theory starts from the Einstein-Hilbert Lagrangian 



with the Ricci scalar R which leads to the Einstein vacuum field equation 
R* — being of second order in the metrical tensor g ik . Its validity is proven 
with high accuracy in spacetime regions where the curvature is small only. 
Therefore, the additional presence of a term being quadratic in the curvature 
is not excluded by the standard weak field experiments. 

Lagrangians with squared curvature have already been discussed by Weyl 
in 1919 [268], by Bach in 1921 [8] and Einstein in 1921. They were guided by 
ideas about conformal invariance, and Einstein [66] proposed to look seriously 
to such alternatives. In Bach [8] , the Lagrangian Lw with the the Weyl scalar 



L EH = R/2k 



(10.1) 



C, 



L w = C/2k, with C=--C ljkl C^ kl 



(10.2) 
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was of special interest. Then it became quiet of them for a couple of decades 
because of the brilliant results of General Relativity but not at least because 
of mathematical difficulties. 

Later, new interest in such equations arose with arguments coming from 
quantum gravity, cf. e.g. [257], Borzeszkowski, Treder, Yourgrau [37], Stelle 
[249] , Fiedler, Schimming [78] . In the large set of possible quadratic modifi- 
cations a linear combination 

L = L EH + l 2 L w (10.3) 

enjoyed a special interest, cf. e.g. Treder 1977 [257]. The coupling constant 
I has to be a length for dimensional reasons and it has to be a small one to 
avoid conflicts with observations: One often takes it to be the Heisenberg 
length, which is the Compton wave length of the proton, = 1.3 • 1CT 18 cm 
or the Planck length = 1.6 • 10~ 33 cm. 

In this chapter we consider the Lagrangian L eq. (10.3) in connection 
with Einstein's particle programme, see Einstein, Pauli 1943 [69]: One asks 
for spherically symmetric singularity-free asymptotically flat solutions of the 
vacuum field equations which shall be interpreted as particles, but this cannot 
be fulfilled within General Relativity itself; it is still a hope, cf. [35], to realize 
it in such 4th order field equations. 

Two partial answers have already been given: Stelle [249] has shown the 
gravitational potential of the linearized equations to be 

$(r) = — m/r + exp(— r/l)c/r for r < I ; (10.4) 

he also gave an expansion series about r = 0, and Fiedler, Schimming [78] 
proved its convergence and smoothness in a certain neighbourhood of r = 0, 
but said nothing about the convergence radius of this expansion. 
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Now we want to join these two local expansions. To this end we write 
down the field equations in different but equivalent versions section 10.2, 
calculate some linearizations section 10.3 and prove the statement of the 
summary in section 10.4 which will be followed by a discussion of Einstein's 
particle programme in section 10.5. 



10.2 Notations and field equations 

We start from a spacetime metric with signature (H ), the Riemann and 

Ricci tensor being defined by R l ' - fci = r* 7 fc — . . ., Rik = R\j^ respectively. 
The Weyl tensor C l j kl is the traceless part of the Riemann tensor and the 
Weyl scalar is given by 2C = CijkiC 1 ^ 1 . Light velocity is taken to be 1 and 
G = k/8tt is Newton's constant; g = dety^. 

Then we consider the Lagrangian 

C = ^g~L = y/=jj(R + 1 2 C)/2k + £ mat , (10.5) 

where £ mat is the matter Lagrangian. For writing down the corresponding 
field equation it is convenient to introduce the Bach tensor Bij, cf. Bach [8] 
and Wiinsch [271], beforehand. 



^ ' v /— rt//" - ' ^ " / ""' • (J(U)) 



It holds B\ = 0, Bj.j = 0, B^ = Bji, and By is conformally invariant 
of weight -1. Variation of eq. (10.5) with respect to leads to the field 
equation 

Rij - \dijR + fBij = KTij . (10.7) 
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Now let us consider the vacuum case Ty = 0. The trace of eq. (10.7) then 
simply reads R = 0, and eq. (10.7) becomes equivalent to the simpler one 

Rij + PBij = 0. (10.8) 

Writing k = l~ 2 and □ • = g %i - one gets also the equivalent system 

R = ; kRij + OR.J = 2R iabj R ab + ^ gij R ab R ab . (10.9) 

The static spherically symmetric line element in Schwarzschild coordi- 
nates reads 

ds 2 = (1 + (3)e- 2X dt 2 - (1 + - r 2 dtt 2 , (10.10) 

where (5 and A depend on r only. The dot means differentiation with respect 
to r, and defining 

a = 2(3 - 2r$ + 2r\{l + (3) + 2r 2 {\ 2 - A)(l + (3) + r 2 {(3 - 3A/3) , (10.11) 
( = ra, )ji = ar" 3 , i] 2 = 3f3r~ 1 , i] 3 = (r~ 3 , (10.12) 

the field equations (10.9) are just eqs. (2.15.a-c) and (2.16) of the paper [78]. 
To avoid the products rfji we define a new independent variable x by r = le x . 
Then rrji = r( i = drji/dx and we obtain 

= ^ + 3^-^, (10.13) 
= (k + rr] 3 /Q)r]' 2 + ^ + rj 3 + r 2 krn + r%r] 3 /Q + rr] 2 r] 3 /2 + r 3 r] 2 /4 , (10.14) 

= (1 + rr} 2 /3)r} 3 - 2r} x + 2kr] 2 + 2r] 3 - r 2 kr]i - 

r 3 r] 2 /2 - r 3 77i77 3 /6 + 2rrj 2 rj 3 /3 . (10.15) 

Conversely, if the system eqs. (10.13), (10.14), (10.15) is solved, the metric 
can be obtained by 

(3 = rr] 2 /3, A = (2rr] 2 + 2r] 2 + r%)/(6 + 2rr] 2 ) , A(0) = 0. (10.16) 
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10.3 Linearizations 

Now we make some approximations to obtain the rough behaviour of the 
solutions. First, for r — > oo, Stelle [249] has shown the following: Lineariza- 
tion about Minkowski spacetime leads, using our notations eqs. (10.4) and 
(10.10), to A = and /3 = 2$, where powers of $ are neglected and the term 
exp(r//)/r in $ must he suppressed because of asymptotical flatness. 

That means, loosely speaking, the right hand side of eq. (10.9) will be 
neglected. Connected with this one may doubt the relevance of the term 
cr -i e -r// j n e q_ (10.4) for m/0 and r — > oo, because it is small compared 
with the neglected terms. Further, a finite $(0) requires c = mm eq. (10.4). 
Then $ is just the Bopp-Podolsky potential stemming from 4th order elec- 
trodynamics, and it were Pechlaner and Sexl [180] who proposed this form 
as representing a gravitational potential. Now we have 

$( r ) = m (exp(-r/Z) - 1) /r = -m/l + rm/(2/ 2 ) + . . . , (10.17) 

but also this finite potential gives rise to a singularity: For one of the curva- 
ture invariants we get RijR 1 ^ ~ m 2 /r 4 as r — * 0, that means, the linearization 
makes no sense for this region. From this one can already expect that in a 
neighbourhood of Minkowski spacetime no regular solutions exist. 

Second, for r — > 0, in [78] it is shown that there exists a one-parameter 
family of solutions being singularity-free and analytical in a neighbourhood 
of r = 0. The parameter will be called e and can be defined as follows: 
neglecting the terms with r in the system eqs. (10.13), (10.14), (10.15) 
one obtains a linear system with constant coefficients possessing just a one- 
parameter family of solutions being regular at r = 0; it reads 

771 = ek 2 r , r] 2 = —5ekr , r] 3 = Aek 2 r . (10.18) 
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The factors k n are chosen such that e becomes dimensionless. Now one 
can take eq. (10.18) as the first term of a power series rji = S n a| n V" , 
and inserting this ansatz into that system, one iteratively obtains the coeffi- 
cients a\ n \ For n even, af^ vanishes: The r 3 -terms are /c 3 (e/14 + 10£ 2 /21), 
-k 2 (e/2 - 10e 2 /3) and A: 3 (3e/7 + 20e 2 /7) respectively. Furthermore, the 
r 2n_1 -term is always a suitable power of k times a polynomial in e of the 
order < n. 

Up to the r 2 -terms the corresponding metric (10.10) reads 

ds 2 = (1 + hker 2 /3)dt 2 - (1 + 5ker 2 /3)dr 2 - r 2 d£l 2 . (10.19) 



Third, we look for a linearization which holds uniformly for < r < oo. 
A glance at eqs. (10.16) and (10.10) shows that r]i = gives Minkowski 
spacetime, and therefore we neglect terms containing products of rji in the 
system. In other words, instead of the linearization used before we addition- 
ally retain the term r 2 r\\. Then we proceed as follows: from eqs. (10.13) and 
(10.15) we obtain 

m = r][ + 3% and (10.20) 

= -y[ ~ 577i + (e 2x - 4)77! + Pe^mim + rj[/G) , . 

V2 2k + le x (r)' 1 + 5r ]l + 6r] 1 )/3 ' 1 ' ' 

Inserting eqs. (10.20), (10.21) into eq. (10.14) one obtains a third order 
equation for r\\ only whose linearization reads 

= + 577'/ + V [ (2 - e 2x ) - 771(8 + 4e 2x ) . (10.22) 

In Schmidt and Miiller [231] the same vacuum field equations were discussed 
for axially symmetric Bianchi type I models. They possess a four-parameter 
group of isometries, too, and the essential field equation is also a third order 
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equation for one function. The difference is that in the present case spherical 
symmetry implies an explicit coordinate dependence. 

The solution of eq. (10.22) which is bounded for x — > — oo reads 



A second solution is 771 = — 12mr -4 leading to the Schwarzschild solution, 
this solution solves both the full and the linearized equations in accordance 
with the fact that it makes zero both sides of eq. (10.9), and the third one can 
be obtained from them up to quadrature by usual methods. A comparison 
of eq. (10.23) with eq. (10.18) gives 7 = 15d~ 3 . 

Now we insert this 771 into eqs. (10.21), (10.16) and neglect again powers 
of 771 ; then the metric reads 



This linearization has, in contrary to eqs. (10.4) and (10.18), the following 
property: to each r > and A > there exists an 5 > such that for 
—5 < e < S the relative error of the linearized solution (10.24) does not 
exceed A uniformly on the interval < r < r$. 

10.4 Non-standard analysis 

In this section we will prove that in a neighbourhood of Minkowski spacetime 
there do not exist any solutions. "Neighbourhood of Minkowski spacetime" 
is in general a concept requiring additional explanations because of the large 
number of different topologies discussed in literature. To make the above 



77! = 7 • (3e~ 4x + e~ 2x ) sinh e £ - 3e~ 3:c cosh e 



(10.23) 



(3 = 5e Ir 1 sinh(r/ x ) — cosh(r/ x ) 
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statement sufficiently strong we apply a quite weak topology here: given a 
5 > then all spacetimes being diffeomorphic to Minkowski spacetime and 
fulfilling \RijR lj \ < 5k 2 form a neighbourhood about Minkowski spacetime. 
In eq. (10.19) we have at r = 0, independently of higher order terms in r, 

RijR i3 ' = 100&V/3 and R 00 = 5ke . (10.25) 

Therefore it holds: the one-parameter family of solutions being regular at 
r = is invariantly characterized by the real parameter e, and a necessary 
condition for it to be within a neighbourhood of Minkowski spacetime is that 
e lies in a neighbourhood of 0. 

Now we suppose e to be an infinitesimal number, i.e., a positive number 
which is smaller than any positive real number. The mathematical theory 
dealing with such infinitesimals is called non-standard analysis, cf. Robin- 
son [191]. The clue is that one can handle non-standard numbers like real 
numbers. Further we need the so-called Permanence principle, which is also 
called the Robinson lemma: Let A(e) be an internal statement holding for 
all infinitesimals e. Then there exists a positive standard real 5 > such 
that A(e) holds for all e with < e < 5. The presumption "internal" says, 
roughly speaking, that in the formulation of A(.) the words "standard" and 
"infinitesimal" do not appear. 

This permanence principle shall be applied as follows: A(e) is the state- 
ment: "Takeeq. (10.18) as initial condition for eqs. (10.13), (10.14), (10.15) 
and calculate the corresponding metric (10.10), (10.16). Then there exists 
an r < / such that RijR lj > k 2 at r = r ." Remark: At this point it is not 
essential whether r is a standard or an infinitely large non-standard number. 
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Proof of A(e) for infinitesimals e with e 7^ 0: the difference of eq. (10.24) 
in relation to the exact solution is of the order e 2 , i.e. the relative error 
is infinitesimally small. For increasing r, becomes arbitrary large, i.e. 
(1 + j3)~ l becomes small and RijR lj increases to arbitrarily large values, cf. 
eq. (10.24). Now take r such that with metric (10.10), (10.24) RijR ij > 2k 2 
holds. Then, for the exact solution, RijR lj > k 2 holds at r = r , because 
their difference was shown to be infinitesimally small. 

Now the permanence principle tells us that there exists a positive standard 
real 5 > such that for all e with < \e\ < 5 the corresponding exact solution 
has a point r such that at r = r the inequality RijR^ > k 2 holds. But 
u RijR % i < k 2 " is another necessary condition for a solution to lie within a 
neighbourhood of Minkowski spacetime. Remark: Supposed this tq is an 
infinitely large non-standard number, then by continuity arguments also a 
standard finite number with the same property exists. 

The statement is proved, but we have learned nothing about the actual 
value of the number 5. Here, numerical calculations may help. They were 
performed as follows: the power series for the functions rji were calculated 
up to the r 6 -term, then these functions taken at x = — 4 i.e. r = 0.018 I 
were used as initial conditions for a Runge-Kutta integration. We got the 
following result: firstly, for e = ±10~ 5 and r < 10 /, the relative difference 
between the linearized solution eq. (10.23) and the numerical one is less than 
2 per thousand. Secondly, for < \e\ < 1 the behaviour u f3 — > — oo • sgne" 
is confirmed. That means, the statement made above keeps valid at least 
for 5 = 1. Remark: For large values e the power series for the rji converge 
very slowly, and therefore other methods would be necessary to decide about 
asymptotical flatness. 
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10.5 Discussion — Einstein's particle programme 

Fourth-order gravitational field equations could be taken as a field theoretical 
model of ordinary matter, the energy-momentum tensor of which is defined 
by 

nT^R^-^R. (10.26) 

For our case one obtains at r = that T*j represents an ideal fluid with the 
equation of state p* = /i*/3 and, cf. eq. (10.25), 

//(0) = hkeK' 1 = 5e/8nGl 2 . (10.27) 

Inserting \e\ > 1 and I < 1.3 • 10~ 13 cm into eq. (10.27), we obtain 

|//(0)| > 1.5- 10 53 g cnT 3 . (10.28) 

Therefore it holds: if there exists a non-trivial static spherically symmetric 
asymptotically flat singularity-free solution of eq. (10.8) at all, then the corre- 
sponding particle would be a very massive one: its phenomenological energy 
density exceeds that of a neutron star by at least 40 orders of magnitude. 

The resulting statement can be understood as follows: For a small cur- 
vature the 4th order corrections to Einstein's equations are small, too, and 
the situation should not be very different from that one we know from Ein- 
stein's theory. Now we want to refer to a problem concerning Schwarzschild 
coordinates: the transition from a general static spherically symmetric line 
element to Schwarzschild coordinates is possible only in the case that the 
function "invariant surface of the sphere r = const, in dependence on its 
invariant radius" has not any stationary point. Here two standpoints are 
possible: either one takes this as a natural condition for a reasonable particle 
model or one allows coordinate singularities in eq. (10.10) like u (3 > — 1, 
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and (3 — — 1 at single points" which require a special care. For (3 < — 1 one 
would obtain a cosmological model of Kantowski-Sachs type. Eqs. (10.5) till 
(10.16) remain unchanged for this case. The discussion made above is not 
influenced by this. 

The statement on the existence of solutions can be strengthened as fol- 
lows: Fiedler and Schimming [78] proved that the solutions are analytical in 
a neighbourhood of r = 0. Further, the differential equation is an analytical 
one and, therefore, in the subspace of singularity-free solutions they remain 
so in the limit r — > oo. Then, there exists only a finite or countably infinite 
set of values e n such that the corresponding solution becomes asymptotically 
flat; the question, whether this set is empty or not, shall be subject of further 
investigation. Furthermore, the e n have no finite accumulation point. That 
means, there exists at most a discrete spectrum of solutions. 

With respect to this fact, we remark the following: as one knows, Ein- 
stein's theory is a covariant one. But besides this symmetry, it is homotheti- 
cally invariant, too. That means, if ds 2 is changed to e 2x ds 2 with constant Xi 
then the tensor Rij — \g%jR remains unchanged, whereas the scalars R and 
C will be divided by e 2x and e 4x respectively. From this it follows: with one 
solution of Einstein's vacuum equation one obtains by homothetic invariance 
just a one-parameter class of solutions. On the other hand, the sum R + l 2 C 
has not such a symmetry and, therefore, one should not expect that a one- 
parameter family of solutions globally exists, and this is just in the scope of 
the particle programme where a definite particle's mass is wanted. 
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10.6 Non-trivial solutions of the Bach equation exist 

Following [204] , we show that solutions of the Bach equation exist which are 
not conformal Einstein spaces. In connection with fourth-order gravitational 
field equations, cf. e.g. Weyl 1921 [268], [37] where the breaking of conformal 
invariance was discussed, the original Bach equation, 

Bij = 0, (10.29) 

enjoys current interest, see e.g. [64], [71] and [91]. Eq. (10.29) stems from a 
Lagrangian 

L=±y/=ijC ijkl C ijkl , (10.30) 
and variation gives, cf. Bach 1921 [8], 

-j= 8L/5g* = B i3 = 2 C% \ ba + C% b R ba . (10.31) 

An Einstein space, defined by 

Rij^Xgij, (10.32) 

is always a solution of the Bach equation (10.29). But eq. (10.29) is confor- 
mally invariant, and therefore, each metric, which is conformally related to 
an Einstein space, fulfils eq. (10.29), too. We call such solutions trivial ones. 

Now the question arises whether non-trivial solutions of the Bach equation 
eq. (10.29) do or do not exist. In this section we will give an affirmative 
answer. As a by-product, some conditions will be given under which only 
trivial solutions exist. Observe that eq. (10.29) is conformally invariant 
whereas eq. (10.32) is not. Therefore, a simple counting of degrees of freedom 
does not suffice. 
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Because the full set of solutions of eq. (10.29) is not easy to describe, let 
us consider some homogeneous cosmological models. Of course, we have to 
consider anisotropic ones, because all Friedmann models are trivial solutions 
of eq. (10.29). Here, we concentrate on the diagonal Bianchi type I models 



with Hubble parameters hi = a^ddi/dt, h = S/i, and anisotropy parameters 
mi = K — h/3. The Einstein spaces of this kind are described in [134] eq. 
(11.52), for A = it is just the Kasner metric Oj = t Pi , = Ep? = 1. 
All these solutions have the property that the quotient of two anisotropy 
parameters, mi/rrij, which equals (3p, — l)/(3pj — 1) for the Kasner metric, 
is independent of t, and this property is a conformally invariant one. Fur- 
thermore, it holds: A solution of eqs. (10.29), (10.33) is a trivial one, if and 
only if the quotients rrii/mj are constants. 

Restricting now to axially symmetric Bianchi type I models, i.e., metric 
(10.33) with hi = h 2 , the identity Em 8 = implies m 1 /m 2 = 1, m 3 /mi = 
m 3 /m 2 = —2 , i.e., each axially symmetric Bianchi type I solution of eq. 
(10.29) is conformally related to an Einstein space. Analogously, all static 
spherically symmetric solutions of eq. (10.29) are trivial ones, cf. [77]. 

Finally, the existence of a solution of eqs. (10.29), (10.33) with a non- 
constant mi/m.2 will be shown. For the sake of simplicity we use the gauge 
condition h — 0, which is possible because of the conformal invariance of 
eq. (10.29). Then the 00 component and the 11 component of eq. (10.29) 
are sufficient to determine the unknown functions h\ and h 2 ; h 3 — —hi — h 2 
follows from the gauge condition. Defining r = (h\ + h\h 2 + h 2 2 ) x l 2 and 




(10.33) 
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p = hi/r, eq. (10.29) is equivalent to the system 

3d 2 {pr)/dt 2 = 8pr 3 + 4c, c = const., p 2 < 4/3 , (10.34) 

9(dp/dt) 2 r 4 = [2rd 2 r/dt 2 - (dr/dt) 2 - 4r 4 ](4r 2 - 3p 2 r 2 ) . (10.35) 

As one can see, solutions with a non-constant p exist, i.e., mi/m 2 is not 
constant for this case. 

Result: Each solution of the system eqs. (10.33), (10.34), (10.35) with 
dp/dt 7^ represents a non-trivial solution of the Bach equation (10.29). 
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11 The Newtonian limit of non- linear gravity 

In this chapter, the weak-field slow-motion limit of fourth and higher-order 
gravity will be deduced. Here we follow [208] for the fourth-order case and 
[187] for the general case. More explicitly: we consider the Newtonian limit 
of the theory based on the Lagrangian 

C= [R+it a k Ra k R) y/=jj. (11.1) 

The gravitational potential of a point mass turns out to be a combination 
of Newtonian and Yukawa terms. For fourth-order and sixth-order gravity, 
p = and p = 1 respectively, the coefficients are calculated explicitly. For 
general p one gets the potential to be 

$ = m/r (l + a exp(-r//^ (11.2) 

with certain coefficients q fulfilling the relation 

Eq = 1/3. (11.3) 

i=0 

Therefore, the potential is always unbounded near the origin, see also [10]. 

11.1 The Newtonian limit of 4th-order gravity 

Let us consider the gravitational theory defined by the Lagrangian 

L g = (8ttG)- 1 (r/2 + (aRijR ij + (3R 2 )l 2 ) . (11.4) 

G is Newton's constant, I a coupling length and a and j3 numerical param- 
eters. Ry and R are the Ricci tensor and its trace. Introducing the matter 
Lagrangian L m and varying L g + L m one obtains the field equation 

Eij + aHij + fiGij = 8nGT ij . (11.5) 
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For a = (3 = this reduces to General Relativity Theory. The explicit 
expressions and Gy can be found in sections 8.4 and 10.2, see also Stelle 
[249]. 

In a well-defined sense, the weak-field slow-motion limit of Einstein's the- 
ory is just Newton's theory, cf. Dautcourt 1964 [60]. In the following, we 
consider the analogous problem for fourth-order gravity eqs. (11.4), (11.5). 
In some cases, the Newtonian limit of the theory defined by eq. (11.5) has 
already been deduced in the literature: For the special case a = see Pech- 
laner and Sexl [180] or Polijevktov-Nikoladze [183]. For the case a + 2(3 = 
see Havas [105] or Jankiewicz [121], and for a + 3(3 = see Borzeszkowski, 
Treder and Yourgrau [37]. Cf. also Anandan [5], where torsion has been 
taken into account. 

The slow- motion limit can be equivalently described as the limit c — > oo, 
where c is the velocity of light. In this sense we have to take the limit G — > 
while G ■ c and I remain constants. Then the energy-momentum tensor Ty 
reduces to the rest mass density p: 

Tij = $qp, (11.6) 

x° = t being the time coordinate. The metric can be written as 

ds 2 = (1 - 2<f))dt 2 - (1 + 2i()(dx 2 + dy 2 + dz 2 ) . (11.7) 

Now eqs. (11.6) and (11.7) will be inserted into eq. (11.5). In our approach, 
products and time derivatives of <fi and ip can be neglected, i.e., 

i? = 4A^-2A0, where A/ /.,,,. • /, /;; • /.,,. (11.8) 

Further R 00 = -A0, H 00 = -2AR 00 - AR and G o = -4Ai?, where 1 = 1. 
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Then it holds: The validity of the 00-component and of the trace of eq. 
(11.5), 

R 00 - R/2 + aH 00 + (3G 00 = 8nGp (11.9) 

and 

-R-4(a + 3l3)AR = 8nGp, (11.10) 
imply the validity of the full eq. (11.5). 

Now, let us discuss eqs. (11.9) and (11.10) in more details: Eq. (11.9) 
reads 

-A0 - R/2 + a(2AA0 - AR) - 4{3AR = 8nGp . (11.11) 

Subtracting one half of eq. (11.10) yields 

-A0 + 2aAA0+ (a + 2(3)AR = A^Gp. (U-12) 

For a + 2(3 = one obtains 

-(1 -2aA)A0 = AnGp (11.13) 

and then ip = <fi is a solution of eqs. (11.9), (11.10). For all other cases the 
equations for and ip do not decouple immediately, but, to get equations 
comparable with Poisson's equation we apply A to eq. (11.10) and continue 
as follows. 

For a + 3(3 = one gets from eq. (11.12) 

-(1 - 2«A)A0 = 47rG(l + 2aA/3)p . (H-14) 

The A-operator applied to the source term in eq. (11.14) is only due to the 
application of A to the trace, the original equations (11.9), (11.10) contain 
only p itself. 
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For a = one obtains similarly the equation 

-(1 + 12/3A)A0 = 47rG(l + 16/3A)p . (11.15) 

For all other cases - just the cases not yet covered by the literature - the 
elimination of tp from the system (11.9), (11.10) gives rise to a sixth-order 
equation 

-(1 + 4(a + 3/3) A) (1 - 2aA)A</> = 47rG(l + 2(3a + 8/3)A)p . (11.16) 

Fourth-order gravity is motivated by quantum-gravity considerations and 
therefore, its long-range behaviour should be the same as in Newton's theory. 
Therefore, the signs of the parameters a, (3 should be chosen to guarantee an 
exponentially vanishing and not an oscillating behaviour of the fourth-order 
terms: 

a>0, a + 3/3<0. (U-17) 

On the other hand, comparing parts of eq. (11.16) with the Proca equation 
it makes sense to define the masses 

m 2 = [2aY l/2 and m = (-4(a + 3/3))" 1/2 . (11.18) 

Then eq. (11 .17) prevents the masses of the spin 2 and spin gravitons to 
become imaginary. 

Now, inserting a delta source p = m5 into eq. (11.16) one obtains for cf> 
the same result as Stelle [249] , 



4> = mr x (l + exp(— m r)/3 — 4exp(— m 2 r)/3) . (11.19) 
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To obtain the metric completely one has also to calculate ip. It reads 

ip = mr _1 (l — exp(— m r)/3 — 2exp(— m 2 r)/3) . (11.20) 

For finite values m and m 2 these are both bounded functions, also for r — > 0. 
In the limits a — > 0, i.e. m 2 — > oo, and ct + 3/3 — > 0, i.e. mo — > oo, the terms 
with m and m 2 in eqs. (11.19) and (11.20) simply vanish. For these cases 
and ip become unbounded as r — > 0. 

Inserting eqs. (11.19) and (11.20) into the metric (11.7), the behaviour 
of the geodesies shall be studied. First, for an estimation of the sign of the 
gravitational force we take a test particle at rest and look whether it starts 
falling towards the centre or not. The result is: for m < 2m 2 , gravitation 
is always attractive, and for m > 2m 2 it is attractive for large but repelling 
for small distances. The intermediate case m = 2m 2 , i.e., 3a + 8/3 = 0, is 
already known to be a special one from eq. (11.16). 

Next, let us study the perihelion advance for distorted circle-like orbits. 
Besides the general relativistic perihelion advance, which vanishes in the 
Newtonian limit, we have an additional one of the following behaviour: for 
r — > and r — > oo it vanishes and for r l/m and r ~ l/m 2 it has local 
maxima, i.e., resonances. 

Finally, it should be mentioned that the gravitational field of an extended 
body can be obtained by integrating eqs. (11.19) and (11.20). For a spheri- 
cally symmetric body the far field is also of the type 

mr _1 ^l + a exp(— m r) + 6exp(— m 2 r)^ , (11.21) 

and the factors a and b carry information about the mass distribution inside 
the body. 
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11.2 Introduction to higher-order gravity 

One-loop quantum corrections to the Einstein equation can be described by 
curvature-squared terms and lead to fourth-order gravitational field equa- 
tions; their Newtonian limit is described by a potential "Newton + one 
Yukawa term", cf. e.g. Stelle [249] and Teyssandier [253]. A Yukawa poten- 
tial has the form exp(— r/l)/r and was originally used by Yukawa [274] to 
describe the meson field. 

Higher-loop quantum corrections to the Einstein equation are expected 
to contain terms of the type Rn k R in the Lagrangian, which leads to a grav- 
itational field equation of order 2k + 4, cf. [97]. Some preliminary results to 
this type of equations are already due to Buchdahl [42]. For k — 1, the cos- 
mological consequences of the corresponding sixth-order field equations are 
discussed by Berkin and Maeda [25], and by Gottlober, Miiller and Schmidt 
[94]. 

In the present chapter we deduce the Newtonian limit following from this 
higher order field equation. The Newtonian limit of General Relativity The- 
ory is the usual Newtonian theory, cf. e.g. Dautcourt 1964 [60] or Stephani 
[250]. From the general structure of the linearized higher-order field equa- 
tion, cf. [218], one can expect that for this higher-order gravity the far field 
of the point mass in the Newtonian limit is the Newtonian potential plus a 
sum of different Yukawa terms. And just this form is that one discussed in 
connection with the fifth force, cf. [85], [242] and [197]. Here we are inter- 
ested in the details of this connection between higher-order gravity and the 
lengths and coefficients in the corresponding Yukawa terms. 
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11.3 Lagrangian and field equation 

Let us start with the Lagrangian 

C = \ R + 12 akRO k R) ■ y/=g, a p ^0. (11.22) 
V fc=0 / 

In our considerations we will assume that for the gravitational constant G and 
for the speed of light c it holds G = c = 1. This only means a special choice of 
units. In eq. (11.22), R denotes the curvature scalar, □ the D'Alembertian, 
and g the determinant of the metric. Consequently, the coefficient a k has the 
dimension "length to the power 2k + 2" . 

The starting point for the deduction of the field equation is the principle 
of minimal action. A necessary condition for it is the stationarity of the 
action: 

' 8ti V" v '~. (11.23) 



5gij 

where T y denotes the energy-momentum tensor. The explicit equations for 

P l] V^9 = -^- (11-24) 

ogij 

are given in [218]. Here we only need the linearized field equation. It reads, 
cf. [97] 

pij = Rij _R g ij + 2 J2 a k [g ij n k+1 R - U k R^] = 8nT i] , (11.25) 

2 k=0 

and for the trace it holds: 

n - 1 p 

9H ■ P ij = 7) — R + 2n^ a k [g ij D k+1 R] = 8nT . (11.26) 

n is the number of spatial dimensions; the most important application is of 
cause n — 3. From now on we put n — 3. 
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11 A The Newtonian limit in higher-order gravity 

The Newtonian limit is the weak-field static approximation. So we use the 
linearized field equation and insert a static metric and an energy-momentum 
tensor 

T l3 = 5°5°p, p>0 (11.27) 

into eq. (11.25). 

Without proof we mention that the metric can be brought into spatially 
conformally flat form, and so we may use 

9ij Vij fij i 

r/ij = diag(l, —1, —1, —1) and 
fij = diag(-2$, -2*, -2*, -2*) . (11.28) 

Then the metric equals 

ds 2 = (1 - 2<$>)dt 2 - (1 + 2^)(dx 2 + dy 2 + dz 2 ) , (11.29) 

where <£> and \1/ depend on x, y and z. Linearization means that the metric 
gij has only a small difference to rjij; quadratic expressions in fj and its 
partial derivatives are neglected. We especially consider the case of a point 
mass. In this case it holds: $ = $(r), * = \I/(r), with 

r = ^x 1 + y 2 + z 2 , (11.30) 

because of spherical symmetry and p = m5. Using these properties, we 
deduce the field equation and discuss the existence of solutions of the above 
mentioned type. 

At first we make some helpful general considerations: The functions $ 
and \1/ are determined by eq. (11.25) for % — j — and the trace of eq. 



11.4 Newtonian limit in higher-order gravity 



167 



(11.25). If these two equations hold, then all other component-equations are 
automatically satisfied. For the 00-equation we need R 00 : 

i?oo = -A$. (11.31) 

Here, the Laplacian is given is as usual by 

d 2 ^ (P_ 

dx 2 dy 2 dz 2 

For the inverse metric we get 

g ij = diag (1/(1 - 2$), -1/(1 + 2*), -1/(1 + 2*), -1/(1 + 2*)) (11.33) 

and 1/(1 — 2$) = 1 + 2$ + h(&), where h(&) is quadratic in <3> and vanishes 
after linearization. So we get 

g ij = rf' - f j . (11.34) 

In our coordinate system, / u equals for all For the curvature scalar 
we get 

R = 2(2A^ - A$) . (11.35) 

Moreover, we need expressions of the type O k R. OR is defined by OR = 
R,ij9 % \ where ";" denotes the covariant derivative. Remarks: Because of 
linearization we may replace the covariant derivative with the partial one. 
So we get 

U k R= ( i) fc 2(-A fc+1 $ + 2A fc+1 ^) (11.36) 
and after some calculus 

-8vrp = A$ + A^ . (11.37) 

We use eq. (11.37) to eliminate \1/ from the system. So we get an equation 
relating $ and p = mS. 

(p \ p 

P + a k (-l) k A k+1 p = A$ + 6 J2 a k (-l) k A k+2 <S> . (11.38) 
fc=0 / k=0 
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In spherical coordinates it holds 

A$ = -$ >r + $ >rr , (11.39) 
because <£> depends on the radial coordinate r only. 

We apply the following lemma: In the sense of distributions it holds 



A fie-/') =^-^-4^. (11.40) 



Now we are ready to solve the whole problem. We assume 

$ = - (l + I> exp(-rA)) , Zi > . (11.41) 



r 



=o / 

Without loss of generality we may assume U ^ lj for i ^ j. Then eq. (11.38) 
together with that lemma yield 

+24.EtE^(-l) J+1 iA^ (11.42) 

k=0j=ki=0 *i 

where ti = If ; therefore also t« ^ tj for i ^ j. This equation is equivalent 
to the system 

i 

E 

i=0 



5> = l/3, (11.43) 

ES = °; s = l,...p (11.44) 

i=0 Z i 

C 1 + 6 E a fe (-l) fc r fe = , 1 = 0,...?. (11.45) 



q c, 



k=0 



11.4 Newtonian limit in higher-order gravity 



169 



From eq. (11.45) we see that the values ti represent q + 1 different solutions 
of one polynomial. This polynominal has the degree p + 1. Therefore q < p. 



Now we use eqs. (11.43) and (11.44) . They can be written in matrix 
form as 

/ 1/3 \ 



' 1...1 
l/t ...l/t q 



Cl 







(11.46) 



v i/«g . . . v<? y \c q J \ J 

Here, the first q + 1 rows form a regular matrix, the Vandermonde matrix. 
Therefore, we get 

9 



fe=0 



(11.47) 



with certain coefficients \jk i.e., the remaining rows depend on the first q + 1 
ones. If \j 7^ then the system has no solution. So A j0 = for all g + 1 < 
j < p. But for q < p we would get 

9 



fc=i 



11.48) 



and this is a contradiction to the above stated regularity. Therefore p equals 
q. The polynomial in (14) may be written as 



I 1 l/fo.-.l/fg \ 



^ 1 llt 9 ...\/V p j 



( 



a 



\ 



V (-i) p s / 



(11.49) 



This matrix is again a Vandermonde one, i.e., there exists always a unique 
solution (a , . . . a p ), which are the coefficients of the quantum corrections to 
the Einstein equation, such that the Newtonian limit of the corresponding 
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gravitational field equation is a sum of Newtonian and Yukawa potential with 
prescribed lengths /j. A more explicit form of the solution is given in section 
11.6. 



11.5 Discussion of the weak-field limit 



Let us give some special examples of the deduced formulas of the Newtonian 
limit of the theory described by the Lagrangian (11.22). If all the vanish 
we get of course the usual Newton theory 



in 

$ = — , A$ = -Aty5 . 
r 



11.50) 



<E> and \I> refer to the metric according to eq. (11.29). For p = we get for 

a < 



$ = 



m 



1 

1 + 



-r/yJ—Qao 



cf. [249] and 



1 p-r/V-oao 

3 

cf. [208]. For a > no Newtonian limit exists. 



m 

* = — 

r 



For p — 1, i.e., the theory following from sixth-order gravity 
we get, see [187] 



and 



where 



m 

$ = _ 

r 



* = — 

r 



1 + c e- r/ '° + Cl e- r/11 



r/ii 



1 - c e~ r/l0 - Cl e- r/h 



a 



(11.51) 
(11.52) 



(11.53) 

(11.54) 
(11.55) 
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and 



l ,i = \/-3a ± ^9a 2 + 6 Gl . (11.57) 



This result is similar in structure but has different coefficients as in fourth- 
order gravity with included square of the Weyl tensor in the Lagrangian. 

The Newtonian limit for the degenerated case l = h can be obtained by 
a limiting procedure as follows: As we already know a < 0, we can choose 
the length unit such that a = —1/3. The limiting case 9a^ + 6a x — > may 
be expressed by a\ = — 1/6 + e 2 . After linearization in e we get: 

h = 1± \[z/2sci = 1/6 ± l/(6"\/6e) (11.58) 

and applying the limit e — > to the corresponding fields $ and \I/ we get 

$ = m /r{l + (1/3 + r/6)e~ r } 
* = m/r{l - (l/3 + r/6)e~ r }. (11.59) 

For the general case p > 1, the potential is a complicated expression, but 
some properties are explicitly known, these hold also for p = 0,1. One gets 

$ = m/r ^1 + £j Cj exp(-r/Z^ (11.60) 

and 

* = m/r ^1 - J2 Qexp(-r/Z^ (11.61) 

where J2 c i — 1/3; S means Y%=o an d ^ an d Q are > U P to permutation of 
indices, uniquely determined by the Lagrangian. 

There exist some inequalities between the coefficients a iy which must be 
fulfilled in order to get a physically acceptable Newtonian limit. By this 
phrase we mean that besides the above conditions, additionally the fields <3> 
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and \& vanish for r — > oo and that the derivatives dQ/dr and d^/dr behave 
like 0(l/r 2 ). These inequalities express essentially the fact that the U are 
real, positive, and different from each other. The last of these three conditions 
can be weakened by allowing the q to be polynomials in r instead of being 
constants, cf. the example with p — 1 calculated above. 

The equality £)cj — 1/3 means that the gravitational potential is un- 
bounded and behaves, up to a factor 4/3, like the Newtonian potential for 
r pa 0. The equation <3> + \P = 2m /r enables us to rewrite the metric as 

ds 2 = (1 - 29) [(1 - 2m/r)dt 2 - (1 + 2m/r)(dx 2 + dy 2 + dz 2 )] , (11.62) 

which is the conformally transformed linearized Schwarzschild metric with 
the conformal factor 1 — 29, where 

0=™5>e~^ (11.63) 

can be expressed as functional of the curvature scalar, this is the linearized 
version of the conformal transformation theorem, cf. [218]. For an arbitrary 
matter configuration the gravitational potential can be obtained by the usual 
integration procedure. 



11.6 A homogeneous sphere 



For general p and characteristic lengths U fulfilling < 1$ < h < . . 
write the Lagrangian as 



. < l p we 



R 



R r 



ill 



l 2 p )R+(l 2 l 2 1+ l 2 l 2 2 ... + l 2 l 2 p )nR+ 



(l 2 lil 2 + ... + l^lUl^R • ... • /;•/-...• l 2 ^R (11.64) 
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the coefficients in front of O l R in this formula read 

E fit- 

0<jo<ji<---<ji<P "i=0 



(11.65) 



Using this form of the Lagrangian, the gravitational potential of a point mass 
reads 



m 

$ = — 

r 



m 

* = — 

r 



i+^B-i^ml-ir 1 ^ 

' i=0 j# ' 



i-i^-i^nil-ir 1 ^ 



i=Q 



11 I 12 



(11.66) 
(11.67) 



For a homogeneous sphere of radius r and mass m we get 



m 

$ = _ 

r 



where 



1 + 4 E e_r/ ' 1 £ 5 *( r o cosh(r /Zi) - k sinh^o//,)) , (11.68) 
r o i=o 

Ci = (-l) <+1 IIlI- i r 1 - (11-69) 
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12 Cosmic strings and gravitational waves 

We consider strings with the Nambu action as extremal surfaces in a given 
spacetime, thus, we ignore their back reaction. Especially, we look for strings 
sharing one symmetry with the underlying spacetime. If this is a non-null 
symmetry the problem of determining the motion of the string can be dimen- 
sionally reduced. Following [232], we get exact solutions for the following 
cases: straight and circle-like strings in a Friedmann background, straight 
strings in an anisotropic Kasner background, different types of strings in the 
metric of a gravitational wave. The solutions will be discussed. 

12.1 Introduction to cosmic strings 

To give detailed arguments for considering strings means carrying coals to 
Newcastle. So we only list the main points: A string is, generally speaking, 
an object possessing a two-dimensional world surface, which is also called 
a world sheet, in contrast to a point particle possessing a one-dimensional 
world line, cf. the review article by Vilenkin [261]. In details we have 

1. One considers strings with the Nambu action in a D-dimensional 
flat spacetime. The theory can be consistently quantized for D = 26 only: 
otherwise the light cone quantization leads to a breakdown of Lorentz co- 
variance, cf. Green, Schwarz and Witten [98]. But we consider a classical 
non-quantized theory only and require D = 4 henceforth. 

2. Cosmic strings are one-dimensional topological defects in gauge field 
theories. They have a large mass and can be seeds for larger objects by 
accretion processes, cf. Zeldovich [276]. 

3. One looks for solutions of the Einstein equation with distribution- 



176 



12 COSMIC STRINGS AND GRAVITATIONAL WAVES 



valued energy-momentum tensor whose support is a two-dimensional sub- 
manifold of indefinite signature. The equation of state is p z = —/i and the 
solution is called cosmic string. These are good candidates for the origin of 
structure, especially for the formation of galaxies, in the early universe. 

4. If one looks for strings according to 2. or 3. in a given spacetime, 
i.e., with negligible back reaction of the string onto spacetime geometry, one 
arrives at the Nambu action, too, see Nielsen and Olesen [171] for 2. and 
Geroch and Traschen [87] for 3. In other words, cosmic strings are extremal, 
maximal or minimal in dependence of the boundary conditions, surfaces of 
indefinite signature in a given spacetime. This approach is justified if the 
diameter of the string is small compared with the curvature radius of the 
underlying manifold and if for the string tension a' <C l/G holds; we use 
units with c = 1. Normally, one thinks in orders of magnitude Go 1 = 10~ 6±2 , 
see Brandenberger [39]. In other words, the string is supposed to possess a 
mass per unit length of 10 22±2 g/cm, if the phase transition is supposed to 
be at the GUT-scale, see Frolov and Serebriany [80]. 

In [81], a stationary string in the Kerr- Newman metric has been consid- 
ered. We use the method developed there and apply it to other cases. In 
this chapter we shall adopt the 4. approach and try to give some geometric 
insights into the motion of a string. To this end we give a sample of closed- 
form solutions for a special class of string solutions: strings which share a 
spacelike or timelike isometry with the underlying background metric. 

12.2 The main formulae 

The string is a two-dimensional extremal surface of indefinite signature. Let 
us take coordinates (r, a) = (y A ), A = 0, 1 within the string and coordinates 
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(x l ), % — 0, 1, 2, 3 for the spacetime V4 with the metric g^. The string is 
given by specifying the four functions x l {y A ). The signature for the metric 
gij is (H ). The tangents to the string are 

e\ = dx i /dy A (12.1) 

and induce the metric 

h A B = e\e j B gij (12.2) 

at the string world sheet. The signature of the string is required to be (H — ), 
thus 

h = deth AB <0. (12.3) 

The action to be varied is 

1 = li 1 1 ^~ hi ° iT ■ < 12 ' 4 > 

Instead of writing down the full equations we specialize to the case we are in- 
terested here: we require that an isometry for both the underlying spacetime 
and the string exists. Let hi be a non-null hypersurface-orthogonal Killing 
vector field, i.e., 

%% ;fc] = 0, % ; ,) = 0. (12.5) 
Then there exist coordinates x l (x° = t) such 

ds 2 — 9oodt 2 + g 0l pdx a dx 13 . (12.6) 

a, P — 1, 2, 3, the do not depend on t, fcj = d/dt. The sign of g 00 is 
determined by the condition g 0Q hik 1 > 0: for timelike ki we have goo > 
and for spacelike ki we have g o < 0. The requirement that h L is also an 
isometry of the string gives us the possibility to specify the functions x l (y ) 
to be t — t, x a depends on a only. Then we get 

e\ = (1, 0, 0, 0) , e\ = (0, dx a /da) (12.7) 
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and h { 



■oi — 



0, hoo — goo, 



hn = g a (3 dx a /da dx 13 /da , 



I 




—goo g a /3 dx a /da dx^/da dadr 



(12.8) 



The integrand does not depend on r, so we can omit the r integration. 
Therefore, extremizing eq. (12.8) is the same as solving the geodesic equation 
for the auxiliary metric f a p of a V 3 defined by 



Remarks: 1. If ki is not hypersurface-orthogonal then we get contains terms 
with go a , and one has to add go a go/3 to the r.h.s. of the definition of f a p. 

2. The dimensional reduction is not fully trivial: the geodesic equation for 
metric f a p corresponds to compare the strings in eq. (12.8) with other strings 
sharing the same isometry induced by ki, whereas the Nambu action has to 
be compared with all other strings, too. But writing down all full equations 
or counting the degrees of freedom one can see that at our circumstances no 
difference appears. 

3. For a timelike ki, f a p is positive definite, and the condition h < 
is automatically fulfilled. On the other hand, for a spacelike ki, f a p is of 

signature (H ), and the validity of h < requires the vector dx a /da to 

be timelike, i.e., the root in eq. (12.8) has to be real. 

12.3 The string in a Friedmann model 

Now we specify the underlying V 4 to be a spatially flat Friedmann model 



fa/3 — —gOOga/3 ■ 



(12.9) 



ds 2 = dt 2 - a 2 {t){dx 2 + dy 2 + dz 2 ) . 



(12.10) 
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12.3.1 The open string 

First we use the spacelike Killing vector d/dx of metric (12.10). We have 
in mind an infinitely long straight string moving through the expanding uni- 
verse. We apply the formalism deduced in the previous section and get the 
following result: we write metric (12.10) in the form 

ds 2 = -a 2 dx 2 + dt 2 - a 2 (dy 2 + dz 2 ) . (12.11) 

The metric f a p then reads 

ds 2 {3) = a 2 dt 2 - a 4 (dy 2 + dz 2 ) . (12.12) 

Without loss of generality the string is situated at z = and moves into the 
^/-direction according to 

ds/ 3 N = a 2 dt 2 — a 4 dy 2 

1 da . 9 da , , 
t + — -t 2 + 2— ay 2 = with • = d/d\ 
a dt dt 

a 4 y = M 

a 2 i 2 - a 4 y 2 = 1 
r dt 

y(t) = M / (1213) 

J a{t)^a\t) + M 2 

where M is an integration constant. The natural distance of the string from 
the origin is 

s(t) = a(t) ■ y(t) . (12.14) 

I. Example. Let a(t) = t n , n > then for t 3> 1, cf. Stein-Schabes and 
Burd [247], 

y(t) w Mt 1 " 3 ", s(t) w Mt x - 2n for n ^ 1/3 (12.15) 
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and y(t) ~ Mint for n = 1/3. Interpretation: n — 0, i.e., the absence of 
gravity, implies a linear motion as it must be the case. We have 

lim y(t) = oo (12.16) 

t — >oo 

if and only if n < 1/3, i.e., for n > 1/3 the string comes to rest with respect 
to the cosmic background after a finite time. A more stringent condition to 
be discussed is that the string comes to rest in a natural frame of a suitably 
chosen reference galaxy. This means 

lim s(t) < oo (12.17) 

t^oo 

and is fulfilled for n > 1/2. Therefore, the most interesting cases n = 1/2, 
the radiation model and n = 2/3, the Einstein-de Sitter dust model, have 
the property that straight open strings come to rest after a sufficiently long 
time independently of the initial conditions. 

2. Example. Let a(t) = e Ht , the inflationary model, H > 0. Then 
y(t) = e~ 3Ht , s(t) = e~ 2Ht . We have the same result as in the first example 
with n ^> 1. 

3. Example. Let 

a(t) = -t 2/3 (l + r 2 cos(mt)) . (12.18) 

This background metric is from damped oscillations of a massive scalar field 
or, equivalently, from fourth-order gravity L = R + m~ 2 R 2 . We get with eq. 
(12.13) 

y(t) pa — l/t — ct~ A sinmt , c = const, for t — > oo , (12.19) 

which is only a minor modification of the result of the first example with 
n — 2/3, therefore, one should not expect a kind of resonance effect between 
the open string and a massive scalar field. 
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12.3.2 The closed string 

We insert dy 2 + dz 2 = dr 2 + r 2 rf$ 2 into metric (12.10) and use the spacelike 
Killing vector <9/<9$: its trajectories are circles, so we have in mind a closed 
string with radius r moving and oscillating in the expanding universe. The 
corresponding geodesic equation leads to 

A(_ a Vr) - a 2 ri 2 + a 4 r(x 2 + r 2 ) = , (12.20) 
dX 

^(aV±) = 0, (12.21) 

i 2 -a 2 {x 2 + f 2 ) = —, (12.22) 
a 2 r 2 

where the dot denotes d/dX, X is the natural parameter along the geodesic. 
We are mainly interested in solutions not moving into the x-direction, to 
understand the oscillating behaviour. Inserting x = into this system of 3 
equations and using the fact that always t ^ holds, this system can be 
brought into the form 

r • rV - 2rr'a 4 a - r'V + 3rr'a 2 a + a = , (12.23) 

where the dash denotes d/dt. 

4. Example. Let a = 1, i.e., we have the flat Minkowski spacetime. Then 
this equation reduces to 

rr" = r' 2 -l. (12.24) 

The solution reads 

r(t)=r cos((f-to)/r ), r' = - sin((t - t )/r ) . (12.25) 

This is the oscillating solution for closed strings. At points t where r = 0, 
we have \r'\ = 1. These points are the often discussed cusp points of the 
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string, where the interior string metric becomes singular and the motion 
approximates the velocity of light, cf. Thompson [254]. 

Let us compare eq. (12.24) with the analogous equation for a positive 
definite back-ground metric. It reads 

rr" = r' 2 + l (12.26) 

and has solutions with cosh instead of cos. This is the usual minimal 
surface taken up e.g. by a soap-bubble spanned between two circles. 

5. Example. Let a(t) = t n , then we have to solve 

rr"t 2n - 2nrr'h 4n - 1 - r'H 2n + Snrr't 2 "" 1 + 1 = 0. (12.27) 

This equation governs the radial motion of a circle-like closed string in a 
Friedmann background. The solutions can be hardly obtained by analytic 
methods. 

12.4 The string in an anisotropic Kasner background 

Now we take as background metric 

ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)dy 2 - c 2 (t)dz 2 (12.28) 

and d/dx as Killing vector. Astonishingly, the anisotropy has only a minor 
influence on the motion of the string, so we get mainly the same formulae as 
in the isotropic case: the geodesic equations are 

^L(a 2 b 2 y)=0 hence a 2 b 2 y = M 1 (12.29) 
dX 

jL(a 2 c 2 z) = hence a 2 c 2 z = M 2 (12.30) 
d\ 

a 2 i 2 - a 2 c 2 z 2 = 1 (12.31) 
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and can be integrated to yield for MiM 2 ^ 

, , f fa 2 b 2 Ml b 2 \' 1/2 dt „ nnn . 

M = l {-W + Wf7 2+1 ) WV (12 ' 32) 

The equation for z(t) can be obtained from this equation by interchanging 
b <->• c and M x <-> M 2 . 

For the Kasner metric we have a = t p , b = t q , c = t r , p + q + r = 
p 2 + q 2 + r 2 = 1. We get 

y( r )= I / , ^ , =k • ( 12 -33) 

The behaviour for t — > oo in dependence of the values p, g, r can be seen 
from this equation. 



12.5 The string in a gravitational wave 

As background metric we use the plane-wave ansatz 

ds 2 = 2dudv + p 2 (u)dy 2 + q 2 (u)dz 2 . (12.34) 

This metric represents a solution of Einstein's vacuum equation if 

qd 2 p/du 2 + pd 2 q/du 2 = 0. (12.35) 

Let us take d/dy as Killing vector. Then the geodesic equations for the 
auxiliary metric are 

^-{p 2 u) = hence p 2 u = M 1 ^ , (12.36) 
dX 

-^L(p 2 q 2 z) = hence p 2 q 2 z = M 2 , (12.37) 
dX 

-2p 2 iiv - p 2 q 2 z 2 = 1 . (12.38) 
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If we take u as new independent variable, which is possible because of ii ^ 0, 
we get the solutions 

M2 f du' 
Z[U) ~ Mj q\u') ' 

As an example let us take 

p{u) = sin it, q(u) = sinh-u (12.40) 

which produces a vacuum solution. With this choice of p and q we get 

,z(w) = Mcothw, v(u) = -M 2 (4cothw + sin2w-2w)/8. (12.41) 

Another Killing vector of the plane wave is d/dv, but it is a null Killing 
vector and so the reduction used above does not work. But there exists a 
further non-null Killing vector of this plane wave: It reads 

ki = (0,y,H,0), where H = -Jp- 2 du. (12.42) 

By a coordinate transformation 

du 



u = t y = w ■ e G G{u) = J 



p 2 H 



z = k v = x + ^-e- 2G (w 2 -l) (12.43) 

2x2 



we get the form 



— 2G 

ds 2 = e~ 2G p 2 dw 2 + 2dxdt + ^r-^dt 2 + q 2 dk 2 (12.44) 

p 2 H z 



and have to solve the equations 



-pYe-^k 2 - -^e~* a i 2 - 2pV 2G « = 1 . (12.45) 
H 
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For the special case k = one has finally 

^4/- 2C (w + w)* 

Let us suppose a spatially flat Friedmann model eq. (12.10) with scale 
factor a(t) = t n , and n — 2/3, representing the Einstein-de Sitter dust model, 
or n — 1/2, representing the radiation model. There we consider an open 
string which is only a little bit curved such that the approximation of a 
straight string is applicable. At time t = to > we can prescribe place and 
initial velocity vq of the string arbitrarily and get for t — > oo the behaviour 

s(t) « const + M(v )t 1 - 2n (12.47) 

where s denotes the natural distance from the origin. That means, for the 
cases n = 2/3 and = 1/2 we are interested in, the string comes to rest for 
large values t at a finite distance from the origin. 

We compare this result with the analogous motion of a point-particle in 
the same background: in the same approximation we get 

s(t) « const + M(v )t l - n , (12.48) 

i.e., \s\ — > oo as t — > oo, a totally other type of motion. It should be 
mentioned that in the absence of gravity, i.e., n — 0, both motions are of the 
same type, but otherwise not. 

The remaining calculations above indicate that the interaction of the 
motion of the string with scalar field oscillations, with anisotropy, and with 
gravitational waves is quite weak, we did not find any type of resonance 
effects. 



(12.46) 
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